ANGELES PUBLIC LIBRARY 
MAR 2 8 1960 


JOURNAL 


A PUBLICATION OF THE AMERICAN ROCKET SOCIETY 
VOLUME 30 NUMBER 3 MARCH 1960 


Becent Advances in Chemical Kinetics of Homogeneous R 
Air . 


Spectroscopic ausdinces Reverse Jet Flame Stabilization . 


Relativistic Theory of Rocket Flight With Advanced Propulsion Systems. 


Ascent From Inner Circular to Outer Co-Planar Elliptic Orbits. . . . L.Rider 254 
Structural Damage fee Other Effects of Solar Plasmas. 


TECHNICAL NOTES 


Surface Recombination Plow of a Dissociated Diatomic Gas Past 
a Catalytic Flat Plate... .. ,Paul M. Chung and Aemer D. Anderson 262 


Vapor-Liquid Equilibria of Ozone-Oxygen Systems . . . 
Charles K. Hersh, Robert I. Brabets, Gerald M. Platz, Raymond J. Swebla and Donald P. Kirsh 264 


Measurements of the Heat Transfer Coefficients for —, Flowing in a Heated Tube . . 
. M. Fowler and C. F. Warner 266 


A Method for Heat From 
, Roger William E. ‘Smith 268 


_William T. Snyder 270 

Jehm L. Webster and David N. Schultz 

Some Characteristics of the Planar Satellite Orbit... . Josef S. Pistiner 275 
Syrmimetry of the Earth’s Figure... .....- _ J. L. Brenner, R. Fulton and N. Sherman 278 


Thermal Decomposition of Liquid Phase a Nitrate and Its Compatibility With Various 
Metals and Other Additives... . . . Adolph B. Amster and Joseph B. Levy 280 


Refraction Correction boa cute im Earth Satellite ae Tracking for Precision Prediction 


Iterative Solution of the N-Body Problem for Real Time AL. ML. Rauch 264 
Mass Considerations in Ring Supported Solar Sails. E. 287 
Prandtl-Meyer Expansion in Equilibrium Air, Kennet 
Convection in o Fluid at § 


New Patents. 
Book Reviews. . 
Technical Literature. Digest. 


CI 
E& TECHNOLOGY, 
4 q 
Walter G. Zinman 233 
= 
a 


” 
® 
— 
© 
” 
© 
° 
2 
> 
= 
— 
= 
> 
2 
” 
” 
— 
© 
2 
® 
£ 
= 
° 
2 
Ww 
2 
G 
> 
© 
” 
Q 
® 
Oo 
© 
” 
— 
3 
O 
® 
= 
= 
= 
= 


ngineer, there is 
more stimulating than to work ina creative 
environment. Space engineering programs 
now in process at Martin-Denver demand 
unusual creativity and may be your ticket 
to the personal and professional achieve- 
ments which you are seeking. Make your 
desires and qualifications known to N. M. 
Pagan, Director of Technical and Scientific 
Staffing, The Martin Company, (Dept. 4) 
P. O. Box 179, Denver 1, Colorado. 


DENVER DIVISION 


— 
» = 5 
» 


Sometime in the 60’s, America may have a 
manned station in space —and it could look a 
lot like this scale model by Lockheed. 


We have much to learn before it can be built. 
Research is our stairway to space—basic 
research that seeks to discover the new rather 
than develop the known. We cannot predict what 
such research will discover, or when —but we 
delay it or curtail it at our peril. 


Today, at Lockheed’s Missiles and Space Divi- 
sion, more than 5,000 scientists and engineers 
are engaged in one of U.S. industry’s broadest 
research and development programs. One group 
is conducting private industry’s largest, most 
diversified program of fundamental research in 
space physics. Already they have made massive 
contributions to America’s space technology— 
particularly in the Discoverer, MIDAS, and 
Samos satellite programs of the U.S. Air Force. 


MISSILES & SPACE DIVISION 


Sunnyvale, California 


i 
= 
— 


JOURNAL | 


A PUBLICATION OF THE AMERICAN ROCKET SOCIETY. 


EDITOR Martin Summerfield 


ASSOCIATE TECHNICAL EDITOR 
MANAGING EDITOR 

ART EDITOR 
CONTRIBUTOR 


Irvin Glassman 
Barbara Nowak 
John Culin 
George F. 


ASSOCIATE EDITORS 


Russian Supplement; Ali Bulent Cambel, North- 
western University, Book Reviews; Charles J. Mundo Jr., Raytheon Manu- 
_ facturing Co., Guidance; Bernard H. Paiewonsky, Aeronautical Research _ 
Associates of Princeton, Flight Mechanics; M. H. Smith, Princeton 
University, Technical Literature Digest = 


ASSISTANT EDITORS 
Julie Hight, Estelle Metz, Carol Rubenstein 


ADVERTISING AND PROMOTION MANAGER 
William Chenoweth 


ADVERTISING PRODUCTION MANAGER 
Walter Brunke 


ADVERTISING REPRESENTATIVES one 
Los Angeles 
James C. Galloway and Co. _—© 

6535 Wilshire Blvd., Los Angeles, Calif. 
Telephone: Olive 3-3223 

Detroit 
R. F. Pickrell and Vincent Purcell _ om 


Cc. and Associates 
Madison Ave., New York, N. Y. 
Telephone: Plaza 9-7460 


: 


_ Jim Summers and Associates 


E. Wacker Dr., Chicago, III. 318 Stephenson Bldg., Detroit, Mich. 
Telephone: Andover 3-1154 Telephone: Trinity 1-0790 
Robert G. Melendy John W. Foster 
17 Maugus Ave., Wellesley Hills, Mass. Sua + 239 4th Ave., Pittsburgh, Pa. 


Cedar 5 Telephone: Atlantic |-2977 


American Rocket Society 


500 Fifth Avenue, New York 36, N. Y. 


1930 
OFFICERS 
President Howard S. Seifert 
Vice-President Harold W. Ritchey _ 
Executive Secretary James J. Harford 
Treasurer _ Robert M. Lawrence 
Secretary and Asst. A. C. Slade 
Andrew G. Haley 


General Counse 


Director of Publications Irwin Hersey 


BOARD OF DIRECTORS 
Terms expiring on dates indicated Sis 
Ali B. Cambel 1962 William H. Pickering 1961 
Richard B. Canright 1962 Simon Ramo 1960 
James R. Dempsey 1961 William L. Rogers 1960 
Herbert Friedman 1962 David G. Simons 1961 _ 
Robert A. Gross 1962 7 Se John L. Sloop 1961 3 
Samuel K. Hoffman 1960 Martin Summerfield 1962 
A. K. Oppenheim 1961 Wernher von Braun 1960 | 
Maurice J. Zucrow 1960 _ 


TECHNICAL COMMITTEE CHAIRMEN 


David B. Langmuir, Ion and Plasma . 


Lawrence S. Brown, Instrumentation 
Propulsion 


and Contro 
Ali B. Cambell, Magnetohydrodynamics 
William H. Dorrance, Hypersonics 
James S. Farrior, Guidance and 
Navigation 
Herbert Friedman, Physics of the 
Atmosphere and Space 
George Gerard, Structures and Materials 
Martin Goldsmith, Liquid Rockets 


Max A. Lowy, Communications 

Irving Michelson, Education 

Peter L. Nichols Jr., Propellants and 
Combustion 

Eugene Perchonok, Ramijets 

Richard A. Schmidt, Test Facilities and 
Support Equipment = 

John I. Shafer, Solid Rockets ~ 


Andrew H. Haley, Space Law and 7 
Sociology C. J. Wang, Nuclear Propulsion 
Samuel Herrick, Astrodynamics Stanley C. White, Human Factors and ae 
Bioastronautics 


Maxwell W. Hunter, Missiles and Space 
Vehicles 
L. Karsch, Systems Operations 
and Sup 


George F. Wislicenus, Underwater 
ropulsion 


Abe M. Zarem, ) Power Systems 


Scope of ARS JOURNAL 

This Journal is devoted to the advance- 
ment of astronautics through the dissemina- 
tion of disclosing new scientific 
knowledge and basic applications of such 
knowledge. ‘The sciences of astronautics are 
understood here to embrace selected aspects 
of jet and rocket propulsion, spaceflight 
mechanics, high speed aerodynamics, flight 
guidance, space communications, atmospheric 
and outer space physics, materials and struc- 
tures, human engineering, overall system 
analysis, and possibly certain other scientific 
areas. The selection of papers to be printed 
will be governed by the pertinence of the topic 
to the field of astronautics, by the current or 
probable future significance of the research, 
and by the importance of distributing the in- 
formation to the bers of the S ty and 
to the profession at large. 


Information for Authors 


Manuscripts must be as brief as the proper 
presentation of the ideas will allow. Ex- 

lusion of ble material and concise- 
ness of expression will influence the Editors’ 
acceptance of a manuscript. In terms of 
standard-size double-spaced typed pages, a 
typical maximum length is 22 pages of text 
(including equations), 1 page of references, 
1 page of abstract and 12 illustrations. 
Fewer illustrations permit more text, and vice 
versa. Greater length will be acceptable 
only in exceptional cases. 

Short manuscripts, not more than one 
quarter of the maximum length stated for full 
articles, may qualify for publication as 
Technical Notes or Technical Comments. 
They may be devoted to new developments 
requiring prompt disclosure or to comments 
on previously published papers. Such manu- 
scripts are published within a few months of 
the date of receipt. 

Sponsored manuscripts are published 
occasionally as an ARS service to the indus- 
try. A manuscript that does not qualify for 
publication, according to the above-stated 
requirements as to subject, scope or length, 
but which nevertheless deserves widespread 
distribution among jet propulsion engineers. 
may be printed as an extra part of the Journal 
or as a special supplement, if the author or 
his sponsor will reimburse the Society for 
actual publication costs. Estimates are 
available on request. Acknowledgment of 
such financial sponsorship appears as a 
footnote on the first page of the article. 
Publication is prompt since such papers are 
not in the ordinary backlog. 

must be double spaced on one 
side of paper only with wide margins to allow 
for instructions to printer. Include a 100 to 
200 word abstract. State the authors’ posi- 
tions and affiliations in a footnote on the 
first page. Equations and symbols may be 
handwritten or typewritten; clarity for the 
printer is pee Greek letters and unusual 
symbols should be identified in the margin. If 
handwritten, distinguish between capital and 
lower case letters, and indicate subscripts and 
superscripts. References are to be grouped at 
the end of the manuscript and are to be given 
as follows. For journal articles: Authors 
first, then title, journal, volume, year, page 
numbers; for books: Authors first, then title, 
publisher, city, edition and page or chapter 
numbers. Line drawings must be clear and 
sharp to make clear engravings. Use black 
ink on white paper or tracing cloth. Lettering 
should be large enough to be legible after 
reduction. Photographs should be glossy 
prints, not matte or semi-matte. Each illus- 
tration must havea legend; legends should be 
listed in order on a separate sheet. 

Manuscripts must be accompanied by 
written assurance as to security clearance in 
the event the subject matter lies in a classified 
area or if the paper originates under govern- 
ment sponsorship. Full responsibility rests 
with the author. 

Preprints of papers presented at ARS 
meetings are automatically considered for 
publication. 

Submit manuscripts in duplicate (origi- 
nal plus first carbon, with two sets of 
rag to the Managing Editor, ARS 
JOURNAL, 500 Fifth Avenue, New York 


36, NY. 

ARS JOURNAL is published monthly by 
the American Rocket Society, Inc. and the 
American Interplanetary Society at — & 
Northampton Sts., Easton, Pa., U. A. 
Editorial offices: 500 Fifth Ave., New _ 
36, N. Y. Price: $12.50 per year, $2.00 per 
single copy. Second-class mail privileges 
authorized at Easton, Pa. This publication 
is authorized to be mailed at the special 
rates of postage prescribed by Section 132.122. 
Notice of change of address should be sent 
to the Secretary, ARS, at least 30 days prior 
to publication. Opinions expressed herein 
are the authors’ and do not necessarily reflect 
the views of the Editors or of the Society. 
© Copyright 1960 by the American Rocket 
Society, Inc 


ARS Journal 


— 


or 
4 
4 
Pe ~ 
> 
> = 
Ex 
> 
* 
> 
: 
4 


There are two ee to the ; STL coin... 


What STL does: 
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technical management of advanced ballistic 
missile and space systems. STL conducts 
advanced space flight experiments under the 
executive management of the Air Force 
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NASA. In addition STL’s leadership in a 
military applications of space technology is 
illustrated by its successful accomplish- 
ments as the contractor responsible for 
over-all systems engineering and technical 
direction of the Atlas, Titan, Thor, and 
Minuteman portions of the Air Force ef * 
Ballistic Missile Program. 


STL offers: 


For scientists and engineers with out- 
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growth opportunities in many areas of 
technical activity, including: 


Electronic and Electromechanical Systems 
Vehicle Engineering and Development 
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Systems Engineering and Technical 
Direction 


j Telecommunications 
+ Airborne Systems 
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positions at ort. now. Please send your resume today to Mr. Richard A. Holliday. 
»>~ Space Technology Laboratories, Inc., P.O. Box 95004 6p 
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The New 


Ramo-Wooldr ‘idge 


Laboratories 


in Canoga Park 


The new Ramo-Wooldridge Laboratories in Canoga 
Park, California, will provide an excellent environment 
for scientists and engineers engaged in technological 
research and development. Because of the high degree 
of scientific and engineering effort involved in Ramo- 
Wooldridge programs, technically trained people are 
assigned a more dominant role in the management of 
the organization than is customary. 

The ninety-acre landscaped site, with modern build- 
ings grouped around a central mall, contributes to the 


RAMO- WOOLDRIDGE LABORATORIES 


academic environment necessary for creative work. The 
new Laboratories will be the West Coast headquarters — 
of Thompson Ramo Wooldridge Inc. as well as house 
the Ramo-Wooldridge division of TRW. 

The Ramo-Wooldridge Laboratories are engaged in 
the broad fields of electronic systems technology, com- 
puters, and data processing. Outstanding opportunities 
exist for scientists and engineers. 

For specific information on current openings write 
to Mr. D. L. Pyke. 


8433 FALLBROOK AVENUE, CANOGA PARK, CALIFORNIA 
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Recent Advances in Chemical Kinetics 
Of Homogeneous Reactions in 
Dissociated Air 


related projects. 


S A MISSILE re-enters the atmosphere, it is surrounded 
by a detached shock wave. The air immediately behind 
the shock front becomes so hot that it partially dissociates 
into atomic oxygen and nitrogen. This dissociated air then 
diffuses toward the skin of the missile through a relatively 
cool boundary layer. The kinetics of possible homogeneous 
reactions in the boundary layer together with the transport 
properties of the gas mixture determine whether this disso- 
ciated gas recombines in the gas phase before impinging on 
the missile skin. The dissociated gas may well react much 
more rapidly with the missile skin than does the corresponding 
undissociated gas. For this reason, the kinetics of certain 
gas phase reactions may strongly influence the heat transfer 
rate to the missile, its aerodynamic characteristics and the 
ablation rate of the missile skin. Rosner (1)? discussed this in 
detail and reviewed several theories which considered the in- 
fluence of kinetic parameters on convection heat transfer. 
This review will supplement Rosner’s survey paper by examin- 
ing in some detail the kinetics of some of the reactions which 
are considered to be of importance. The influence of chemical 
kinetic parameters also should be considered in many other 
aerodynamic problems; for example, in the design of a rocket 
nozzle (2), or in the interpretation of experiments conducted 
in hypersonic low density wind tunnels. 

The predominant mechanism for the combination of atoms 
at ambient density and temperature is a three-body reaction. 
Two atoms simultaneously collide with a third particle and 
recombine; the third particle carries away some of the energy 
liberated during this recombination process and thereby pre- 
vents the reverse dissociation from occurring immediately. 
At sufficiently low densities, these three-body collisions be- 
come so uncommon, compared with two-body collisions, that 
radiative recombination becomes a predominant mechanism 
in the gas phase. In this case, two atoms collide and some- 
times emit part of the energy of recombination as a light 
quantum instead of again flying apart. In a gas containing 
only one element, these would be the only reactions of im- 
portance. In mixtures, bimolecular reactions (e.g., N + Oo—> 


Received Jan. 11, 1960. 

1 This research is based on work performed under the auspices 
of the U. 8. Air Force, Ballistic Missiles Division, Contract no. 
AF 04(647)-269 at General Electric Co., Missile and Space Ve- 
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NO + O), in which the total number of particles is conserved, 
couple the concentrations of the different dissociated gases 
to each other and are, therefore, of importance. Each of these 
processes is discussed in detail. 

Trotman-Dickenson (3) prepared an outstanding book re- 
viewing the research work on gaseous kinetics done prior to 
July 1954. However, he did not discuss radiative recombina- 
tion. Since that time the application of the shock tube to 
chemical kinetic studies has led to a new understanding of 
reactions at elevated temperatures. Several publications have 
also appeared recently on various reactions of dissociated nitro- 
gen and oxygen. For these reasons, a more current review ap- 
pears highly desirable. This review is not intended to be ex- 
haustive, but covers the work published prior to October 1959, 
which the author considers most pertinent in a boundary layer 
containing dissociated air. 

Since the terminology of chemical kinetics has been defined 
in several texts (4-9), it will not be covered in this paper. 


Three-Body Recombination 


For most simple atom recombination reactions, k, is of the 
order of 10'4-10!5 em® per gm-mol? see (2.6 X 10'° to 2.6 X 
101! ft® per lb-mol? sec) at room temperature. For the recom- 
bination of a given pair of atoms, k, depends on the choice of a 
third body. However, for most common third bodies, k, is in 
the range indicated at room temperature as will be illustrated. 
The exceptions to this rule are not likely to be encountered in 
problems of technical interest. 

As the temperature increases, k, decreases slowly. The dif- 
ference between k, at room temperature and k, at 1000 K is 
generally less than one order of magnitude. This is best il- 
lustrated by studies of iodine atom recombination and is 
corroborated by analogous work on bromine atom recombina- 
tion (10,11), which is not discussed further in this paper. For 
oxygen atoms, it will be shown that sufficient work has been 
done to indicate that the magnitude of the recombination co- 
efficient is about the same at 3500 K as at room temperature. 
The recombination of atomic nitrogen has been studied only 
at room temperature. However, it appears to be possible to 
predict the effect of raising the temperature by means of an 
analogy with iodine recombination. The effect of changes of 
k, with temperature is illustrated also by some observations 
on flames (12). This w ork ilh illustrates some applications of the 
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concepts developed herein and is, therefore, reviewed at the 
end of this paper. 

Several experimental determinations have been made of k, 
for simple atoms. At room temperature, k, can be measured 
directly by using flash photolysis, e.g., (13) or discharge tubes, 
e.g., (14). At very high temperatures, k, is determined in- 
directly by means of a shock tube, e.g., (15). 

With the shock tube, kz is determined directly from the 
initial temperature and pressure, shock velocity, and measured 
density or concentration ‘profile; k, is then determined by 
dividing ka by an appropriate equilibrium constant. Because 
ka and the equilibrium constant are both much more sensitive 
to temperature than k,, it has not been possible to determine 
experimentally the functional dependence of any given k, on 
temperature. 

Two theories have been proposed which predict that the re- 
combination coefficient will decrease with temperature, and 
which are in qualitative agreement with experiment. One 
theory (10) is that when two atoms recombine, some of the 
energy of recombination will be dissipated through the inter- 
nal degrees of freedom of the newly formed molecule and the 
third body. This theory predicts that the recombination co- 
efficient will be inversely proportional] to a half-integral power 
of the absolute temperature. A second theory (15-17) is that 
an atom forms an intermediate compound with the third body. 
When this compound collides with a second atom, it reacts to 
yield a molecule and the original third body. This theory 
predicts that the recombination coefficient will have a low 
negative energy of activation and will be exponentially de- 
pendent on temperature. A more general theory (18) is that 
the rate limiting step in recombination changes with tempera- 
tures. At low temperature the reaction rate is determined 
by the collision frequency and, therefore, increases with the 
square root of the absolute temperature; on the other hand, 
at very high temperatures, the rate limiting step will be one 
of those previously listed so that the recombination coefficient 
will tend to decrease with temperature. This theory is 
plausible; however, it has not been verified experimentally. 


Iodine Atom Recombination 


The recombination of iodine atoms has been studied ex- 
tensively because of its experimental convenience. The low 
dissociation energy of molecular iodine makes it feasible to 
dissociate iodine by either photolysing it with visible light or 
by subjecting it to relatively weak shock waves, and the spec- 
trum of molecular iodine permits one to monitor its concen- 
tration by absorption of the visible light. The results of this 
work give insight into reactions of greater technical interest. 

At 298 K for iodine atoms with argon as a third body, k, is 
(2.9 + 0.3) X 10" cm per gm-mol? sec (13). k, is shown to 
have the same order of magnitude for most atoms with argon 
as the third body. At this relatively low temperature, iodine 
is at least 250 times more effective than argon as a third body. 
(L.e., the rate constant with iodine as a third body is more than 
250 times greater than with argon as a third body.) 

This reaction has been studied also at elevated tempera- 
tures in a shock tube (15). The results are tabulated in Table 
1. At 1300 K, it was found that iodine was less than 35 times 
as effective a third body as nitrogen. 

Two very important general results which have been found 
to be applicable to all atom recombination reactions are il- 
lustrated by this work: 

1 The recombination coefficient with argon as a third 
body decreases very slowly with temperature. The difference 
between the recombination coefficient at 1300 K and that at 
297 K is less than one order of magnitude. 

2 From Table 1 and from the effect of temperature on the 
efficiency of iodine as a third body, it can be seen that the 
differences between the efficiencies of various third bodies 
tend to vanish at elevated temperatures, as suggested in (15). 
Because of this, molecular oxygen, molecular nitrogen and 
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argon probably can be considered equivalent third bodies at 
elevated temperature. To within one order of magnitude, 
atomic oxygen and atomic nitrogen are probably equivalent 
to the respective undissociated molecules as third bodies. A\- 
though the theoretical explanation for these results has not 
been established, all of the current plausible theories dis- 
cussed previously predict that the recombination coefficient 
becomes less dependent on the absolute temperature as the 
temperature is elevated (10,15-17). Thus, the recombination 
coefficient for atoms probably can be considered independent 
of temperatures above 1000 K for the third bodies encountere: 
in a normal atmosphere. 

In Table 1 as regards N» it has been assumed as it was als» 
by Britton et al. (15) that the reaction of interest occurre:| 
so rapidly that the nitrogen remained vibrationally unr - 
laxed during this time. At 1300 K, the relaxation time for 
iodine dissociation under the conditions of Britton’s exper: - 
ments is 75 microsec. By contrast, Blackman’s results (1‘ 
indicate that at the temperature and pressure at which th:< 
experiment was conducted, pure nitrogen has a vibration«! 
relaxation time of 1700 microsec. The possibility must b» 
considered that vibrational relaxation is much more probab)» 
when a nitrogen molecule collides with an iodine atom or 
molecule than when it collides with another nitrogen molecul: . 
Since the atomic weight of iodine is greater than the molecula: 
weight of nitrogen, the parameter (x) will be even smaller for 
the above collision than for a collision between two nitroge:: 
atoms. Therefore, according to the theory of Landau an:| 
Teller (20), it is even less likely that the vibrational energy of 
an excited nitrogen molecule will be dissipated as translational! 
energy when it collides with an iodine particle than when it 
collides with another nitrogen molecule. 

A vibrational quantum in molecular iodine (21) has 214.6 
of energy whereas that in nitrogen has 3394 cm~!.  Be- 
cause of this large difference, resonance transfer of vibrational! 
energy between nitrogen and iodine is improbable. There- 
fore, the assumption of Britton et al. seems plausible. 4 q 


Oxygen Atom Recombination 


At room temperature, the recombination of oxygen atoms 
might be studied under two extreme conditions. Either the 
ratio of oxygen atoms to oxygen molecules is much greater 
than 1, or much less than 1. In the former case, the following 
mechanism of recombination will predominate ‘yeh tole 


Golden and Myerson (22) have studied the recombination of 
atomic oxygen under this condition by using an electric dis- 
charge to dissociate oxygen in a static system and then follow- 
ing the recombination with flash spectroscopy starting milli- 
seconds after their discharge is extinguished. Their pre- 
liminary results indicate that k, for the reaction described is 
3 X 10'4 em® per gm-mol? see at room temperature with oxy- 
gen as a third body. This constant is not unusually large for 


Table 1 Recombination coefficients for iodine atoms 


© with different third bodies 
kr, per k; argon 


kr gas 


Third em®/gm- at temp. per k; argon 
body mol sec at room temp. 
A 2.9 X 10% 298 ied Bo 

A 4.5 X 10 1300 1 1 

He 1.8 X 1400 0.36 

4.4 10% 1300 1.0 1.25 

O. 5.3 X 10% 1275 1.2 1.8 

I, 1100 less than 35 250-600 
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Table 2 Rate constants for oxygen atom recombination 
by formation of ozone 


Rate constant 

(4.61 + 0.25) X 

1015¢— 24,000/R Tee /gm-mol 

sec 
(6.00 + 0.33) XK 1013 X 

e+ 600/RTom*/gm-mol? sec 
(2.96 + 0.21) X 1013 X 

e— 6000/RT ee /gm-mol sec 


Reaction 
1 O; + 0;— O; +0.+0 


2 O;+0.+0—-0; + 0; 
3 O + O; — 202 


r-combination. It indicates that the reaction O + O. —~ Os; 
followed by O; + O — 20, is unimportant under their ex- 

perimental conditions. 
By contrast, if the ratio of molecular oxygen to atomic 
oxygen is much greater than 1, collisions of the type O + O2+ 
/f — Os + M will be much more common than collisions 
the typeO +O + M. Under this condition, 


ti ie mechanism for recombination will be iv Fo - 
ka 
0+0,+ M 
A 
@ 


k 


These reactions can be studied through the thermal decom- 
position of ozone. Benson and Axworthy (23) made such a 
study; their results are given in Tables 2 and 3. 

Schiff (24) and his co-workers have recently studied the re- 
combination of oxygen atoms under these conditions. They 
began their measurements at least } sec after the oxygen had 
been dissociated by a discharge; under their experimental 
conditions, the O, to O ratio was at least 20 to 1. Their re- 
sults were in accord with the described mechanism and with 
the rate constants of Benson and Axworthy. 

At elevated temperatures, this mechanism will become less 
important both because O; becomes more unstable and be- 
cause O will probably be present in larger concentration. 
Matthews (25) has recently found that k, for oxygen atoms, 
with molecular oxygen as a third body, is (8.4 + 3) X 10" 
cm® per gm-mol? sec at 3500 K. This is an outstanding piece 
of work which illustrates how effectively a shock tube can be 
used in kinetic work. Because oxygen in contrast to iodine 
does not have a convenient absorption spectra in the visible, 
it was necessary to obtain the profile of the oxygen atom con- 
centration from the density profile. This profile was de- 
termined from the profile of the index of refraction which 
was measured by interferometry. Matthews was able to 
show that the experimental density profile agreed with the 
calculated profile at seven points. 

This work is corroborated by an earlier study of Chesick 
and Kistiakowsky (26) who measured the density profile of 
two shocked 75 per cent oxygen-25 per cent xenon mixtures 
by x-ray absorption. Within the experimental error their rate 
constants agree with those of Matthews, provided that the 
reasonable but very tacit assumption of Chesick and Kistia- 
kowsky is accepted, i.e., that xenon cannot be much more 
effective than oxygen as a third body. This agreement is 
especially impressive because the technique of Chesick and 
Kistiakowsky does not require knowledge of the index of re- 
fraction of a partially dissociated gas at elevated temperature. 

Byron (27) has recently published a shock tube study of 
oxygen dissociation similar to, but less precise than, that of 
Matthews. Byron’s resolving time was 1 microsec, whereas 
Matthews’ resolving time was 0.25 microsec. Instead of 
using the complete density profile, Byron determines the rate 
constants from the initial slope of the density profile im- 
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mediately behind the shock front, and from the time (é:,) it 
took for the density to reach the average between the initial 
value immediately behind the shock and its ultimate limiting 
value at a point far enough behind the shock for thermody- 
namic equilibrium to be achieved. From the initial slope of 
the density profile, Byron determines k, with molecular oxygen 
as a third body. This k, lies at the lower limit found by Mat- 
thews. Then from t:;,, the effect of atomic oxygen is found. 
Byron reports that atomic oxygen is 10 times more effective 
than molecular oxygen as a third body between 2000 and 5000 
K. However, at equilibrium, Byron calculates that under his 
most favorable condition, 11 per cent of oxygen was disso- 
ciated. At t:y,, there is only about 5 per cent dissociation at 
most. Therefore, collisions involving molecular oxygen occur 
10 times as frequently as those involving atomic oxygen. The 
author feels that Byron’s error in the rate constants is larger 
than Matthews’, and that it is completely reasonable to as- 
sume that Byron’s value of k, with molecular oxygen as a 
third body should be doubled, bringing it into agreement with 
Matthews’ result. Then the same t:y, would be found by as- 
suming that molecular and atomic oxygen were about equal 
as third bodies. 


Nitrogen Atom Recombination 


A recent study (14) shows that k, for nitrogen atoms with 
molecular nitrogen as a third body is (5.7 + 0.8) X 10" cm® 
per gm-mol? sec between 195 and 450 K; with argon as a 
third body, this coefficient is (2.8 + 0.3) X 10" ec? per mol? see 
at 297 K. This is in fair agreement with an earlier piece of 
work (28) which is less precise as will be shown. 

In both studies, nitrogen atoms were generated in a flow 
system by an electric discharge, and the concentration of 
atoms was then determined as a function of distance. Know- 
ing the flow velocity, the rate constant could then be deter- 
mined. 

In the more recent study, the N atoms could be titrated 
with nitric oxide (28,29) at five points in the apparatus. This 
titration was followed by a mass spectrometer. By contrast, 
in the earlier work, the nitrogen atoms were titrated at only 
one point and the end point was determined visibly. The 
author has found in his own work that the end point deter- 
mined in this way is very subjective. The intensity of the re- 
sultant afterglow was then used to measure the atom concen- 
tration. In the more recent work, total pressures between 1 
and 10 mm Hg were used. It was shown directly that at 
higher pressures the nitrogen atoms disappeared by a third- 
order reaction, and therefore heterogeneous recombination at 
the wall of the apparatus was negligible; at lower pressures, a 
first-order reaction due to recombination at the wall had to be 
considered. By contrast, the earlier workers used pressures 
between 0.1 and 1 mm Hg and just assumed that wall recom- 
bination could be neglected. 

No direct study has been made on the recombination of 
nitrogen atoms at elevated temperatures. However, note that 
the value of the recombination coefficient of nitrogen atoms 
at 297 K with argon as a third body is the same as the cor- 


receanbination coefficient for iodine atoms. It is 


ozone in reactions 1 and 2 in Table 2 


Relative efficiency 


1 
0.44 
0.41 
1.06 
0.34 
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believed, therefore, that the values of the recombination co- 
efficient for nitrogen atom recombination at elevated tem- Table 4 Calculated and experimental activation energies 
perature are probably equal to the corresponding recombina- for several reactions 
tion coefficients for iodine atoms. However, the author feels — _ ait 
that it is desirable to confirm this experimentally. = | — Activation energy, 

Radiative Recombination Reference values values Source 

NO+N 8.3 very small (29) 

It has been pointed out, using a theoretical argument, that O. +N 6.5 6.2 (29) 
radiative recombination will occur in less than one out of one i 4-H, 6.2 7s 4 40 (39) 
million collisions between atoms in the most favorable case Cl + H: 6.2 5.5 (40) 
(30). Since the collision rate is about 1014-7 ce per gm-mol sec NOCI + Cl 2.0 1.2 (41) 

(31), the rate constant for radiative recombination has an O + NO; 3.9 very small (42) 
upper limit of about 10° ce per gm-mol sec. It can be seen 0 +0; 1.25 6 (23) 
by appropriate substitutions that radiative recombination 

might become important when the particle density is less than 

one thousandth of its value at sea level. 

However, various restraints which can be derived, at least 
in theory, by application of quantum mechanics generally will the atoms can be lumped together and recombination can | 
lower the radiative recombination coefficient by many orders considered to occur by the following mechanism 
of magnitude. For example, the radiative recombination of 
oxygen atoms and carbon monoxide to form carbon dioxide A+A+M—~A:+M 
must occur through a transition not favored by selection rules Relations [6 and 7] replace two kinetic equations. Of course. 
derived through quantum mechanics (32). The rate constant by multiplying Equations [6 and 7], it is found that the rela- 
for this reaction is 1.6 X 10° cc per gm-mol sec (33). tion 

At present, it has not been shown that this mode of recom- 
bination is of sufficient technical importance to an aerody- K, = = (Po)*(Px.)( (Px) ~? 
namicist to warrant a more detailed review. 

is also satisfied. This was previously indicated. 


= . . . 
- To determine which bimolecular reactions are rapid enoug): 
Bimolecular Reactions mp Phen a ao to achieve equilibrium, the appropriate rate constants must 


certain equilibria may be achieved even in a nonequilibrium ~ 
(94,38). For ‘example alt exponential factor is generally between 10!* and 10" ce per 
the temperature is changing very rapidly, the equilibrium gm-mol mites However, because of ynanee mechanical con- 
relation (Px,)(Po)*(Px)-*(Po) = K (T) ney wager eee siderations reactions between materials of low molecular 

Ne 2 1 


even though the equilibrium relations (Po)? (Po,)-! = K3(T) weight, in which electron spin (angular momentum) is not 
and (Py)?(Px.)~! = K;(T) are not satisfied. This illustration conserved, have much lower pre-exponential factors (37). 
If the reaction is exothermic, the activation energy is 


Consider a mixture of nitrogen atoms and oxygen atoms in generally less than 10,000 cal per gm mol and can be approxi- 


equilibrium at high temperature (7;). Assume that the tem- mated by the Hirschfelder rule (38) which states that the 
perature is rapidly reduced to a lower temperature (72) so activation energy is about 5.5 per cont of the energy of the 
that the degree of dissociation at equilibrium is now negligible. bond between 6 and he (To obtain the rate constant for 
The problem is to find out how fast the atoms will recombine the endothermic reaction in the opposite direction, divide the 
A priori, one way of doing this is to assume iakdeneiitnidion previously mentioned rate constant by the equilibrium con- 


mee the following reactions Several comparisons between experimental activation ener- 
Pit Ss N+04+M-—-+NO+4+M [1] gies and activation energies calculated from Hirschfelder rules 
are given in Table 4. Most of the dissociation energies neces- 
Toes [2] sary to calculate the activation energies and heats of reaction 
N+N+M—-NL+M [3] in Table 4 were taken from Gaydon (43). For the reactions 
NOC! + Cl and O + NO., the data were taken from the 
It is then necessary to solve three kinetic equations simul- Bureau of Standards Table (44). Several other comparisons 
taneously. The previously illustrated treatment ignores the between experimental results and results computed from 

following two shuffling reactions aul Hirschfelder rules are given in (45). 
ae wi > A close examination of the table shows that some very exo- 
[4] thermic reactions actually have lower activation energies than 
[5] predicted by Hirschfelder rules. This is expected from 


theoretical considerations discussed in the literature (46). 


Under many conditions which can be determined from rate 4 : é : : 
Some specific reactions are discussed in greater detail: 


constants next described, reactions [4 and 5] will be so rapid 


compared with reactions [1, 2 and 3] and with the rate at N +NO—>N.+0 
which the temperature is changing, that the following equi- 
librium relations will be obeyed in the recombining mixture Because this reaction is extremely rapid, only a lower limit 
of 5 X 10" ce per gm-mol sec (47) can be set on its rate con- 
(PxoXPo) _ x (Ts) (6 stant at room temperature. 
N 
(PuoXPx) = [7] This reaction has also been studied experimentally (29). 
It was found that k equals 2 X 10'* exp (—6200/RT) 
If the first three reactions occur at about the same rate (as cc per gm-mol sec between 394 and 515 K. 
indicated, this is probably a reasonable assumption), then all Reactions of the type Ae + B, —~ 2AB are generally very 
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These rates slowly decreased with temperature. 


slow compared with other shuffling reactions. This is il- 
lustrated by the reaction below 


N2 + O2. + 2NO 


This reaction is so slow that its rate constant has never been 
determined. An effort (48) to do this by measuring the rate 
of formation of NO in air between 2000 and 3000 K showed 
that oxygen atoms formed by the thermal dissociation of 
molecular oxygen reacted with molecular nitrogen so rapidly 
to yield NO and atomic nitrogen (i.e.,O + Ne—~ NO +N), 
that the effect of the described reaction could not be observed. 
Therefore, the reaction rate can be considered virtually zero in 
all aerodynamic problems. This conclusion has been verified 
by direct studies (49,50) of the reverse reaction whose activa- 
tion energy was found to be between 63 and 64 kcal per gm 
mole. By addition of the heat of the previously described 
reaction to the activation energy of the reverse reaction, the 
energy of activation of this reaction is found to be between 106 
and 107 keal per mole. This is extremely high. 


Application 


Some of the ideas expressed in this paper were originally 
suggested by Sugden and his colleagues who have utilized 
them in their work on flames. He has burned fuel rich mix- 
tures of hydrogen, oxygen and nitrogen (12), and found that 
the gas rapidly reacts exothermally to yield a product con- 
taining an excess of H and OH radicals which slowly recom- 
hine as the gas streams from the burner. In this mixture, the 
following reactions are responsible for recombination 


H+H+M—H:.+M 
H+0H+M—>HO+M 


The reaction 
OH +0H 


is unimportant because H.O. is unstable at the temperatures 
used (1600-2400 K). 

Because of two-body reactions, the concentrations of H 
and OH are coupled and the following equilibrium is main- 
tained 


20H +H, 

This principle was utilized to measure the concentration of H 
in the flame. Flame gas containing less than 1 ppm of NaCl 
was used. This was not sufficient to perturb the flame. How- 
ever, the sodium atom concentration would be measured 
from the intensity of its resonance line. This concentration 
was coupled to the hydrogen atom concentration because the 
equilibrium Na + HCl = NaCl + H was established in the 
flame. 

At 1650 K, the following rate constants were found 


ky = 5.4 X 107 em®/gm-mol? sec 


ks = 2.2 X em*/gm-mol? sec 


Sugden (12) felt that this work showed that H,O was 10 
times as efficient as H, or Nz in promoting the recombination 
of H and OH and that he could, therefore, neglect the effect of 
three-body collisions with H, or N2 as third bodies in the burnt 
gas. He supported this argument by showing that the prod- 
ucts (k2)(H,O) and (k3)(H:O) were independent of the com- 
position of the burnt gas. The burnt gases used contained 
from 18 to 36 mole per cent H.O. The remainder was various 
mixtures of nitrogen and hydrogen. Assuming that hydrogen 
and nitrogen had the same efficiency as water in promoting 


Marcu 1960 


these recombinations, (k2) (H,O) would change by a factor of 2 
over the composition range used. Unfortunately, the authors 
do not show that they can detect this difference, which is 
commonly the experimental error in kinetic work. Therefore, 
all that can be said is that H,0 is at least as efficient as H, and 
N2 in promoting the illustrated recombination reactions. The 
rate constants given may, therefore, be as much as five times 
too large because the effect of N, and He was neglected. 


Summary 


The kinetics of the homogeneous reactions involving mix- 
tures of diatomic gases and their dissociation products are 
important in many aerodynamic problems. There are two 
general types of reactions of importance, atom recombination 
(A + A + M— A, + J) and its reverse reaction (dissocia- 
tion), and two-body reactions (A + BC—~ AB+(C). 

For many atom recombination reactions of technical im- 
portance, k, is about 10‘ em* per gm-mol? sec at room tem- 
perature and decreases slowly as the temperature is increased. 
Its value at elevated temperature generally differs by less than 
one order of magnitude from its value at room temperature. 

For many two-body reactions as illustrated, k = 10'%e~ 2/#7 
ce per gm-mol sec where F is less than 10 keal per gm mol for 
an exothermic reaction. These two-body reactions couple the 
concentrations of different dissociated gases to each other so 
that it is frequently possible to assume that thermodynamic 
equilibria of the type A» + 2B 22A + Bare achieved even 
though the equilibria 2A = A, and 2B = B, are not reached. 
An application of this to a study of flames is given. 
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Nomenclature 
® 
E = activation energy 
K = equilibrium constant a 
P = pressure 
R = universal gas constant (cal/gm-mol K) 
T = absolute temperature 
a = pre-exponential factor 7 
e = base of natural logarithms 
k = rate constant ar 
kr = termolecular recombination rate constant 
ka = dissociation rate constant 
x = ratio of the time for a molecule to make one vibration to the 
time two molecules remain in contact during a collision 
A = sum ofall atomic species in a gas 


A» = sum of all diatomic species in a gas 7 
M = third body in a recombination reaction © 
N = atomic nitrogen 
Ne. = molecular nitrogen 
NO= nitric oxide 


O = atomic oxygen 
O2 = molecular oxygen 
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N ORDER to stabilize a flame in a high velocity stream of 
combustibles, it is necessary to provide a continuous igni- 
tion source. The recirculation zone associated with an ob- 
stacle placed in the flow is a low velocity region of hot, burned 
gases which ignite the reactants continuously. The obstacle 
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Conditions in the critical zone have a 
stationary by a small jet directed counter to the main stream. Using the spectral intensity ratio 
method, the local equivalence ratio in the critical zone was measured. The results suggest that 
there is one stability curve (blowoff velocity as a function of the equivalence ratio in the critical 
The equivalence ratio in a reverse jet stabilized flame was de- 
termined locally as a function of duct radius. Jet composition affects the observed equivalence 
ratio radially as far as the duct walls although only weakly beyond about one third of the duct 


zone) regardless of jet composition. 
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strong influence on the stability of a flame which is held 


may be a bluff body, an array of bluff bodies, a recess in the 
duct wall or a small jet opposing the main flow. A flame 
stabilized by a reverse jet has been studied using the spectral 
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intensity ratio method which presumably gives the local 
equivalence ratio. The restrictions on the application of this 
method will be discussed. 

Reverse jet flame stabilization introduces several new 
variables and hence new design possibilities. Among the new 
variables are jet momentum, composition, temperature and 
mass flow rate. Since the jet has a pronounced effect on flame 
stability, the jet variables may be optimized to fit the widely 
varying operating conditions in ramjets and afterburners. 


Experimental Evidence Indicating the Existence 
of a Critical Zone (1-6)* 


Calculations show that the reverse jet contributes only 
alout 0.01 of the mass flow of the main stream. Yet changing 
the jet composition from pure air to an air-fuel mixture pro- 
duces a large effect in the stability curves (blowoff velocity as 
a function of equivalence ratio). 

In order to account for the pronounced change in the sta- 
bility curves, it seems reasonable to assume that there is a 
critical zone which determines flame holding properties. The 
jt flow, although only a very small fraction of the main stream 
fi»w, must contribute a large fraction to the flow entering the 
“«ritical zone.’’ Numerous investigators (3,6) have proposed 
a model in which the critical zone is a region of turbulent 
n.ixing located directly behind the nose of the flame. Agoston, 
Noon and Witherly (11) agree with the critical zone concept 
but feel that it is a torroidal recirculation region with the jet 
as the axis. 

Optimum flame stability occurs for a stoichiometric jet and 
niain stream. When the main stream is not stoichiometric, 
Noon (5) found that the equivalence ratio in the critical zone 
¢, at the maximum blowoff velocity V*,, is unity, regardless 
of jet temperature. Penner and Williams (6) have extended 
the concept by stating that at V*,., ¢, is a constant character- 
istic of the fuel-air system. 

Experimental results will be shown which suggest that 
actually there is one stability curve relating V;,, to ¢. for a 
given fuel and duct geometry. 


Spectral Intensity Ratio Method 


The ratio of the apparent peak intensities of two rotational 
lines belonging to different emitters, such as CH and C; can 
be correlated with the local equivalence ratio. To apply the 
spectral intensity ratio (SIR) method for determining local 
equivalence ratios, it is necessary to measure SIR as a function 
of » using a known fuel-air input to a Jaminar flame. In 
order to apply the calibration curve with some degree of con- 
fidence, the manner in which changes in operating conditions 
affect SIR should be examined. 

If the emission from optically thin gases is observed with a 
spectroscope having a slit width large compared to the line 
width, then it can be shown that SIR is proportional to the 
ratio of the population of emitters in the excited states. 

In a flame the excitation is not necessarily thermal but may 
be the result of chemical reactions which produce molecules in 
excited states (chemiluminescence). In particular, the radia- 
tion from Cz. and CH must be associated with the combustion 
reactions. The SIR is assumed to be a function only of the 
local equivalence ratio. This same assumption has been 
made without proof also by several other investigators. 

John and Summerfield (7) have measured the effect of tur- 
bulence on the radiation intensity from propane-air flames. 
Turbulence may reduce the specific intensity (radiation in- 
tensity divided by fuel flow rate) as much as 35 per cent for 
fuel-lean mixtures. Since the relative changes of C. and CH 
intensity for a turbulent flame enclosed in a duct are nearly 
equal, turbulence does not alter the correlation of SIR with ¢. 

A test for the application of the calibration curve obtained 
from a laminar flame to the reverse jet flame may be made by 


8’ Numbers in parentheses indicate References at end of paper. 
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Fig. 1 The observed equivalence ratio as a function of radial 

distance from the duct axis for a stoichiometric main stream and 

stoichiometric jet. Duct: 1}-in. standard pipe. Jet velocity: 
480 fps. Main stream velocity: 42 fps 


using a stoichiometric main stream. The results of these tests 
are shown in Fig. 1. Each point shown is the result of several 
different SIR measurements. The average is shown as a 
circle and the range of values are indicated by the lines. An 
inference from these data is that SIR does not necessarily 
yield correct absolute values for the reverse jet flame. How- 
ever, relative values of ¢ are probably given correctly, be- 
cause changes in jet composition produce the expected changes 
in the measured values of ¢. Fortunately, the interpretation 
of the data does not depend in an important way on the ab- 
solute value of ¢. 

Radiation measurements were made by Clark and Bittker 
(8) to determine the extent to which radiant intensity could be 
used to measure flame velocities. There is a relation between 
surface area and radiation intensity for laminar flames burn- 
ing at fixed equivalence ratio; the radiation per unit area of 
a laminar flame depends only on the equivalence ratio. Using 
the preceding conclusion, John, Wilson and Summerfield (9) 
correlated spectral intensity ratio with the equivalence ratio 
of a laminar flame stabilized by a bluff body in a duct. To ob- 
tain the correlation, interference filters with a band width of 
about 100 A were used. One filter was centered to pass the 
emission from CH produced by the transition A — *II(0, 0) 
and the other filter transmitted the emission from C2 associated 
with the transition *I] — ‘II (1, 1). Using the correlation of 
SIR with ¢, John, Wilson and Summerfield (9) investigated 
preferential diffusion in flame stabilization. Later Wilson, 
John and Summerfield (10) studied the effect of boundary 
layer flow on flameholding by using the SIR method. 


Experimental Determination of Local 
Equivalence Ratio in a Reverse Jet Stabilized 
Flame 


In this section the apparatus and procedure for obtaining 
the calibration curve and the local equivalence ratio are de- 
scribed. 


Experimental Determination of Spectral Intensity Ratios 
as a Function of Equivalence Ratio 


A laminar flame from a welding torch served as the calibra- 
tion flame. An image of the inner cone was formed on the slit 
of a diffraction grating spectroscope. Masking the slit pre- 
vented radiation from the tip and the base of the inner cone 
from entering the spectroscope. Spectral intensity ratios de- 
termined from any part of the straight portion of the inner 
cone were the same for identical fuel-air input to the flame. 
(See Table 1.) The lines given in Table 1 were selected be- 
cause they are well isolated, and because, for a given equiva- 
lence ratio, the spectral intensity ratios based on these lines 
were reproducible to an accuracy of about 5 per cent. 
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Fig. 2 Apparent emission intensity as a function of wave length 
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Fig. 3 Spectral intensity ratio as a function of equivalence 
ratio for the line with center at 4 A to the line with center at 
4324 


240 


Table 1 Rotational lines selected for the computation of 
intensity ratios 


Transitions Emitter 


2A — *11(0, 0) 
2A 11(0, 0) 


Wave length, A 
; 
4252 


4324 = 0) 
311 3{1(1, 0) 


TIC, 1) 


1) 
6191 511(0, 2) 


In the experimental studies steady-state conditions were 
established at a fixed equivalence ratio. A wave length in- 
terval including the line was then scanned. A typical record- 
ing of apparent intensity as a function of wave length is show: 
in Fig. 2. These data refer to the electronic transition *II, t) 
3]I, and the vibrational transition 1 — 0 of the radical C 
The height of the peak was measured above the continuun 
(See Fig. 2.) 

Typical results of the calibration from the laminar flame 
are shown in Figs. 3 and 4. It is important to recogniz: 
that the calibration curves can be applied only to measurc- 
ments of SIR made by the same instrument. The spectro- 
scope in this case was a Jarrell Ash diffraction grating spectro- 
scope using an RCA 1P28 photomultiplier tube for detection. 
The spectroscope has an f number of 6.5 and a linear disper- 
sion of 18 A/mm. Wyre 


equivalence ratio, 


Fig. 4 Spectral intensity ratio as a function of equivalence 
ratio for the line with center at 4737 A to the line with center at 
4252 
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Description of the Apparatus and Procedure Used for Ob-— 
taining Local Equivalence Ratios 


In order to view the flame, particularly the critical zone, 
and to measure emission intensity from limited portions of the 
flame, it was necessary to use a mirror system as shown in Fig. 
5. A mirror was located on an adjustable mount within the 
plenum chamber so that the flame was viewed along the axis 
of the duct. 

Commercial grade propane (95 per cent propane with higher 
hydrocarbons as impurities) was used as fuel. Facilities were 
provided for a jet of air, oxygen, nitrogen or fuel-air mixtures. 
Ignition was accomplished by means of an automobile spark 
plug. Two duct sizes were used: 1}-in. diameter and 2-in. 
diameter standard pipe. Mixing of the fuel and air in the 
main stream was assured with a mixing chamber which had 
a porous metal (Poroloy) sheet across the flow. The diameter 
of the steel tubing used for the jet was 7g in. 

Steady-state operating conditions were established in the 
duct. By adjusting the mirror and masking off a portion of 
the spectroscope slit height, the radiation from a small portion 
of flame surface was selected for study. The size of the seg- 
ment of the flame surface was an area about % in. long by 
several hundred microns wide. The radiation entering the 
slit was scanned with the spectroscope; a large enough wave 
length interval on either side of the line was taken so as to 
determine the magnitude of the continuum radiation. 


Experimental Results 


Two types of experiments were conducted. The equiva- 
lence ratio in the critical zone at blowoff was determined, 
and the local equivalence ratio as a function of radial position 
was measured. 


Equivalence Ratio in the Critical Zone at Blowoff 


With an air jet, the observed equivalence ratio in the critical 
zone ¢, is less than the main stream equivalence ratio gms. 
For the rich blowoff limit, g, is considerably smaller than the 
main stream equivalence ratio. (See Fig. 6.) 
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Fig. 5 Schematic illustration of the apparatus used for studies of a reverse jet stabilized flame — 
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Fig.6 Observed equivalence ratio in the critical zone at blowoff 

for an air jet. Open points are the observed equivalence ratio 

in the critical zone; the solid points are main stream values. 

Duct size: 1j-in. standard pipe. Jet velocity : 480 fps 


Addition of fuel to the jet shifts the blowoff curves toward 
the lean side. Substitution of pure oxygen for air in the jet 
shifts the blowoff curves toward richer mixtures. This agrees 
with the results which other investigators have reported (1-5). 
When the jet is a stoichiometric fuel-air mixture, then ¢, is 
larger than ¢», for lean fuel-air mixtures in the approach 
streams. This result is illustrated in Fig. 7 by the points 
represented as squares. As with an air jet, the pure oxygen - 
jet gives ¢, less than g»;. No data were obtained for an oxy- 
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gen jet or an air-fuel jet at the rich blowoff limit. When the equivalence ratio of 0.58, for any air jet gms is 0.78, and for an 
main stream is rich and the jet is pure oxygen, the flame tem- oxygen jet gm, is 0.82. But for all three jet compositions, the 
peratures are very high, making it impossible to achieve a local equivalence ratio in the critical zone ¢, is approximately 


steady-state temperature without damage to the equip- 0.67. The flame geometry in the region of critical zone de- 
ment. pends primarily on the ratio of jet velocity to main stream 


The conclusion was previously stated that in the critical velocity. This is a fixed ratio for the data presented in Fig. 7. 
zone the local equivalence ratio at V*, was the same for any Under these flow conditions the minimum equivalence ratio 


jet composition. An extension of this idea is the statement at which the critical zone can maintain a flame is g, = 0.67. 
that at blowoff, for any jet composition, there is one stability 
curve provided V, is plotted as a function of g,. In Fig. 7 sf Radius 


the values of Vi. have been plotted as a function of gm, and 


yg; All points for ¢, fall nearly on the broken line curve. The equivalence ratio was measured as a function of radius. 


This is consistent with the idea of one stability curve Vio = These data yield a comparison of the equivalence ratio in th 
Vio(y.) regardless of jet composition. This statement merits critical zone with the equivalence ratio in other regions of th 
additional discussion. As an example, consider Vs, = 50 fps. flame. When the jet stabilized flame is viewed parallel to th 
In Fig. 7 for an air-fuel jet blowoff occurs at a main stream duct axis, the depth of luminous gas increases with increasin:: 
radius. (See Fig. 5.) Experiments were performed t: 


establish the equivalence ratio gradient for a stabilized flam 
and also the gop.-gradient (spectroscopically determine: 


values) for a pilot flame. 


In Fig. 8 curves of the observed equivalence ratio as a func 
tion of radial distance from the duct axis are presented. Fo: 


oxygen jet LLL G 
a rich main stream ¢,», remains constant to a radius R/3. 
Here is the duct radius. At larger radii it decreases. 
a / 4 air jet and the oxygen jet give results which are practically in- 
sieinmatete / y distinguishable. When fuel is added to the jet, gobs is higher 


ir-fuel jet 
aa : by about 5 per cent near the center of the duct; near the duct 
portion of stability curves ° 
plotted here walls, gobs merges with the other curves. The fact that gobs for 
iiiniitienk Wade the air-fuel jet is constant from the axis to R/3 may be ex- 
plained if the assumption is made that the jet tip, and hence 
the critical zone, precesses about the duct axis. The preces- 
100} 


sion of jet tip has been observed by Agoston, Noon and 
Witherly (11). 
In Fig. 9, ¢ovs(7) has been plotted for a lean main stream. 
As in the case of the rich main stream, the observed equiva- 
lence ratio decreases slightly near the walls. ’ ; 


Pilotflame 
air jet It was not possible to maintain a flame when the sake con- 
sisted of pure fuel and the main stream of air only; however, 
it was possible to burn a “pilot flame” by having a fuel-air 
: jet with equivalence ratio gj = 2.8 and gn, = 0. The ob- 

i=, served equivalence ratios gos for the pilot flame are plotted 
wt = in Fig. 10. The luminous portion of the pilot flame did not 

fill the duct. 

met - ' When there was no fuel in the main stream, the flame burned 
with uniform equivalence ratio. Addition of fuel to the main 
stream (Gm. = 0.26 and gy», = 0.47) resulted in larger ap- 
equivalence ratio parent equivalence ratios for the center of the flame. Addi- 
tion of more fuel to the main stream causes an upward shift 
of the curve. Near the center of the duct, the flame burns 


air fuel 
jet 


60 


blow off velocity, fps 


Fig. 7 Blowoff velocity as a function of main stream and critical 
zone equivalence ratios. Duct size: 1j-in. standard pipe. 


l 2 0 
° stoichiometric air fuel jet 1. 
z 1.1 ° 0.9 stoichiometric air fuel jet 
0.97 oxygen jet air jet oxygen jet 
‘3 > 
0.8 
center wall 0.5 center : wal 
Fig. 8 Local equivalence ratio as a function of radius for a rich Fig. 9 Local equivalence ratio as a function of radius for a lean 
main stream (g = 1.16) and various jet compositions. Jet main stream (g = 0.92) and various jet compositions. Jet 
velocity: 480 fps. Mainstream velocity: 42fps velocity: 480 fps. Main stream velocity: 42 fps 
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Fig. 10 Local equivalence ratio in a pilot flame as a function of 


voms and radial distance. Jet velocity: 480 fps. Main stream 
velocity: 32 fps 


Equivalence Ratio for an Open Flame Stabilized by a 
Reverse Jet 


For the series of experiments involving an open flame, the 
main stream, which was an air-fuel mixture, emerged as a 
jet from a nozzle with an opening of 2 in. by 3 in. The flame- 
holding jet was located 16 in. vertically above the main stream 
nozzle. The open flame was stabilized between the outlets 
of the nozzle and the jet at a position which was determined 
by the jet pressure. It should be emphasized that the physical 
processes occurring in an open flame may differ from a flame 
stabilized in a duct. A new variable, the amount of en- 
trained air, is introduced. Since there is no wall to provide 
constraint, the flame spreads laterally. The jet composi- 
tion does not have a detectable influence on the combustion 
occurring at the tip of the flame. For an open flame the con- 
cept of a critical zone appears to lose its significance. 

Air from the surroundings was entrained by the main 
stream and by the jet, diluting the combustibles. In order 
to maintain a flame, the nozzle flow had to be considerably 
richer than stoichiometric. The results are summarized in 
Table 2. 


Conclusions 


The spectral intensity ratio method for determining equiva- 
lence ratio must use a calibration curve which is preferably 
obtained under operating conditions similar to those of the 
test flame. When applied to flames burning under conditions 
unlike those of the calibration flame, the absolute value for ¢ 
may be in error; however, relative values of ¢ are probably 
given correctly. 

The equivalence ratio is nearly independent of radial dis- 


Table 2 Observed equivalence ratio for an open flame 
stabilized by a reverse jet 

Velocity Equiva- 

at outlet Jet Average Input lence 
of pres- Fluid jet equiva- ratio in 

nozzle, sure, in velocity, lence critical 
fps psig jet? fps ratio zone 
42 14 Air 330 2.12 1.02 
39 12 Air 306 2.19 1.08 
38 12 Air 306 2.32 1.01 
40 14 Air 330 2.29 1.05 
40 14 Air 330 2.29 1.06 
40 9 Air 264 2.49 1.05 
53 10 Air 278 1.57 0.76 
53 10 No 286 1.56 0.79 
53 10 O» 254 1.57 0.79 
40 10 O» 254 2.18 1.06 
40 10 N2 286 1.05 
40 10 Air 278 2.18 1.04 
45 10 Air 278 1.68 0.83 
@ The inside diameter of the tubing used for the jet was 
zs in. 


tance for lean flames. For rich flames the equivalence ratio 
decreases near the duct wall. Jet composition not only in- 
fluences g near the axis of the duct but also affects the 
equivalence ratio at larger radii. 

The importance of the critical zone is illustrated by the ex- 
perimental result that there appears to be a single stability 
curve relating the blowoff velocity Voo to ¢., regardless of 
jet composition. 
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N THE past decade, as the result of a great deal of work, 

significant advances were made in all phases of rocketry. 
However, relatively little attention was given to the problem 
of nozzle configuration. This is easily understood, as there 
were many problems of greater significance which required 
attention. The conical nozzle, which had proved satisfactory 
in most respects, was used almost exclusively. The conical 
nozzles with a 15-deg half-angle (hereafter referred to as a 15- 
deg nozzle) has become almost a standard as it is a good 
compromise on the basis of weight, length and divergence 
loss. 

The literature sporadically has had suggestions on the 
possibilities of configurations other than conical for rocket 
nozzles. The Prandtl-Buseman design, also known as a con- 
tour or bell, has usually been suggested. However, most of 
the work in this field has been concerned with wind tunnels. 
In order to obtain uniform straight flow in the test section, 
these nozzles are designed with rather low expansion angles 
which result in long nozzles impractical for use in rockets. 
Clippinger (1)? did design a number of axisymmetric nozzles 
for air, y = 1.4, using initial expansion angles up to the allow- 
able maximum. Unfortunately, these results are not applica- 
ble to the rocket nozzle problem, and designs based on the 
rocket exhaust gases must be made. There has been some 
experimental work with one of Clippinger’s designs (2). 

Whereas the problem of nozzle configuration has been con- 
sidered, little design information has been presented. Re- 
cently, indicating a renewed interest in the subject, optimiza- 
tion of nozzle configuration has again been considered (3-5) 
and a few designs were presented (4,5). It may be presumed 
that this renewed interest is the result of the separate prob- 
lems of continuously trying to increase performance, and the 
large nozzles of large thrust engines. As it is becoming more 
difficult to increase performance by making major improve- 
ments, more attention must be given to minor details, such as 
nozzle configuration. 

Large thrust engines present a slightly different problem. 
As an example, the size of a rocket engine delivering 1,000,000 
Ib of thrust from one nozzle is considered. A chamber pres- 
sure of 300 psia and an exhaust gas of y = 1.22 are assumed. 


Presented at the ARS 13th Annual Meeting, New York, N. Y., 
Nov. 17-21, 1958. 
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A number of contour nozzles for an ideal gas with y = 1.22 are presented. The variation of gammi 
is considered. The design is based upon the assumption of a radial flow section in the initial diver- 
gent section which is turned to uniform axially directed flow at the exit. For each specified exit Mach 
number, about 95 per cent of the maximum expansion angles were used. The contours obtained were 
cut back in order to yield nozzles of various exit areas. The surface area of the divergent section, 
thrust coefficient and length of these cutback nozzles are compared with those properties of 15-de. 
half-angle conical nozzles with the same exit area. 
designs are also presented. These designs result in moderate savings in length, surface area an: 
heat flux over a 15-deg nozzle with the same thrust coefficient. Greater savings are possible at the 
expense of reducing the thrust coefficient. 


Simple heat flux calculations for some of the 


If the engine is to deliver this thrust at sea level, optimun 
expansion in a 15-deg conical nozzle requires that the diver- 
gent section be 91.5 in., or almost 8 ft, long. If the engine is 
to operate in a vacuum, an area ratio of 20 requires a nozzk 
length of 306 in., or 25.5 ft. Obviously, shorter nozzles woul 
be of interest. The other pertinent dimensions are listed in 
Table 1. The use of multiple nozzles is immediately sug- 
gested; however, these nozzles present new problems, and 
any thorough study must consider other configurations. 

This paper considers length, weight and thrust coefficient 
for a consistent set of contour nozzles for a gas with y = 1.22. 
These axisymmetric nozzles were designed on the same basis 
as wind tunnel nozzles; however, much larger expansion 
angles were used. The resulting contours were then cut off, 
resulting in nozzles which were practical for use in rockets. 
Some of the designs result in savings in length and weight over 
a 15-deg nozzle of the same thrust coefficient. Because a 
weight calculation requires a detailed design, the surface area 
of the expansion cone was used as a weight parameter. Heat 
flux calculations were made for several contour designs and 
are compared with conical nozzles of the same area ratio or 
thrust coefficient. Three contours with the same exit Mach 
number and different value of y are also presented. 


The determination of the contours presented herein is 
based upon the assumptions which have been previously used 
in the calculation of wind tunnel nozzles. The procedure, 
originally presented by Foelsch (6), assumes that downstream 
of the throat there exists a radial flow region resulting from a 
source flow. This flow is then turned in order to give uniform 
axial flow at the exit. Although most wind tunnel nozzles are 


Contours 


Table 1 
Parameter 
Pa, psia 
Pe, psia 
A ty in.? 
€ 

é Cr 


P. = 300 psia; y = 1.22 
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Table 2 Contour nozzle design data 
v(M), deg Om, deg deg Ri Calculated Theoretical? 
1.22 : 3 61.872 15.468 14.70 1.627 1.356 3.85 1.4247 1.4248 
122 4 85.372 21.343 21.00 1.975 1.612 5.24 1.4765 1.4781 
ee —— «4.6 96.38 24.095 22.90 2.416 1.850 7.03 1.5170 1.5194 
1.22 5 102.724 25.681 24.40 2.586 1.978 7.90 1.5309 1.5337 
1.22, 6 115.808 28.952 27.50 3.018 2.295 10.38 1.5545 1.5601 
1.22 a i 125.408 31.352 29.80 3.344 2.610 12.75 1.5692 1.5790 
a 8 133.800 33.450 31.75 3.579 2.884 14.87 1.5742 1.5895 
1.22 — 9 140.164 35.041 33.30 3.950 3.141 17.83 1.5739 1.5963 
4.6 92.93 23.233 22.10 2.353 1.797 6.71 1.5071 1.5090 
1.26 4.6 89.75 22.438 21.30 2.314 1.755 6.49 1.4985 1.5009 
See discussion under “‘Design Considerations.” 


two-dimensional, there have been some axisymmetric designs 
based on this method (7,8). 

Remarkably, there has been little experimental verification 
of this assumption in the literature. Apparently, the attain- 
ment of uniform flow was deemed proof of the assumption; 
whereas, nonuniform flows were discovered to be the result of 
other problems. Foster and Cowles (9) used an actual rocket 
motor and made pressure measurements on the centerline and 
wall to within 2 in. of a 1.5-in. throat. Oliver and Nagamatsu 
(10) measured the distribution across a section of a two- 
dimensional hypersonic wind tunnel. Both investigations 
showed good agreement with the radial flow assumption. 

Secondary proof might be gained from the agreement be- 
tween calculated and measured values of the thrust coefficient. 
There have been several experiments in which the variation of 
thrust coefficient with divergence angle has been measured. 
The agreement with the theoretical values of Tsien (11) is 
good in some instances (11,12) and only fair in others (13,14). 

The radial flow assumption, although it may not be exact, 
appears to be reasonable. The assumption of a distribution 
at, or near, the throat (1,3-5) does not appear to be any 
better. The necessity of further experimental evidence is 
obvious. The calculations based on the radial flow assump- 
tion are presented in the following section. 


Calculation of Contours 


The flow regions of a typical contour nozzle are pictured in 
Fig. 1. The contour in Region I, the entrance section, is 
arbitrarily drawn. Region II is the radial flow region. Its 
boundaries are the characteristic lines AB and BC. The 
flow is turned in Region III in order to give axially directed 
flow. The flow field is calculated by the method of charac- 
teristics in this region, the contour being essentially identical 
with the streamline through point B the inflection point. 
The flow is uniform and axial in Region IV. The nozzles, dis- 
cussed herein, were obtained by cutting off the aft end of the 


rather long nozzles obtained in the calculation, and therefore | 


they do not have uniform axial flow at the exit. 

The flow field along BC was determined by the method pre- 
viously outlined by Foelsch (6). Starting with the flow along 
BC and knowing that the flow along CD is uniform and axial, 


the flow in Region III is determined by the characteristics — 


method. The details of the calculation have been examined 
many times (1,7,8) and are not discussed here. A thorough 
discussion of the design techniques used here is currently 
under preparation.‘ Certain aspects of the problem are of 
interest and are mentioned briefly. 


* The programming and computation of the contours was done _ 


by Ying-Nien Yu of the Sandberg-Serrell Corporation, Pasadena, — 
Calif. ‘A Summary of Design Techniques for Axisymmetric 
Hypersonic Wind Tunnels’’ by Ying-Nien Yu is to be published. 
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The problem was originally programmed on an IBM 704 
for the calculation of axisymmetric nozzles for wind tunnel 
use. The program works off the centerline Mach number 
distribution, a procedure previously suggested for two-dimen. 
sional nozzles (15). The change in Mach number at C is 
rather abrupt. Therefore, the nozzle is lengthened by 2 per 
cent and a more gradual change made at C. A check of the 
M;.. contour, eliminating this added section, showed a negli- 
gible change in nozzle length over the region of interest. 

Since there was some concern about the functioning of the 
program with maximum expansion angles, angles about 95 
per cent of maximum were used. The data are presented in 
Table 2. Subsequent analysis showed that these fears were 
unnecessary, and the calculation of some contours at maxi- 
mum expansion angles are planned. The results of the 
present data are given in the following section. 


Design Considerations 


The design of a nozzle configuration must consider factors 
other than length. In general, length is probably the least 
important parameter. Among the other factors to be con- 
sidered are weight, thrust coefficient and heat flux. 

Since a 15-deg nozzle has only about a 1.7 per cent diver- 
gence loss, the possible increase in thrust coefficient is small. 
However, it would usually be desirable to keep the same thrust 
coefficient while making improvements in the other factors. 
A shorter length will not necessarily result in a lower weight- 
expansion section. However, as previously pointed out (3), 
a shorter length will probably result in a reduction of the 
weight of the attachment section and, in the case of swivel 
nozzles, the actuating mechanism. The heat flux may be 


smaller with a lower area; however, it must be investigated in 


ane. 


A c 7) 
REGION I — THROAT 
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Fig. 1 Contour nozzle flow regions — 
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order to be sure that the heat-transfer coefficient has not the throat diameter and tangent to the radial flow line extend- 
increased. ing to the inflection point. Recent work (13) indicates that 
a aur | this entry radius will probably be most efficient commensurate 

swith weight considerations. Smaller radii have been used at 

this laboratory. The distance from throat to inflection point 
The results of the calculation gave the contour downstream is then 
of the inflection point (B, Fig. 1). The contour between ‘ ; 
throat and inflection point was not computed, but drawn L; = [K + (Ri — 1 — K) cos 6;]/sin 6; (1) 
arbitrarily. This section has a radius of curvature equal to where K is the ratio of entry radius to throat radius and is 2 
for these nozzles. The values of LZ; used are given in Table ». 
The resulting contours are plotted in Fig. 2, in which the con- 


Contours and Length Comparison 


tour of a 15-deg nozzle is also plotted for comparison. Th: 
2 Me contours are designated by the full contour exit Mach numbe: 
Le The wall angles as a function of radius are plotted in Fig. 3. 
ca The foregoing method of obtaining the initial contour ma 
i 6 give a somewhat longer section than is necessary. In th 
2° respect, this design is conservative. A shorter nozzle will | 
2 | 59 obtained by assuming that the change in Mach number wit | 
K = theass length is constant from the throat to point A (Fig. 1) and equa 
—s to (dM/dL) at A. The method of characteristics may the 


be used to compute the contour in Region I (Fig. 1) (7,8) « 
the general method of this paper may be used (6). Althoug! 
wind tunnels use an arbitrary entrance section, their lon: 
- | radius of curvature results in a long entrance section, neces 
. sary to fit their special requirements for a flexible throat (16 
A contour for a particular ¢ is then obtained by cutting of! 
the aft end in order to obtain the desired area ratio. I 
_ order to facilitate the length comparison, the ratios of th: 
lengths of these nozzles to that of a 15-deg nozzle of the samy 
ay exit area are plotted in Fig. 4. Obvious savings in length ar 
possible; however, weight, thrust coefficient and heat transfer 
must also be considered. ace 


} | | 


A comparison of nozzle weights is difficult without a de- 
tailed design. As previously mentioned, some weight sav- 
ings will result from merely a shorter length, since the weight 
of the attachment section and, in the instance of swivel 
nozzles, the actuating mechanism will be reduced. The 
weight of a cooled nozzle constructed of tubing should be 
closely related to the surface area of the nozzle. For an 
uncooled nozzle, a true comparison requires more details. 
RaowS However, the surface area of the divergent section will bear 
some relation to weight for either type of nozzle. This sur- 
face area is presented in Fig. 5. The values up to the inflec- 
tion point are listed in Table 2. A comparison with a 15-deg 
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Fig. 2 Contours of nozzles, y = 1.22, varying M, 
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Fig. 3 Flow angles along wall, y = 1.22, varying M, 
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Fig. 5 Surface area of expansion cone 
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nozzle is made in Fig. 6. The surface area of a conical nozzle 


is 
S = 2K (11+ K)w — Ksinw] + (e—«)/sinw [2] 


Here w is the expansion angle, and ¢; is the area ratio at the 
tangency point and is related to the other variables by 


é; = [1 + K(1 — cos w)/]? [3] 


The first term in Equation [2] is the area to the tangency 
point, and the second term is the area of the conical section. 
I: appears that large reductions in nozzle weight are possible, 
but this will be at the expense of decreasing the thrust coeffi- 


cient. 


Thrust Coefficients _ 


The vacuum thrust coefficients obtained are plotted against 
e in Figs. 7 and 8. The values for the 15-deg nozzle are 
0.983 = [(1 + cos 15 deg)/2] of the theoretical one-dimen- 
s onal vacuum thrust coefficients. 

The vacuum thrust coefficients were calculated by subtract- 
ing the pressure integral along the wall from the value at the 
nozzle exit 


Pd 


‘The vacuum thrust coefficient at the exit of the nozzles, where 
uniform axial flow exists, can be related to the exit Mach 


number 
f.2 ‘ 
Cain yM.? + 1 x 
y-1 
— 
M, (: + ue) 


The data included the Mach number distribution along the 
wall, and a simple numerical integration was carried out in 
order to obtain Cp,. 
The thrust coefficient at any point can also be determined 
by integrating the flow conditions occurring along the plane of | 
constant L. The required equation is 
(yM? cos? 6 + 1) de 


This requires a rather tedious interpolation of the points in the 
Mach net which may be eliminated as follows. The net 
force across any two control surfaces meeting at the nozzle — 
wall will be zero. The integration can then be made along 
one of the lines in the net for which the data are available. 
The required equation is 


f cos R (cos + sin + RdR 
{1 + Cy — 1)/2) 


in which the correct sign and appropriate limits are taken. 

The flow-net data were used to calculate several vacuum- 
thrust coefficients for each of the nozzles for purposes of com- 
parison with the wall calculation. Significant differences 
were not expected. The net calculations gave higher results 
than the wall calculation, increasing from 0.1 per cent for the 
nozzle to about 4 per cent for the My nozzle. A small part 
of the difference was the result of the net spacing and, thus, 
the simple trapezoidal numerical integration used gives 
slightly higher results. In some instances values higher than 
theoretical were obtained. To date, these differences have 
not yet been resolved. 

The thrust coefficient at the inflection point calculated using 
{quation [4] may be compared with the theoretical value ex- 
pected from radial flow considerations. This comparison is 
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Fig. 7 Thrust coefficients, = 1.22, varying 
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Fig. 8 Thrust coefficients, M,.., varying y 7 


made in Table 2. Initially, the agreement is exceedingly 
good, the difference increasing as the initial angle increases. 
The calculated values are lower than the theoretical values. 
The theoretical values were obtained from Seifert and Crum 
(17), with the necessary correction for divergence loss being 
made (18). 
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LOCAL HEAT FLOW ( ¢/¢,) 


It is believed that these differences are the result of the net 
spacing. Future calculations will use a smaller net. For- 
tunately, the differences in those contours of interest are 
negligible. 

The thrust coefficient plot, Fig. 7, indicates that, for each 
area ratio, there will be a design from this consistent set which 
has the same thrust coefficient as a 15-deg nozzle. Reference 
to Figs. 4 and 5 shows that these nozzles will have shorter 
lengths and lower surface areas than their respective conical 
nozzles. The M,.. nozzle was designed after the others, as it 
appeared likely that it would cross the 15-deg thrust coefficient 
line at an area ratio of 20, a design currently of interest at this 
laboratory; it very nearly does fit this condition. 

Of greater interest would be a plot of surface area and 
length against thrust coefficient for each design in order to 
determine the minimum values of these parameters for each 
thrust coefficient. The minimums are rather flat, and the 
graphs are too crowded to be of practical use. For each 
thrust coefficient of interest, the surface area or length must 
be plotted against contour number (/,) to give the optimum 
contour. This has been done in Fig. 9. The minimums are 
rather flat and the MW, = 5 contour gives approximately mini- 
mum surface-area for expansion-area ratios from 6 to 28. 
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Fig. 9 Surface area comparison for constant thrust coefficient 


Further analysis shows that, although the contour obtained 
from the crossing point on the thrust coefficient plot is near 
minimum, further reduction in surface area is obtained with a 
contour of slightly greater exit area, shorter length and a 
larger divergence loss. This is only for vacuum thrust; if 
there is any back pressure the process must be repeated for 
those conditions. 


In order to obtain an idea of the extent of variation of he::t 
flux in the divergent section, calculations were made for sev- 
eral contour and conical nozzles at either a constant ¢€ of 20 
or a constant Cr, = 1.775. The simplified equations of Bariz 
(19) were used to calculate the heat transfer coefficient.> \ 
constant wall temperature was assumed. Calculations were 
made for the M4, Ms. and M; contours and 13-, 15- and 17- 
deg conical nozzles. The results are presented in Table 3 an 
Figs. 10 and 11. 

The local heat flux as a function of the length along tlie 
nozzle wall L’ is given in Fig. 10. Initially, the heat flux 
decreases as the expansion angle increases. Later, the curves 
for some of the contour nozzles cross that of the 15-dez 
nozzle. Although the total flux may be lower for some con- 
tours, the local flux at the exit cone may be higher than a 15- 
deg nozzle for some designs. This might affect the design of 
an uncooled nozzle. 

The total heat flux was obtained from an integration of the 
curves in Fig. 10. For the contour nozzles 


qt 


The total heat flux to the conical nozzles was calculated at this 
laboratory® 


Q Ine 
— =0.8+ - [9] 
quel: sin 


Inspection of Table 3 reveals that, in general, the contour 
nozzles will have lower Q than the conical nozzles. The sur- 
face areas and lengths will also be lower. The Q for the con- 
tour nozzles is plotted in Fig. 11. The M,. = 5 nozzle has 17 
per cent lower total heat flux than the 15-deg nozzle of the 
same thrust coefficient. The surface area is only 12 per cent 
less, indicating a reduction also in the average heat transfer 


coefficient. 


The effect of y was considered. An J; nozzle was de- 
signed with y = 1.24 and 1.26. The contours and other 
associated data are presented in the graphs and figures pre- 
viously mentioned. This is not the best way of studying the 
effect of y. Studying the variation of y in a fixed contour 


Table 3 Total heat flux 
—Cr, = 1.775~ —e = 
Nozzle € Cr, Q/qAt 
M, 16.4 12.4 1.801 14.7 
Mis 20.0 TI 1.775 11.1 
M; 22.8 10.3 1.759 9.6 
13 deg 18.8 13.8 1.782 14.1 
15 deg 20.3 12.4 1.774 12.3 
17 deg 22.4 11.4 1.765 Eee 


3 6 8 9 10 " 3 16 7 


5 Donald R. Bartz, this laboratory, performed the heat transfer 
calculations. 
6 Private communication from Donald R. Bartz, Aug. 1958. 
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would giye more useful results. However, this problem had 
not been programmed on the computer. 

Varying y while keeping the exit area and initial expansion 
angle constant will probably affect the contour only slightly. 
There is one limitation: The initial expansion angles cannot 
be too close to the theoretical maximum. However, if 
approximately 6; < 0.95 6,,, then the initial design conditions 
can be met, at least over the range of y and ¢€ of interest. The 
effect on thrust coefficient is unknown, but will probably only 
be that resulting from the change in y. - 


Comparison of Contour With Other Designs | 


The only contour designs in the literature with which those 
of this paper can be compared are those presented by Rao (4). 
In this reference, there are two contours for y = 1.23, € = 
19.36 and « = 4.973. There is insufficient information in 
this paper to enable a comparison with the smaller nozzle; 
however, the M4, data can be interpolated between y = 1.22 


and 1.24 in order to make a comparison with the larger 


nozzle. These data are presented in Table 4. 
Table 4 Nozzle comparison, y = 1.23 
JPL data*’—— 
Parameter I II Rao data (4) 

€ 19.4 15.3 19.36 

Cr, 1.7676 1.737 1.7676 

L 10.75 8.7 8.19 

S 64.5 48.1 48.1 


* Interpolation of M4.. contour. 
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If the nozzles are compared at the same thrust coefficient, 
the design presented by Rao has a shorter length and smaller 
surface area. If a nozzle of comparable surface area is 
selected, it has a lower exit area and thrust coefficient and 
slightly longer length. The length difference is the result of 
the shorter radius of curvature (K = 0.45) used by Rao. 
The difference in thrust coefficients is 1.75 per cent. The 
thrust coefficient in Rao’s design may easily be that much 
lower in actual use because of the small radius of curvature 
used at the throat (13). Based on the actual vacuum thrust 
coefficient to be expected, it is believed that both methods of 
design will result in essentially the same savings in weight and 
length. 


Conclusions 


The information presented in this paper has been abstracted 
from a Jet Propulsion Laboratory Memorandum (20). Addi- 
tional points discussed in the memorandum are the pressure 
ratio and Mach number distribution along the wall. The 
pressure ratios are larger than the theoretical one-dimensional 
pressure ratios allowing a slightly larger expansion ratio for a 
particular atmospheric pressure. The variation of Mach 
number with ¢€ is approximately the same for all the contours. 
The memorandum also contains, in tabular form, all the 
design information discussed herein. 

This paper has presented a number of bell-shaped nozzle 
contours designed on reasonable assumptions about the initial 
flow conditions in the nozzle. Some of the designs result in 
savings in length, weight and heat flux over a 15-deg nozzle with 
the same thrust coefficient. It is believed that the design is 
conservative and that these nozzles will exhibit the perform- 
ance indicated. 


Nomenclature 

A 
A = area 
Cry = thrust coefficient 
K = radius of curvature at throat/R, 
L = length along centerline/R, 
= length along wall/R; 
M = Mach number 
= pressure 
Q = total heat flux, Btu/(time) ; 
q = local heat flux, Btu/(area)(time) 7 
R = radius (from centerline)/R, 7 
S = surface area of expansion section/A; : 
¥ = specific heat ratio 
€ = expansion area ratio 
0 = flow angle with centerline 
v(M) = Prandtl-Meyer angle 
w = expansion angle of conical nozzle a 
Subscripts 
€ = chamber 
e = exit 
m = maximum 
t = throat 
v = vacuum 
w = wall 
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Relativistic Rocket Flight 
‘With Advanced Systems 


HEMICAL rockets, of producing high 

thrust-to-weight ratios, are essentially mass limited. 
That is, no standard chemical reaction (as opposed to a free 
radical reaction) is sufficiently energetic to yield specific im- 
pulses of over a few hundred seconds. Asa result, the ratio of 
initial to burnout rest mass becomes very large. To surmount 
this limitation, advanced propulsion systems attempt to utilize 
more energetic sources of energy. These are listed in Table 1, 
with matter annihilation as a point of reference. Nuclear 
fission is the most energetic source which we currently know 
how to utilize on a large scale. 

Generally speaking, however, it is not possible to use the 
more energetic sources in as direct a manner as chemical pro- 
pellants in a rocket. In the latter case, the products of the 
chemical reaction are themselves exhausted directly. In the 
former, reaction energy is usually transferred from a separate 
energy source (such as a reactor) to the jet. Although this 
adds some flexibility, the greater bulk involved severely limits 
the acceleration capabilities of such systems. This paper de- 
velops a general approach to a wide variety of propulsion sys- 
tems, clearly indicating the performance limitations imposed 


by both fuel energetics and specific power. oe ae 


In developing a general theory of rocket propulsion, it is 
convenient to first separate fuel energetics from considerations 
of powerplant weight. For this purpose let us imagine a 
rocket engine fixed on a test stand. (The test stand itself will 


Powerplant Performance 
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% This paper presents a general relativistic theory of rocket propulsion, applicable to either adiabatic 
or diabatic jets, and to several methods for handling fuel residual. 


Emphasis is given to the per- 


formance limitations imposed by both fuel energetics and power supply weight. The discussion is 
limited to constant-thrust vehicles. 


be assumed fixed in an inertial coordinate system.) Define: 
my; = total fuel consumption per unit time 


In a heterogeneous reactor, for example, m; would represent 
the fresh fuel element mass divided by the length of time be- 
tween shutdowns for recharging. 
€ = fraction of my converted to energy 

In accordance with Table 1 € could never exceed about 107? in 
a fission reactor. 

= mass of inert? working fluid ejected per unit time in the 

jet 


Table 1 Primary energy sources 


Fraction of 
Energy release, annihilation 
Type erg/gm energy* 
Chemical 
(hydrogen-oxygen ) 
Free radical 
(H + H > 2.2 101 2.45 X 10-° 


Nuclear fission, with 

complete “burnup’’? 7.4 7.9 xX 
Nuclear fusion 3.6 X 10% £0 
Matter annihilation 9.0 < 10% 1.0 


@ Clearly, this is also the fraction of the fuel mass con- 


verted to energy. 
>’ Only 1 per cent burnup is not uncommon in current 


reactors. 


2 Inert is here used in the sense that the working fluid does not 


enter into the energy producing reaction. 
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In a nuclear rocket, m; would thus represent the mass flow of 
hydrogen, ammonia or other working fluid through the nozzle. 


n = fraction of unconverted fuel mass (residue) ejected with 
the working fluid, at the jet velocity 


Thus 7 would include any fission products or unconverted = 


uranium caught up in the stream of working fluid. For an | 
ordinary chemical rocket, 7 = 1, since the jet is made up en- 


tirely of reaction products. 


jet velocity relative to the test stand 


c = light velocity 

a= m,;/my 

w = fraction of the power developed which is wasted, does not 
add to the thrust 

g = acceleration of gravity 

Now allowing relativistic jet velocities, conservation of 


mass and energy can be written as (1) of. *t 
mjc? (l—e)nme? 
+ 
V1 — (v/c)? (v/c)? 
(1 — €)(1 — )myc? + wemyc? 


myc? ++ m;c? 
[1a] 


This can be solved generally to yield 


(:) -| | [1b] 
c (1 + a) — (1 — e)(1l — 7) — we 


with the resulting specific impulse* 


g l+a 


[1c] 


= — w) 


We now wish to consider three special cases. For the first, 
= w= O0and7 = 1. There is no separate working fluid; 
instead, the reaction product is used in the jet. An ideal 
chemical rocket would be characterized by this choice of 


parameters. From Equations [1] a 


This formula has been previously given by Sanger (2). The 
specific impulse appears as 


I = (c/g)V 2e — & 


Therefore the maximum attainable specific impulse using pure 
hydrogen and oxygen would be about 500 sec. An ideal rocket 
of this sort will be referred to as adiabatic, since there is 
neither a separate powerplant (such as a reactor) to add 
energy to the jet stream, nor radiation losses. 

The second case will fit either a heat exchanger nuclear 
rocket or an ionic propulsion device. For this case a separate 
working fluid is used, so that » = 0. Now the jet itself is 
nonadiabatic even though w = 0. We shall refer to propulsion 
units with separate working fluids as ‘‘diabatic’”’ rockets. 
Again from Equations [1] 


(v/c) = V2 — 


[3] 


> 
_¢ V2a6 + 6? 

9 (l+a) 


where the specific impulse takes into account consumption of 
fuel as well as working fluid. It should be noted that as 
a— ©,]-—v/g. However, for relativistic jet velocities, this 
simplified form in which we neglect fuel consumption is quite 
different. For example, if we take «, = 10-4 and v/c = 1/2 
then J is approximately equal to (c/2g) X 107%. 


3’ Numbers in parentheses indicate References at end of paper. 
‘ The specific impulse as here defined takes into account con- 
sumption of fuel as well as working fluid. In the discussion, we 
will justify our preference for this definition. 
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The quantity a may be varied to maximize J. The par- 
ticular value of aopt = 1 — & (which for ¢, << 1 means equal 
consumption of fuel and working fluid), and 


v 


These formulas illustrate clearly the limitations imposed by 
fuel energetics. If we try to increase (v/c) indefinitely, J will 
actually decrease due to increased fuel consumption. For a 
heat exchanger nuclear rocket €, is about 10~‘ or less. There- 
fore (v/c)opt iS approximately equal to 107, and Jo: is ap- 
proximately equal to 105 sec. 

The third case is the photon rocket. Applicable formulas 
can be obtained from Equations [4] by merely allowing 
Then (v/c) ~ 1 and 


| 
T photon = : w) = 
9 g 


For = 10~‘, the photon rocket would yield an J photon only 
about 10~? of Zope (Eq. [5]). The photon rocket will give 
optimum performance only for ¢; = 1 (100 per cent efficient 
mass annihilation). 

In passing, one might comment that this general approach 
can also be applied to solar powered vehicles, for example, solar 
heating of a gaseous working fluid such as hydrogen. It is 
only necessary to replace m,c? by the collected solar power, and 
¢, by the fraction of this power which appears in the jet. The 
factor 7 of course would equal zero. 

In this section we have illustrated the test stand limitations 
imposed by fuel energetics. Next let us turn to the per- 
formance of a rocket vehicle in gravity-free space. This will 
bring into focus the additional limitations imposed by power- 
plant weight as well as some other interesting points. 


Straight Line Burnout Velocities in Gravity-Free 
Space 
For a diabatic photon rocket, Singer (2) gives the formula © 


= (7) 
c 1+ 
where 
U = burnout velocity in relation to an inertial coordinate 
system at rest relative to the initial position of the 
rocket 
vy = ratio of final to initial ship rest masses in the same co- 
ordinate system 
€ = the fraction of fuel mass converted to energy _ 
On the other hand Tsien (3) derives a form 7 
U 1 — 
e i+ 


which can apply to either an adiabatic or a diabatic rocket. 
However, for a diabatic rocket, setting v = c does not deter- 
mine ¢ (see Eq. [4]). Therefore, Equations [7 and 8] do not 
agree except for the special case of €, = 1. 

To resolve this apparent discrepancy and elucidate the 
problem, we shall derive a general result, including both 
Tsien’s and Singer’s as special cases. Let us define the follow- 
ing symbols: 
instantaneous ship velocity, relative to a coordinate 


Uo = 
system at rest with respect to the initial position of 
the rocket 

m = instantaneous ship rest mass in the same coordinate 
system 

u. = rocket exhaust velocity relative to the same coordinate 
system 
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at 
5) 
a 


dm, = rest mass of the mass ejected from the rocket during a 
small time interval in the same coordinate system 
dm = rest mass of the mass of unconverted fuel (residue) pro- 


duced per unit time 
 g = fraction of dm dumped overboard at zero velocity (in 
the U direction) relative to the sete during the same 


time interval 


f+ns1 
Then from these definitions 


edm 


dmy = —dm, — 
(1 — e€) fi 
=" 
Therefore 
= d 


+ 1 -©) 


Conservation of mass and energy® requires 


a( moc? ) 
V1 = (w/c)? V1 = (w/e)? 

fce*dm 

= 


conservation of momentum requires® 


(1 — €)V 1 — (u/c)? 


d ( MoUo ) ae, udm bs 
— (o/c)? V1 - (u:/c)? 
fucdm weugdm hel 
In addition we have the relativistic formula 


where v represents the exhaust velocity relative to the ship.’ 
The formulas derived previously for v/c can now be applied’ 
so that for a diabatic jet way [4] would yield 


c 


(We are, of course, assuming a constant thrust throughout 
this discussion.) Equation [9] can first be substituted into 
Equations [10 and 11] to eliminate m, and m. The factor 
Vi — (u/c)? can easily be eliminated between the resulting 
two equations, and (u/c) itself by application of Equation 
[12]. The final formulation appears as 


5In the coordinate system at rest with respect to the initial 
position of the rocket, and in gravity-free space. 

6 Actually only one of Equations [10 and 11] need be used. 
However it seems simpler to carry out the derivation as indicated. 

7y is now measured in an inertial system instantaneously mov- 
ing with the ship. 


Soil 2 


P = 
99) 


_ me 
mo (1 €) 
(2) | (2) 
E — (u/c)? (*) | 
Integration yields 


n=1, f=0, w=0; 6=c/v adiabatic rocket 


Tsien 
v 


diabatic rocket 


Here I is the specific impulse, including fuel consumption, «is 

previously defined. For the particular case of a diabatic 
photon rocket, J = (c/g)e, and Equations [14] reduce to 
result. 

From the table of special cases, it is clear that the difference 
(with w = 0) between Tsien’s and Siinger’s results for diabatic 
rockets lies in the treatment of fuel residual. Tsien keeps it 
aboard, whereas Singer has it discharged (at zero velocit. 
relative to the ship) as formed. The latter course is to bw 
desired, but may not always be possible. 


Discussion 


When fuel residual is discharged along the way, v may be 
written as 


ke*em; + P 


7 +P) + (mp + mtg 


# 


payload (Earth weight)® 
powerplant weight per unit of useful power output 
time interval during which constant thrust is applied 


Under some circumstances, such as delivery of a powerplant 
(ke*eym;) might correctly be included in the payload. In 
any event, if one can take an unlimited time to attain a pre- 
scribed velocity increment, the powerplant weight can be 
made vanishingly small. Under these conditions, vy = P/G, 
where G is the gross initial weight. The specific impulse / 
should then be maximized to achieve a maximum pv. This 
explains our preference for a definition of specific impulse 
including fuel consumption; it represents a quantity which at 
least under these limiting conditions should be made as large 
as possible. If one were to take instead the usual definition 

= v/g, then J would appear in the exponent of Tsien’s case. 
However, the specific impulse with this definition should not 
be maximized, since one must also consider the fuel residual 
carried through to burnout. 

Let us now turn finally to the problem of attaining a pre- 
scribed velocity increment in a given finite time interval. 
Taking again the special case of a diabatic rocket with con- 
tinuous discharge of spent fuel, the appropriate formulas 
drawn from the previous text appear as 


+ a)/2V2ae + at 

1+ U/e 


(Z) - 
G 


1 + (P/kc*em;) 


From Equations [16], with «,, t, g, c, k and U/c given, a@ will 
determine (P/G). There will be one value of @ which will 
maximize (P/G). These values have been computed numeri- 
cally for a range of the variable t,, and the appropriate results 
are presented in Figs. 1 and 2. 


8 To be more consistent, one could of course express the right- 
hand side of Equation [15] entirely in terms of masses. 
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Fuel residual discharged as formed 
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| 
498 107 108 109 10'° 
tp (sec) 
: - Fig. 1 Optimum specific impulse and exhaust velocities as functions of mission duration 
60 


As ty ~ ©, a approaches 1 — & corresponding to maximiz- 
ing the specific impulse (Eq. [5]). The role played by « is then : 
dominant. For lower values of t, the effect of powerplant 
weight becomes increasingly important. When a becomes 
large, J approaches v/g, and the results begin to behave as if 
the fuel were massless. An approximate analysis made on the 
basis of a massless fuel and nonrelativistic velocities® has been 
given by Langmuir (4). Specific powerplant weight is then 
the controlling factor. The results are cut off when é be- 
comes so small that P/G = 0. The required velocity incre- 
ment could not be attained, with the prescribed parameters, 
in any shorter time interval. When P = 0, kc*e, is about atyg; 
i.e., powerplant weight equals propellant weight. 

= 


@ 
Conclusions 


We have given an analysis sufficiently general to cover: 


1 Fuel mass limitations. 
2 Specific power limitations. 


In our numerical example, U/c is also much less than 1. 
However Equations [16] would hold even if this were not true. 
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21074, 5,=10°5, k= 10°? sec/ft #10 1b/kilowatt 
Fuel residual discharged as formed al 
Initial acceleration =0,5x107> "g" 
107 108 109 
tp (sec) 


Fig. 2 Optimum value of P/G and a vs. mission duration 


3 Relativistic velocities. 
4 Various methods for handling fuel residual. 

One specific example has been carried through, illustrating 
the shift in importance from fuel energetics to powerplant 
weight as the allowable mission time shrinks. 

Long flight times and high efficiency lightweight energy 
conversion equipment both tend toward increasing the rela- 
tive role played by fuel mass. A photon rocket cannot be the 
optimum choice unless we have a complete mass annihilation 
energy source and massless conversion equipment. 


i’ 
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ARLY investigations of co-planar transfer trajectories be- 

tween co-planar orbits using instantaneous velocity in- 

crements were conducted by Lawden (1-3)? and indicated 
that: 

1 Approximation methods are required for the general 
optimal solution. 

2 The optimal number of impulse velocity changes is 
generally not known. 

3 For a two-impulse transfer there exists a graphical 
method for choosing a particular trajectory (2). 

4 For a two-impulse transfer between circular orbits, the 
optimal method is Hohmann’s tangential ellipse (3). 

5 For transfer between orbits of similar size and shape, 
the best tangential ellipse (two-impulse) gives a good ap- 
proximation to the optimal solution. 

The success of a cotangential type of transfer in providing 
an optimal or near optimal transfer between certain types of 
co-planar orbits with respect to characteristic velocity re- 
quirement has motivated the present investigation of two- 
impulse cotangential transfers between nonintersecting cir- 
cular and elliptic co-planar orbits. Between circular and 
elliptic co-planar orbits there are, in general, an infinite num- 
ber of possible cotangential transfer ellipses. In the present 
analysis all such transfers are examined as a function of the 
location of the injection into (or departure from) an outer 
elliptic orbit. It is found generally that injection (or de- 
parture) at outer elliptic orbit apogee is the most economical 
and injection (or departure) at perigee the least economical. 
The most economical transfer to an outer elliptic orbit of 
particular energy (or period) is found to occur when the 
eccentricity of the orbit is such that perigee is located on the 
circular orbit. Minor modifications of the analysis can be 
made indicating reversal of the choice of apogee and perigee 
injection or departure for the case of transfer between inner 
elliptic and outer circular co-planar orbits. The present 
analysis makes use of the simplified two-body assumptions. 

Subsequent to the origina] submission of the present paper, 
personal communication with Silber brought to the au- 
thor’s attention the existence of two reports (6,7) not availa- 
ble in the open literature but related to the present paper. 
In the first report (6) the authors have made the interesting 
discovery that for co-planar transfers between circular orbits 
whose radii are at a ratio larger than 11.94, there exist three- 
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Co- Eliptic Orbits 


Vo 2 Veo + AV, a 
-( Va ) 


Convair-Astronautics 
San Diego, Calif. 


The general energy requirement measured in terms of a characteristic velocity for the ascent from 
7 .. - an initial circular parking orbit to any outer co-planar elliptic orbit by a two-instantaneous-impuls-, 

Se cotangential transfer is determined as a function of the location of the second impulse on the fin} 

: orbit. Of the infinity of such transfers to any particular final orbit, it is demonstrated general! y 

_ that injection into the final orbit apogee is the most economical whereas injection at final orbit 

-. perigee the least economical. For the special case of final orbit perigee located on the initial parl. - 

_ ing orbit, the most economical transfer requires a single impulse applied at perigee. 


impulse, “‘bi-elliptic’’ transfers more economical than the 
Hohmann two-impulse transfer between such orbits. In the 
second report (7) the authors investigate the special cases cf 
perigee and apogee transfer from an inner circular into an 
outer co-planar elliptic orbit and determine that apogee in- 
jection is the more economical. 

More recently, Munick (8) has demonstrated that the two- 
impulse cotangential transfer from inner elliptic orbit perigee 
to an outer co-planar circular orbit is the a two-impulec 
transfer from perigee. 


General Energy Requirement 


Derivation 
First impulse magnitude 
4 


By the conservation of angular momentum a 7 
[1] 


c? = V*r? cos? 8 = constant 


For the transfer ellipse, we have ? 
2 1 


= [3] 
bo = (Vo/Veo)? 
atr = ro. Therefore at P’ in Fig, 1 
2 2 b 
= = cos? B [4] 
if To To 


and generally 
_ 2(r/ro(r/ro — 1) cos* B 
(r/ro)? cos? B — 1 [5] 


At the point of tangency F, let r = rr, o = rr/roand B = Br. 
Now, since 
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se 
L. RIDER! 
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‘ 
> 
= 
his ~ ar 
2 — bo 
» 
2 


-1 [7] 
V a? cos? Br — 1 


then 
TRANSFER 


where Br is to be determined. Fd 
Consider the characteristics of the final orbit at the point of © jf 
tangency F 7 
c? = Vr*rr* cos? Br . 
r= \ 
a \ 
\ 
x 
Substitution gives q ' 
6 
cos*Bp (9] Fig. 1 General configuration 
However, the relation between the kinematics and geometry ‘ 
of an elliptic orbit is such that FINAL ORBITAL PERIOD = 24 HOURS ; 
e€=04 
which provides the attitude at any point on the final orbit as . iialicniataalin ia 
a function of the final orbit characteristics and the point of in- A eee Oe AT APOGEE OF 
terest r = rp. 
To satisfy the condition of tangency at F it is required that 
8r in Equation [7] be equal to Br in Equation [11]. There- 
fore, substituting Equation [11] in [7] we have 
Veo o(1 — — 2 *) £8 
a\o a ae 
Fig. 2 Variation of characteristic velocity with second impulse 
Determination of second impulse magnitude location ss 
The second impulse magnitude AV, will be that velocity . - 
added to the velocity in the transfer ellipse at / such that the Ch ne : a 
terist locity 
total velocity equals that in the final orbit at F. a lnaninaiaiieines 
The velocity in the final orbit at F is given by a The total velocity requirement is given by 
2 1 Vr 2 To Vcu AVo AVr 
2 a Veo a Veo Veo Veo 
The velocity in the transfer ellipse at F is ws ‘ xe ie [ 
To {2 ro | 2 
Substituting [5 and 11) results in 
7% It is of interest to note that for a final orbit which is circular, 


Xo — (*) °) Oand 


2 


| 26 (: - 4) 


To _ 


Co a 
| Xo—1)](/n\ (2. 1 which is exactly the optimal solution for the two-impulse 
| (: *) transfer between circular orbits (Hohmann transfer). 
The variation of characteristic velocity with second im- 
eT. (*) ( - =) pulse location on the final orbit is shown in Fig. 2 for various 
_ \ a/\o a eccentricities and final orbital period of 24 hr. From the 


= 


AV in 
~ 
ae A' ~ 
FINAL 
ORBIT 
y 
ge 
: 
Thus, secon 
Vex 
1 8 
V 
\ 


figure it can be seen that as much as a 24 per cent savings in 
characteristic velocity (depending upon eccentricity) may be 
achieved by apogee rather than perigee injection. (The initial 
circular orbita] altitude is taken as 100 statute miles.) 

Fig. 3 illustrates the variation of characteristic velocity with 
final orbital period for perigee and apogee injections into final 
orbit. In this figure ¢ in the range 2 ie 


refers to an injection at final orbit iy aaa 


9 
-1) 
To To 
refers to injection at final orbit apogee. For o = a/ro 


Veu/Veo is the Hohmann characteristic velocity. 


Behavior of the Function Vc, /V-) Over the Possible Range 
of the Independent Variable 


Fig. 2 illustrates that for a final orbital period of 24 hr, 
apogee injection is more economical] than perigee injection. 
Since the curves were calculated at finite intervals, it is 
possible that at some intermediate value of the independent 
variable o a minimum over the closed interval may occur; 
i.e., application of the second impulse at some point of the 
final orbit between apogee and perigee may be more eco- 
nomical than for second impulse at final orbit perigee. 

Additionally, even if Fig. 2 represents the true variation of 
characteristic velocity with location of the second impulse on 
the final orbit, the character of the variation has been shown 
only for a final orbit period of 24 hr. It is true generally? 

That the previous questions are pertinent and not of aca- 
demic interest only is illustrated by the solution to the Hoh- 
mann problem which provides the minimum energy ascent 
from an initial circular orbit to an outer co-planar circular 
orbit. The results of the Hohmann analysis (5) indicate that 
the two-impulse incremental velocity requirement is greatest 
for a final orbit located at approximately 15.6 initial orbit radii 
from Earth center; i.e., the energy requirement to a final 
circular orbit located at any distance greater or less than 15.6 
initial orbit radii is less than that required at 15.6 initial orbit 
radii. 

The generality of minimum energy ascent occurring for 
second impulse at final orbit apogee for any particular final 
orbit will be demonstrated by showing: 

1 The maximum and minimum values of Vex/V-o occur at 


2 


05 


o4 


RIOD 


03 


02 


24 HRs 
10 12 
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Fig. 3 Variation of characteristic velocity with orbital — 


the end points of the range of the independent variable o for 
any final orbit configuration where a 


ew — 119) 
To To 


i.e., re ranges between values at perigee and at apogee. 


2 Second impulse at final orbit apogee is always equally 


or more economical than second impulse at final orbit perigee. 


Statement 1 is accomplished by demonstrating that 1:0 


- maximum or minimum of Veu/V<o exists except at the end 
points of the of o, for any final orbit 


E=1-e70 
m= V/2E 
A = n/a #0 
n V2A 
B=E+A? 
and let 
y? = (o — 1)/(Bo — 2A) 
2? = (2 — Ao)/o 
The equation for Vex/Veo then becomes 
Ven /Ven = my — nyz +2 


and 


do 


= = (m — nz)jy + (1 — ny) (23) 


Since y? and z? are given functions of o, we can eliminate o 
and find 


> 
2 — (22 + A) 
m2 — 


Now, the problem of finding the o values which correspond 
to maximum and minimum values of Vcen/V.o where y(o) 
and z(c) can also be considered a problem of constrained 
maxima and minima. That is, it can be formulated as: 

Find the extremal values of 


y= 


Veu 
= my — nyz+z-—1 
Veo y y 


subject to the condition 


ym? — + 224+ A — 2 [25] 


_ The solution is supplied by the following equations (4) 


oy oy 
oz oz 
where is the undetermined Lagrangian multiplier. 
The Equations [25 and 26] are equivalent to 


" (m — nz) + 2Ay(m? — n2z2) = 0 


(1 — ny) + 2A2(1 — = 0 
ym? — n’z*) + 227+ A—2=0 
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which provides the solution 


= | [28] 


+A?+ — A(2 — A)] 


the o values giving extremal values of Vcu/V ce. 

There are two nondegenerate cases which can be interpreted 
by realizing that for final orbit perigee on or outside the initial 
parking orbit we must have EF > A(2 — A). 

Now, for E > A(2 — A a i 


2A 
E+A?+ VBE[E — A(2 — A)] 
2A 
29 
E + A? — VE[E — A(2— A)] A [29 
Algebraic reduction shows this statement to be true if 
+ (1 + €) > (3 + [30] 


Since we consider only those situations where 0 < € < 1, the 
statement is satisfied. Consequently, for final orbit perigee 
outside the initial circular parking orbit, no maxima or minima 
of occur over the open interval! 


1+e 
o = 04 
A A 


For the limiting case where final orbit perigee is located on 
the parking orbit, we substitute H = A(2 — 4A), ice. 


in Equation [17] and find Oo 


Vi (31). 


Thus for final orbit perigee on the initial parking orbit, 
the value of Vex/Vco is independent of ¢, and all possible o 
values between apogee and perigee are equally economical (in 
fact the most economical to a final orbit of particular period). 
However, in this particular case there are only two possible 
values—those corresponding to apogee and perigee of the 
final orbit. This follows since the only possible elliptic co- 
tangential transfers between the parking orbit and final 
orbit coincide with these orbits and either first or second im- 
pulse is zero. Both possibilities reduce to a single impulsive 
transfer with impulse applied at the final orbit perigee. 

To demonstrate Statement 2 that apogee injection is 
equally or more economical than perigee injection, we make 
use of the expression for characteristic velocity when second 
impulse is at final orbit apogee or perigee. The general ex- 
pression reduces to 


_ ( 1 
Veo Co 


It is required to show that 


o+1 


20 


That is, show that 


1 
(1 Aop _15 
op op + 1 Oop 
+ 
Oa +1 CA 


1.0 
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(ON 
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Fig. 4 A = {V(1 — &)[(1 + A)? — &][(1 — A)? — &]} 
{(1 + A)(1 — A)? + [A2(1 — A) — 2]e® + (1 + Ade} 


where 
op = —— 
A 
l+e 
Cs, = 
A 


_ Algebraic manipulation provides the following necessary 
and sufficient condition 


A= {V0 — + A? — — A)? - 
{(1 + Al — A)? + [A%1 — A) — Qe? + (1 + Ade} 
> 0 [33] 


subject to the side condition that e < 1 — A (ie., that final 
orbit perigee is at a radial distance equal or greater than 
initial orbit radial distance). 

Graphical analysis of the A function indicates that all 
possible combinations of € and A consistent with the side 
condition provide values of A 2 0. 

A contour representation of the three-dimensional surface 


defined by Equation [33] is shown in Fig. ne 


Nomenclature 

r = radial distance from Earth center to point of interest 

V = scalar velocity 

V. = circular velocity a 

a = semimajor axis of the final orbit ; 

= semimajor axis of the transfer orbit 

€ = eccentricity of the final orbit 

= 

AV = impulse magnitude 

Venu = characteristic velocity = AV) + AVr 

B = angle between tangent to the orbit and the local hori- 
zontal 

b = (V/V-)? square of ratio of velocity to local circular 
velocity 

Subscripts 

0 = conditions at the initial circular orbit 


) 
| 
A=0 
0 0.2 0.4 0.6 
i 
+ (; CH > Vou ah” 
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final orbits 
= apogee of final orbit 
= perigee of final orbit 
= the transferellipse 


Sum 


conditions at the point of tangency between tcanster and 
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_ For orbits or trajectories that carry large area lightweight structures outside of protected regions 
: ot space defined by planetary magnetic fields, the damaging effects of solar plasma streams are 
shown to be potentially serious and may result in low durability or high payload penalties. Esti- 
mates given depend directly on plasma stream densities and velocities which are only very ap- 
_ proximately known. However a reduction by one order of magnitude or more in widely held current 

estimates of the properties of the solar plasma would still imply a considerable effect on particular 


 ilarly affected. 


> 


ECENTLY there has been increasing interest in the 
possibility of erecting structures of large area on or near 
the surface of the moon or planets or in interplanetary space. 
Large area structures are of potential utility as marker bal- 
loons, as passive components in communications systems, as 
sampling devices and as focusing collectors for power produc- 
tion from solar radiation. They have often been suggested as 
means for production of thrust from light pressure. Because 
of the great cost per pound of payload, all such devices are 
necessarily envisioned as extremely lightweight and extremely 
thin. The possibility of using very fragile structures derives, 
of course, from the very low forces of acceleration to which 
they are to be subjected during the performance of their mis- 
sion. 

It is of interest to examine limitations on the use of such 
structures imposed by certain aspects of the interplanetary 
environment. Although there are still many gaps in our 
knowledge of the environment of space, it is reasonably cer- 
tain that, outside of limited regions protected by the presence 
of planetary magnetic fields, space is permeated by an ionized 
plasma whose effect on very thin structural members could be 
considerable. In some cases what would ordinarily be thought 
of as a surface effect from such a plasma might become, in 
actuality, a bulk volume effect to the detriment of important 
physical properties. Any such effects are additional to those 
produced by solar electromagnetic radiations, exposure to 
high vacuum, and other more frequently considered aspects 
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structural materials. Substances used as thin coatings on more massive structures would be sim- 


of the environment of space. The ionized plasma presumed 
to be issuing from the sun, especially during periods of high 
solar activity, can also produce mechanical loadings far in 
excess of light pressure effects. We shall consider some of 
these effects of solar plasmas after a brief review of current 
ideas on the nature of corpuscular emission from the sun and 
the particle content of interplanetary space. 


Quiescent Solar Plasma and Solar Streams 


There exist a variety of observations which lead to the con- 
clusion that interplanetary space out to somewhat beyond | 
astronomical unit (A.U.) (1 A.U. = 1.5 X 10'% cm = radius 
of Earth’s orbit) of the sun is occupied by a tenuous highly 
ionized gas. Details as to the composition, temperature, and 
the possible presence of non-Maxwellian components are not 
completely certain, and such knowledge must presumably 
await the results of space probe experiments, some of which 
are now underway. Strong evidence for the presence of an 
outward streaming component from the sun has, however, 
been presented by Biermann (1,2)? in early papers and more 
recently (3). He shows that the accelerations and radial 
alignment observed in type I comet tails containing COt, N2*, 
etc., cannot be explained in terms of solar radiation (light) 
pressure. These accelerations ordinarily appear to be 100 to 
200 times as large as solar gravity but may reach 1000 times 
solar gravity under very active conditions on the sun. A 
streaming solar hydrogen plasma would produce the observed 
accelerations because of electric fields set up by the greater 


* Numbers in parentheses indicate References at end of paper. 
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stopping power of the cometary material for electrons relative 
to the protons. Neutral ionized streams having velocities of 
about 1000 km/sec and densities of the order of 10% to 104 
particles/em$ at the orbit of Earth were suggested to account 
for these observations. Biermann points out that the highest 
accelerations can be related to interactions with individual 
solar corpuscular beams of still higher density presumably 
connected with discrete geomagnetic storms. On occasions of 
high solar activity, the suggested density of the streams (4) 
may reach 10° protons and electrons per cm’ at 1 A.U., and 
velocities might be upwards of 2000 km/sec. Similar values 
nay be obtained from observations on changes in solar radio 
noise (5). 

For a quiet sun, Behr and Siedentopf (6) have examined the 
polarized component of the zodiacal light and find an electron 
density at 1 A.U. of 700 particles/cm*. Electrical neutrality 
would then require an equal number density of protons. Ina 
discussion of Biermann’s early papers, Kiepenheuer (7) 
cerives a lower limit for the particle density of about 6 
j articles/em’ assuming a velocity of 1000 km/sec for the solar 
plasma. 

Recently, Kuiper has suggested (8) that a reasonable 
average density for quiet conditions might be about 600 pro- 
tons/em. 

In a series of papers (9) Parker has summarized various 
geomagnetic effects which are assumed to be associated with 
discrete solar streams from an active sun. He goes on to sug- 
gest possible mechanisms for suprathermal particle generation 
as a possible source of auroral effects and considers cosmic ray 
modulation by solar plasmas as well as magnetic storm effects. 
lrom these considerations, densities for quietest conditions 
appear to be at least 10? ions/cm’ with a minimum velocity of 
500 km/sec. Under active conditions, densities of 105 
ions/em? and velocities of 1500 km/sec or more are obtained. 

Russian workers (10) studying the dissipation of a high 
temperature nonstationary solar corona in the presence of a 
directed stream of matter also conclude that the density of 
particles in solar streams near the Earth is between 10* and 
10° particles/em’. Summaries of observations together with 
some views on the nature of the solar plasma prior to the 
launching of the Russian space probe of Jan. 2, 1959, have also 
been given (11). Velocities of between 1000 and 3000 km/sec 
are suggested with stream concentrations ranging from 10 to 
10° particles/em*. With the low density stream, a background 
gas at about 10‘ deg K is also considered to be present. 

It should be emphasized that the properties of the solar 
plasma under discussion pertain to conditions away from the 
perturbing influence of planetary magnetic fields. Such fields 
will shield out or markedly change the influence of the solar 
plasma in nearby regions of space. In the case of the Earth, 
we may compute the approximate distance beyond which solar 
plasma effects to be discussed will be active. Roughly, the re- 
quirement (9) is that the plasma kinetic energy density exceed 
the magnetic energy density, and we may write this condition 
as 


nmy? = — 
8 


In the equatorial plane of the Earth, assuming a pure dipole 
field, we have ; 


B = 
from which the minimum distance for onset of the effects is 
r = ro(Bo?/4anmv?)'/* 


Assuming the Earth’s surface field (horizontal component) 
as By = 0.4 gauss, setting 7p = 6.4 X 108 cm as the terrestrial 
radius and taking the stream density n to be 10° particles/em* 
with v = 10% km/sec and m appropriate for protons, we find 
R= 291. For this illustrative case, we would expect effects 
beyond 8000 miles from the Earth’s surface in the plane of the 
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geomagnetic Equator. This is to be taken only as an order of 
magnitude figure for a variety of reasons; particularly be- 
cause of the perturbation of the Earth’s dipole field by the 
plasma itself, because of the anisotropy of the pressures in- 
volved, and because of the probable presence of a radial solar 
magnetic field. It does appear that “‘stationary”’ satellites 
orbiting at 22,000 miles might be affected. Note that r varies 
slowly with n and ». 


Mechanical Forces Due to Solar Plasma Streams 


It is of interest to compare the force per unit area exerted 
on a material body by solar plasma streams with that produced 
by light pressure. The latter, in general, is given by the sum 
of momenta carried by incident solar radiation and the anti- 
parallel component of reflected or absorbed and subsequently 
reradiated flux minus the parallel component of momentum 
from reradiated flux. If we assume a perfect reflector, no re- 
radiation can occur, and the force per unit area at 1 A.U. 
reaches a maximum value of about 9 X 10~° dynes/cm?. 

If we assume a solar plasma containing, say, 600 particles per 
cm} and traveling at an average velocity of 1000 km/sec, we 
obtain a force of 10-5 dynes/cm? if the beam is simply stopped 
in the material. For an active sun during a severe storm, with 
stream densities of 10° particles/em* and velocities of, say, 
2000 km/sec, we obtain 6.7 X 107% dynes/em*. Thus, for 
quiet conditions the plasma produces 0.11 times the maximum 
light pressure, but during a storm it may exert in ‘“‘gusts’’ per- 
haps 75 times the solar radiation pressure on a perfect reflector 
at 1 A.U. from the sun. 

The very difficult problem of locating and tracking a small 
instrument package in the vicinity of the moon, for example, 
has led to the very interesting suggestion that highly reflective 
balloons be periodically discharged from the package (12). 
According to this suggestion, dispersion of these balloons by 
light pressure would then produce an advancing line of objects 
of high visibility. If such a system were indeed employed, it 
would be of interest to attempt to observe additional ac- 
celerations above the light pressure value which could be as- 
cribed to the solar plasma. Conceivably, some information 
on the fluctuations in momentum transfer due to changes in 
solar plasma velocities and densities could also be obtained 
from such observations. No telemetry of data would be re- 
quired. Consideration should, however, be given to possible 
destruction of the balloons or their reflective coatings by the 
plasma. This will be discussed more fully in the next section 
since these effects are also germane to questions of long-lived 
passive communications reflector satellites and the durability 
of balloon markers (13) for the lunar surface assuming the ab- 
sence of a lunar magnetic field (14)* and to power generation 
equipment using thin films. 

Before leaving the subject of the forces exerted by solar 
plasmas, mention may be made of an apparent revival of in- 
terest in the old idea of using light pressure for propulsion. 
Originally suggested in the literature of science fiction (15),‘ 
it has received attention in the serious press (16) and in this 
Journal (17). The practicality of these ideas has already been 
questioned on other grounds by Greenwood (18), and there is 
no need to repeat his criticisms here. To the list of problems 
may be added the thought that in the presence of a strong 
solar plasma stream, the forces on the shrouds of a “solar 
sail” will be orders of magnitude greater than the results 


3 Recent Russian observations on gaseous emission from the 
moon have been interpreted by D. H. Menzel as giving evidence 
that the solar plasma can indeed reach the lunar surface. There 
are a number of widely accepted reasons to expect a very small 
lunar magnetic field. 

4 The author is indebted to a number of his colleagues for gen- 
eral information on the history of this idea in science fiction, and 
to the Editor of Astounding Science Fiction for providing the 
most pertinent reference ( 15). In spite of the fictional nature of 


the magazine, this particular article is a detailed technical eval- 
uation of light pressure propulsion. 
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given in (17), for example. In the latter paper the author, 
who opines that light pressure propulsion is practical, pro- 
poses the use of commercially available metallized plastic film 
of 0.1 mil thickness and suggests, as an improved design, that 
film of thickness 2 X 10-* cm be used. Assuming such a 
structure could be built, if it were struck by a solar plasma 
stream, considerations to follow show that its function might 


be seriously impaired inashort time. 


Sputtering by SolarPlasmas 


When high velocity ions or atoms impinge on a solid surface, 
momentum transfer processes or in some cases chemical proc- 
esses lead to the ejection or sputtering of atoms from the sur- 
face. For the case of physical sputtering of metals, the thresh- 
old energy Vo for particles of atomic weight M impinging nor- 
mal to a surface of atomic weight M, is given by 


1.6 X 10°(M + M,)¢d 
Yo = 


2 
) electron volts 


where 


heat of sublimation of the surface material, on 


© 
_ This semi-empirical expression was obtained by Wehner 
(19)° and fits experimental observations obtained with thick 
samples at room temperature or above. Applying these re- 
sults to the sputtering of very thin structures or coatings by 
the hydrogen-rich solar plasma stream can only give approxi- 
mate results because of: 

1 The possibility of chemical sputtering through metal 
hydride formation which would require much lower threshold 
energies. 

2 Possible differences between thin film properties and 
bulk mechanical properties which can change the kinematics 
of momentum transfer—empirically, this might be accounted 
for by changes in ¢ and W. 

3 Differences between ionized and neutral atom sputtering 

mechanisms (20, 21). 
For protons on aluminum, the calculated threshold for per- 
pendicular incidence, using ¢ = 75 kcal/mol and W = 5.1 X 
10° cm/sec, is Vo = 560 ev. The energy of solar protons 
may be readily computed and, for 1000 km/sec streams 
is ~5100 ev. Thus, even for considerably slower streams 
than suggested by some of the observations, sputtering thresh- 
olds at normal incidence are likely to be exceeded by large 
margins. 

The yield of sputtered atoms per incident particle is a dif- 
ficult question, and observations as well as theoretical predic- 
tions vary widely. Depending on the circumstances, yields 
ranging from 0.05 to well over 1.0 might be appropriate for 
our situation on the basis of the measurements reported by 
Weiss and by Wehner and the calculations of Whipple (22). 
Still higher yields would presumably be observed at other than 
normal incidence; Goldman and Simon (23) suggest a sec 6 
dependence in the theoretically simple case of high energy 
sputtering with the yield going with energy as In V/V. It 
should be mentioned that the sputtering properties of some 
compounds, especially oxides, are not well understood. This 
gap in our knowledge is important both from the standpoint 


bulk sound velocity for the material, cm/sec 


5 The author is indebted to Dr. Wehner for comments on the 
present state of knowledge on certain aspects of the sputtering 
process. 

Note added in proof: Wehner has recently used spectroscopic 
methods to measure sputtering thresholds and no longer feels the 


in the present context will be discussed in a separate paper. 
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of understanding the relationship between sputtering and 
meteorite longevity, a subject to be discussed further else- 
where, and also in the present situation. Sputtering studies 
in the laboratory are frequently done by methods such that 
much, if not all, of the available experimental data reall, 
refer to oxide films. This is especially true in the case of 
aluminum. In space, atomically clean metallic surface films 


could result from a variety of processes involving both irra- 
diation and chemical reactions. 
_ involved, we shall assume a sputtering yield of unity and note 


In view of the uncertaintie: 


that the destructive effects to be computed can be readil; 
scaled to other values. 

Consider now a reflective structure consisting of metallize:! 
plastic. In order to achieve high visible opacity, a meta! 
thickness of the order of 300 A will be required which, in the 
case of aluminum, would weigh about 10-' gm/cm? and contai: 
about 2 10” atoms/cm?. Supporting the metal coatin: 

would be a plastic film which, for marker balloons or focusin:; 
— eollectors, might be 0.1 to 0.25 mil thick and which woul 


weigh between 25 X gm/cm? and 60 X gm/cm: 


- Thus, metal thicknesses in excess of 10 times the assume: 


A value (e.g., 3000 A of aluminum) begin to cost heavily i: 


payload for fixed structural area. For light pressure propul- 
sion where plastic film thicknesses of 2 X 10-5 cm have been 
proposed (17), even the 300 A coating adds one third to the 
mass per unit area. 

Complete stripping of a coating by sputtering would occur 
in a time ¢ given by 


Vi t = NLd/Mo<nv> 
N = Avogadro’s number 
L = coating thickness 


0 sputtering yield assumed independent of L but 
averaged over angle and weighted according to the 


incident velocity distribution 


M coating atomic weight ~ 
d = density 4 


<nv> = the solar flux averaged over the exposure time 


Optical properties will, of course, begin to change well before 
complete removal is effected. This expression does, however, 
ignore atoms removed from one part of the structure and de- 
posited on another region, either because of the physical de- 
sign of the structure or by virtue of image charge forces. 
Most sputtered particles are neutral. We also ignore atoms 
ejected in the forward direction through the structure. Mass 
losses or atomic displacements through an interface between 
materials from this latter mechanism could give quite interest- 
ing effects especially for very thin films. 

The time required for solar protons to destroy a 300 A coat- 
ing, on the basis of an average yield of unity, a density of 600 
particles/em* and a velocity of 1000 km/sec would be less than 
one month of exposure. For densities and velocities cor- 
responding to an intense storm (v = 1500 km/sec and n = 
10°/cm*) the coating would be completely removed in about 
104 sec. This time is within an order of magnitude of that re- 
quired for an average beam from a single solar event to pass 
the Earth at a velocity of 1000 km/sec. A beam length in 
space of about 10'' em is suggested by Kiepenheuer (24) to 
account for the time dependence of geomagnetic effects, whereas 
Unsold and Chapman (4) give a beam length of 1.2 X 10!* cm. 
Thus, depending on how literally one takes the assumptions, 
a single encounter with a beam ejected from an active sun 
could conceivably destroy the reflective coating. Scaling to 
other assumptions about proton densities, velocities and 
yields is direct and obvious. 

We turn now to an aspect of the solar plasma usually ig- 
nored—the presence of elements heavier than hydrogen. 
Except for Unsold and Chapman’s attention to possible Ca II 
content of a stream, a plasma consisting purely of protons and 
electrons is generally assumed. However, since the corpuscu- 
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lar beams appear to be merely segments of the solar corona 
blown out bodily, one should expect both helium and higher Z 
elements to be present. If one assumes that no fractionation 
occurs during the explosive acceleration, an assumption which 
might be more applicable to individual streams from an active 
sun than to quiet conditions, it would be reasonable to expect 
a plasma composition similar to that of the corona itself, by 
weight roughly 75 per cent hydrogen, 23 per cent helium and 2 
per cent heavy elements (25). For the last we may take an 
average atomic weight of around 32 (26). In the plasma a 
population ratio, H:He:Heavy, of about 75:6:0.06 might 
then be expected. Charged heavy species would travel at the 
same velocity as the general proton stream because of initially 
high collision frequencies and polarization fields and would 
therefore possess correspondingly higher energies and mo- 
menta. From the viewpoint of any impulse delivered to a 
structure, the influence of all but the helium is negligible (and, 
in fact, is partially taken into account by the methods used 
to calculate stream properties from observations in the first 
place). From the viewpoint of sputtering, however, the pres- 
ence of these components could lead to appreciable effects. 
In comparison to protons, momentum transfer from such 
particles to atoms of a surface is more efficient, differences be- 
tween charged and neutral components not so marked, anda 
variety of chemical sputtering mechanisms possible. For ex- 
ample, the work of Weiss et al. (21) for silver shows a sputter- 
ing yield for He~ higher than for H* by a factor of about 9 
in the 10-kev region. This alone would make the sputtering 
by solar He* about equal to that by solar protons. 

The heavy component in a 5-kev proton stream might have 
an average energy in the vicinity of 160 kev, whereas the 
sputtering threshold for such particles incident normally on an 
Al surface would drop to the order of 94 ev because of higher 
momentum transfer efficiency. It might be noted that if such 
particles were present even in the form of a cold stationary gas, 
perhaps the remnants of previous solar streams, and the 
structure under consideration were moving along with the 
Earth’s orbital velocity (but outside the magnetic shielding 
limit, e.g., near the moon), the equivalent energy of the gas 
impinging on the surface would be about 180 ev—well above 
the sputtering threshold. 

In passing it should be noted that considerations similar to 
the foregoing may be applicable to questions of longevity of 
thin coatings of prescribed properties applied, for instance, 
with the purpose of temperature regulation of extraterrestrial 
structures whether or not these are constructed of especially 
thin members. 


Radiation Damage by Solar Plasmas 


If the surface of the plastic film itself, rather than an over- 
lying metal coating, were exposed to the solar plasma, radia- 
tion damage effects over and above those caused by solar 
ultraviolet and x-ray emission are to be expected. Hydrogen 
evolution leading to charring and carbonization at exposures 
of more than 10” rad (1 rad = 100 erg/gm) should occur. 
If the reflecting film were behind the plastic, a marked drop in 
albedo would presumably precede final destruction by some 
sort of sputtering away of the carbon skeleton of the polymer. 

Radiation damage produced by relatively slow, heavy 
particles of the kind present in solar plasmas has not been ex- 
tensively studied. As a first approximation, we shall assume 
the dose to be delivered uniformly over the range of the inci- 
dent particles. This is justified because multiple scattering 
and similar effects will smear any strong dependence on details 
of the primary energy loss mechanisms. Since hydrogen 
evolution results from molecular excitation and ionization, 
direct conversion from energy loss to chemically effective 
dosage should be permissible down to a limit below which 
elastic collisions predominate as a result of the incident 
particle velocity dropping below the pertinent orbital electron 
velocities. Seitz (27 
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which the incident particle velocity is that of an electron with 
about one eighth to one quarter of the required excitation 
energy L,. If H,~ 4 ev, the proton energy required for excita- 
tion would be 1 to 2 kev. For the heavy element component, 
the lower limit would be around 30 to 60 kev. In both cases, 
the velocities associated with the solar plasma are well above 
the required limits. For a rough estimate we shall compute the 
range of the particles and then convert the energy loss sus- 
tained up to these limits to rads. 

For solar protons, the range measurements of Cook et al. 
(28) may be used along with a range-energy relation of the 
form R = kE*/* to obtain a range of ~10-* gm/cm? for a 
velocity of 1000 km/sec. For a density of 600 particles/cm? 
the dose would be about 5 X 10° rad/sec. Carbonization re- 
quires of the order of 10” rad and hence would occur in about 
2 X 10‘ sec. For a beam from an active sun, the time re- 
quired would be about 100 sec or less than the exposure time 
provided by a single major event. 

For the heavy elements discussed earlier, the range may be 
computed from the expression (29) 


(Z,2/8 + 7,2/3)1/2 (* + As gm 
= 6 1077 — 


where 


Z;, Ai, Z2, Ag = atomic number and weight of the absorber 
and stream particles, respectively 
= E = stream energy, kev 


The range of this component in plastics is thus calculated to 
be 5 to 10 X 10-*gm/cm? for velocities presumed appropriate 
to a quiet sun. For storm conditions, with velocities of 1500 
km/sec or more, the range of the particles would exceed 25 < 
10-5 gm/cm? allowing them to completely penetrate a 0.1- 
mil structure. The dose delivered to these thicknesses of ma- 
terial is of the order of 3500 rad/sec for the quiet sun, as- 
suming 8 X 10~‘ heavy particles per proton, and might rise 
to perhaps 5 X 10°rad/sec asa result of anintense beam. The 
corresponding time to deliver 10’ rad from the heavy com- 
ponent alone is about 40 days in the former case. The dosage 
from the heavy component is naturally small compared to the 
proton component but is effective over larger thicknesses. It 
does not require any sputtering away for total penetration of 
sufficiently thin sections. 


Conclusion 


The detrimental effect of solar plasmas on the performance 
of various thin lightweight structures which might be erected 
in space outside of magnetically protected regions has been 
shown to be potentially serious. The precise magnitude of 
the problem cannot be established with certainty, but assump- 
tions which appear reasonable in the light of present knowledge 
lead either to low durability or to high payload penalties. For 
certain effects a reduction of plasma density by one or more 
orders of magnitude, as a result of more complete experi- 
mental data collected in the future, would still imply a con- 
siderable problem. It will be of interest to review these sug- 
gestions when such information becomes available for both 
quiet and active conditionsonthesun. 
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Compressible Flow of a Dissociated 
a Diatomic Gas Past a Catalytic 
Flat Plate’ 


- PAUL M. CHUNG? and AEMER D. ANDERSON? 
Ames Research Center, NASA, Moffett Field, Calif. 


HE BOUNDARY layer concept was first applied to 
catalytic reaction problems by Chambré and Acrivos 
(1,2,3).4 At the stagnation region, affine closed form solu- 
tions were obtained for catalytic surface recombination by 
Goulard (4) and by Chung (5). Solutions for regions other 
than the stagnation point are not so simple. In (1), two 
solutions, one exact and one approximate, were obtained for 
the incompressible flow of a reactant past a catalytic flat 
plate, assuming a first-order surface reaction. Both solu- 
tions are in seriesform. The exact solution has a small radius 
of convergence, making its applicability limited. Rosner (6) 
presented a rather comprehensive treatment of the incom- 
pressible case, including a solution for a first-order reaction 
on a flat plate. 
This note concerns the compressible frozen flow of a disso- 
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ciated diatomic gas past a flat plate upon which a first-order 
surface recombination. occurs. The surface concentration 
distribution, which determines the amount of recombination 


_ of the atomic species, is to be determined. 


The equation of state for such a dissociated gas is 


= (1 + me.)p RT/M 
The boundary layer equations are 


oz 


pu 


or _ oT ou \? 
cpu dy («27) + [4] 


O 
0% [5] 
oy oY oy 


The usual boundary conditions apply to Equations [2, 3 and 
4]; however, the boundary condition at the wall for Equation 
[5] depends upon the chemical kinetics of the surface reaction 
and is 


dc PM 
=& = k{ — 6) 


Equations [1-5] can be solved by the von Kaérmén-Pohl- 
hausen integral method described in (7). The ratio pu/p*u* 
the Schmidt number, the pressure and the wall temperature 
are considered to be constant. This method is valid for Sc 
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of the order of magnitude of unity. The numerical calcula- 
tions reported here are for Sc = 0.72. 

With the previous assumptions, solution of the energy 
equation is not necessary to the solution of the concentration 
equation. 

The momentum equation was solved in (7) using a sixth 
degree polynomial for the velocity profile. The results, which 
are used here, are 


= u/Ue = 2n — + — 278 (7) 
Sr\* 
= (5) Cue Ag [8] 


where A = 36,036/985. When Equation [5] is integrated 
with respect to y, the resulting integro-differential equation 


is, after a change in variables 
dF Fdx 1 

— +-— == 9 


A fifth degree polynomial is used here to approximate the 
concentration profile. The boundary conditions on this pro- 


file are: 
Atn =0 
an c's: 1 + me. on? 
Atyn=1 
2 


The reference state ( )* is chosen at the wall temperature 
with zero concentration, so that po = p*/(1 + moc.). 
Upon evaluation of F, Equation [9] becomes 


Exact solution 


(Ref. 1) 
0.6 
x Eqs. (17) and (18) 


Ref. 6 

Stagnant film solution 
4 

(Ref. 1) 
0.2 

—— 
0 | 
z 


Fig. 1 Solutions for moc. = 0 


the boundary condition that my is finite at £ = 0, is 
1-—(1 + K,£'/2)— (1+495/164 Sc) 
(1 + 495/16ASc)K,£'/2 


This solution is valid for all values of K;. It shows that mp is 
1 at € = 0. This was not an a priori assumption as it was 


m = [15] 


+ Kig*6*) (dB/dg) + — (Kig — K;)B* — B = 0 in (1 and 6). 
[12] In terms of the parameter used in (1) 
| 2 
K 128 A’/? SckLp* 
i= 
495 the solution is 
= (2 + 495/16ASc)K, 
1 — (1 + — (1+0.846/8e) 
: B = 1/(1 + moce) "The approximate solution given in (1) is 
= 1 — 0.7312 + 0.4532? — 0.25223 ++... [18] 


It can be seen at this point that solutions to Equation [12] 
will be in terms of three dimensionless parameters, K,, Sc 
and K;, with K; being the free stream compressibility. The 
parameter K, varies from about 10~* for almost inert sur- 
faces at low pressures to about 10? for highly catalytic sur- 
faces at high pressures. 

When the compressibility at the wall, 1 + moc., is suffi- 
ciently near 1, Equation [12] simplifies to a first-order linear 
equation 


dmg + + 2) 
dg (1 + 2 (1 + 


[14] 


This corresponds to the problem treated in (1 and 6) where 
compressibility was not taken into account. 
The closed form solution to Equation [14] which satisfies 
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The present solution, Equation [17], is compared to Equation 
[18] and other known solutions in Fig. 1 for Sc = 0.72 and 
C = 1.0. The present solution gives results almost identical 
with the exact solution, which is for Sc = 0.72, in the interval 
of convergence of the exact solution. 

For finite moc., no closed form solution of Equation [12] 
can be obtained analytically. A solution based on a power 
series expansion about the point ¢ = 0 was found to con- 
verge rapidly for small values of K;. If all but the first two 
terms are neglected, the solution is 


= (2K;? — K.y)/[2K;3? + — 1)K2¢] [19] 


This solution agrees quite well with the exact solution for 
< 0.2. 

Numerical solutions to Equation [12] were obtained using 
an IBM 704 digital computer for the conditions not covered 
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concentration integral, Equation (10) 
specific catalytic reaction rate constant 
parameters defined by Equations [13] 
length of flat plate 
C/Ce 

molecular weight of undissociated diatomic gas 
0.23 + 0.77mo 

pressure 

universal gas constant 

= Reynolds number, p *u,L/u* 
Schmidt number, 
absolute temperature 

x component of velocity 

y component of velocity 
coordinate along plate 
coordinate normal to plate 
dimensionless variable, defined by Equation [16] 
boundary layer thickness 


=3 


e@— Typical numerical solutions 


it 


4 
° Eq.(20) 


ll 
° o 
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KS r = dimensionless variable, defined by Equation [8] 
= dynamic viscosity 
Fig. 2 Correlation of Equation [20] with numerical solutions g = 2/L é 
p = density 
1 
by Equations [15 and 19]. A study of these numerical solu- ” all ide 
tions revealed that for each K, the solutions were congruent Subscripts 
with Equation [15] when m was plotted vs. ¢/K3", where - 
n = 0.23 + 0.77m. The form of Equation [15] was then : ne - ot boundary layer 
utilized to bring the solutions together for all values of Ki. * me elena ahielie 
The semi-empirical closed form solution thus obtained is 
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fia Vapor-Liquid Equilibria of Ozone- oY, HE VAPOR-LIQUID equilibria of ozone-oxygen systems 


at 78 to 93 K and at comparatively low total pressures 


Oxygen Sy stems’ have been the subjects of several recent publications (1-4).’ 
addition, the vapor pressures of the individual components 
_ CHARLES K. HERSH,? ROBERT I. BRABETS,? in their pure state have been determined (5,6). The ob- 


GERALD M. PLATZ,4 RAYMOND J. SWEHLA® and jectives of this investigation have been to verify the 1- 


DONALD P. KIRSH* atm vapor-liquid equilibria curve and to extend the data to 
f 20 

Armour Research Foundation, total premureo 

Illinois Institute of Technology, Chicago, III. _ Experimental Procedure 

Received Aug. 27, 1959. The experimental procedure consisted of determining the 
1 The authors are grateful to C. Brown and M. J. Klein for vapor and the liquid curves separately and by different tech- 
niques. When a pressurized vapor of known concentration 
3 Pr Member ARS. is slowly cooled, the temperature at which liquefaction first 
‘Supervisor, Propellant Research. Present address: U.S. occurs can be determined from a study of the temperature- 
ARS. pressure graph. When the vapor volume in a closed system 
ociate Engineer, Propellant earch. resent address: i iqui thi ized 
Oo. Chichen, Th. Menaber ARS. is kept small, the amount of liquid which must be vaporize 
6 Assistant Propellant Numbers i in indicate at end of paper. 


ARS JourRNAL 


he) 
08 
= 
0.6 
~ 
| 
| 
~ 
tr 
fo 
in 
an 
ch 
of 
wa 
ex] 
cor 
Val 
Fo 
020 
J 
lary 
apy 
nee 
8 
Ma 


cor 


to provide 20 atm of pressure will also be small. However, 
when a relatively large liquid volume is used, the change in 
liquid composition can be considered negligible. In both 
of these techniques, it is essential that the composition be 
accurately known. 

Therefore, the test equipment was divided into two sec- 
tions. The first section, which was constructed entirely of 


glass, was used for the collection, distillation and accurate __ 
High concentration 


measurement of the ozone and oxygen. 
ozone was separated from the effluent of a Welsbach ozonizer 

in a trap maintained at 90 K and 100 mm Hg absolute pres- 
sure. After a sufficient quantity had been collected, the 
ozone was concentrated to 100 per cent by vacuum pumping. 
The desired amount of ozone was distilled into the precali- 
brated measuring tube, where the volume was accurately 
measured by means of a cathetometer; then the ozone was 
distilled into the test chamber which was held at 77 K. 
Similarly, the oxygen was measured and transferred into the 
test chamber. 

The test chamber, constructed of stainless steel and con- 
nected to the manifold by a ball and socket joint, consisted of 
a }-in. pipe tee, a short length of small ID tubing, a hypoder- 
mic tube leading to a strain-gage type of pressure pickup and 
a high pressure valve. All joints were sealed with silver 
solder, A copper-constantan thermocouple was _ soldered 
to the tee. A thin coating of Thermon,® a high thermal 
conductivity material, was applied to the unit to equalize 
the temperature. A 1-in. insulating blanket of glass wool 
surrounded the chamber to lower the rate of temperature 
change. Three chambers were constructed and used inter- 
changeably to check the validity of the procedure. Each 
chamber had a different volume: 8.1, 12.0 and 14.0 ce, 
respectively. 

During a run the chamber was supported at only one 
point so that it could be completely detached from the mani- 
fold by removing the double ball and socket bridge. Thus, 
agitation of the ozone-oxygen solution could be maintained 
by swinging the chamber in a small arc. A slow rise in 
temperature was obtained by removing the refrigerant Dewar 
and allowing the chamber to warm in a draft-free location. 
A slow and steady decrease in temperature was obtained by 
placing the chamber in the vapor portion of a partially filled 
Dewar and allowing the cold gas to flow over it. 

The following procedure was used to obtain the vapor data: 

1 A sufficient quantity of material of the desired concen- 
tration was completely vaporized in the chamber, and a uni- 
form gas concentration was obtained by maintaining the gas 
at 240 K for 10-15 min. Since most of the gas was in the 
pipe tee and the connecting line, the concentration of the 
gas was assumed to be the same as that of the original liquid 
in the chamber. 

2 The chamber was slowly cooled and the temperature 
and pressure were recorded at intervals. 

3 Temperature vs. pressure points were plotted, and a 
change in the slope of the resultant curve indicated the point 
of initial condensation at which the known gas composition 
was in equilibrium with a liquid of unknown concentration. 

Because of the similarity of the pressure-time lines in a 
vapor-liquid equilibrium diagram, a minimum number of 
experimental points at different temperatures, pressures and 
concentrations are required. 

The vapor pressures of various ozone-oxygen solutions at 
various temperatures were measured in a different manner. 

For liquid data, a large quantity, usually about 2.5 cc, of an 
ozone-oxygen solution was prepared. The amount was 


large enough so that the small quantity vaporized did not a 


appreciably change the liquid concentration. In all tests 
below 175 K, the vapor pressure data showed that essentially 
all of the vaporized gas was oxygen. The amount of oxygen 
necessary to produce 1-20 atm of pressure in each apparatus 


* Thermon Manufacturing Co., Houston, Texas. 
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Fig. 2 Temperature-pressure relationships for ozone-oxygen 
systems 
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Table 1 Latent heat of vaporization and integral heat of 
solution of the ozone-oxygen system 
Integral heat of solution 
Mole percent Latent heat, keal/mole _ kcal/mole 
ozone kcal/mole solution ozone 
100 3.630 (3.626 )* 
90 1.422 +2.01 +2.23 
80 1.499 78 +2.16 
1.586 +2.07 
1.642 +0.79 +1.98 
1.642 0.39 +1.95 
10 1.642 +0.19 +1.90 
0 1.642 (1.636)? 
Reference (7). 
> Reference (8). 


was experimentally measured. Thus, the liquid concen- 
tration at any point could be calculated accurately. Data 
were obtained by slowly warming the test chamber and re- 
cording the temperature and pressure. The solution was 
agitated continually during the run to insure a uniform 
liquid composition. The accuracy of the procedure used to 
measure the temperature and pressure was established by a 
duplicate run on one sample. The data were reproducible 
to within 1 per cent or } deg K. 


Results 


A total of 29 runs was considered sufficient to obtain the 
vapor-liquid equilibria relationship. These runs produced 
approximately 100 points, and each run was plotted indi- 
vidually. The data are presented graphically in Fig. 1 which 


is the typical plot of vapor-liquid equilibria for a two-com- 
ponent system. The consistency of the data is shown in Fig. 
2 by the straight line produced from the plot of logio P vs. 
1/T K for various ozone concentrations (mole per cent). 
The calculated slopes of the constant composition lines from 
Fig. 2 are shown in Table 1. These slopes are proportional 
to the heat of vaporization of the solutions at constant com- 
position. Thermodynamically, this latent heat is equal to 
the sum of the heats of vaporization of the pure components 
and the integral heat of solution. The integral heat compu- 
tations were based on the assumptions that the vapor is ideal 
and that the sensible heat of the liquid and vapor is negligi}le 
compared to the latent heat. 

The heat of solution at infinite dilution is estimated to be 
+1.90 keal per mole of ozone (heat evolved). In Table 1, 
the values for the latent heats of pure ozone and pure oxygen 
in column 2 can be compared with the literature values of 
3.626 kcal per mole for ozone (7) and 1.636 kcal per mole ‘or 
oxygen (8). In both cases the error is less than 0.5 per ceut. 
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_ Measurements of the Heat Transfer 
_ Coefficients for Hydrogen Flowing 
in a Heated Tube’ 


J. M. FOWLER? and C. F. WARNER: 


Jet Propulsion Center, Purdue University, Lafayette, Ind. 


HIS note presents a portion of the results obtained dur- 

ing a recent investigation of forced convective heat trans- 
fer and the related frictional pressure drop for gaseous 
hydrogen flowing through a smooth electrically heated tube. 
During the experiments the ratio of the average absolute wall 
temperature to average absolute gas bulk temperature 7/7, 
was varied from 1.17 to 2.47 for average gas bulk Reynolds 
numbers ranging from 3670 to 62,200. The gas pressure in 
the test section was varied from 40 to 100 psia. The ex- 
periments were conducted in an electrically heated test sec- 
tion of 0.306 in. ID having a length of 40 diam. The flow 
conditions at the entrance to the heated test section were 
characterized by a fully developed turbulent velocity profile 
and a uniform temperature profile. 

Thermocouples attached to the exterior surface of the test 
section permitted the determination of values of the local 
heat flux input to the gas. In Fig. 1 is presented the vari- 
ation in local heat flux along the axial length of the test 
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section for three different gas flow rates. As can be seen 
from Fig. 1, the boundary condition for the gas at the tube 
wall was essentially that of constant heat flux. The slight 
variations in heat flux at the ends of the test section are 
attributed to heat transfer by conduction from the test sec- 
tion to the flanges welded to the ends of the tube for electrical 
connections. Heat balances for all of the experiments re- 
ported in this note were within the limits of +5 per cent. 
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Fig. 1 Variation in local heat flux with dimensionless axial 
distance zx/d for hydrogen; uniform initial temperature, fully 
developed velocity profile at entrance to test section 
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Fig. 2 Correlation of average heat transfer coefficients with 
constant heat flux; properties evaluated at film temperature 7, 


The average heat transfer coefficients were found to corre- 
late well, as can be seen from Fig. 2, by an equation presented 
by Desmon and Sams (1),4 = 0.023 (N 
over the range of modified Reynolds numbers from 8000 to 
30,000. Good agreement was also found by comparing 
the data for hydrogen with those of a recent and similar 
experiment performed with helium gas by Taylor and Kirch- 
gessner (2). In Fig. 3 is presented the experimental data 
for hydrogen together with the line of the correlating equation 
recommended by Taylor and Kirchgessner. It should be 
pointed out that the correlation for L/d developed by Taylor 
and Kirchgessner was obtained for test section lengths having 
length to diam ratios of 60 and 92, whereas the test section 
employed in the experiments with hydrogen had a length to 
diam ratio of 40. 

The one-half Fanning friction factors obtained from the 
experiments with hydrogen are shown in Fig. 4. The line 
suggested by von Karmén is presented in Fig. 4 for com- 
parison. From Fig. 4 it can be seen that the data for iso- 
thermal flow are in good agreement with the von Kaérmén line; 
however, the data for diabatic heating are in poor agreement. 
It should be pointed out that in computing the friction factor, 
the calculated momentum pressure drop was subtracted from 
the measured pressure drop to obtain the frictional pressure 
drop. 


Nomenclature 
f = Fanning friction factor 
T.. = absolute average temperature of test section wall, R 
T, = absolute average bulk temperature of gas, R 
T; = absolute average film temperature (7, + 7;)/2, R 


Nre = Reynolds number VDp/u 

Nrem = modified Reynolds number V;Dp;/n; 
Np, = Prandtl number Cpy/k 

vi = Nusselt number hD/k 


Numbers i in pare an omens References at end of paper. 
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Fig. 3 Correlation of average heat transfer coefficients with 
constant heat flux; properties evaluated at film temperature 7’, 
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Fig. 4 Correlation of average half Fanning friction factor for 

hydrogen. The average gas density was employed in evaluating 

the Reynolds number and the viscosity was evaluated at the 
average bulk temperature 


Subscripts 


b bulk, indicates properties evaluated at the bulk tem- 
perature 


f = film, indicates properties evaluated at film temperature 
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A Method for Heat Rejection From 


ROGER C. WEATHERSTON? and WILLIAM E. SMITH? 


Cornell Aeronautical Laboratory, Buffalo, N. Y. 


UNDAMENTAL to the design of closed-cycle power- 
plants providing high specific impulse for space propul- 
sion motors is the means by which waste heat is dissipated. 
Heat can be rejected to a space environment only by radia- 
tion. The conventional radiator for high power levels is 
very massive; above the megawatt power level, in fact, its 
weight is about half of the weight of the powerplant. <A dif- 
ferent type of radiator, described here, should weigh con- 
siderabiy less than conventional radiators. In addition, its 
operation should be relatively unaffected by meteoritic im- 
pacts, a serious design consideration for the conventional 
radiator. This new radiator concept may permit significant 
reduction of the weights of other nuclear powerplant com- 
ponents as well, through modification of the thermodynamic 
cycle. 

A conventional radiator relies on circulation of the power- 
plant working fluid or a secondary heat exchange fluid through 
a system of supply and return headers which serve as mani- 
folds for many smaller tubes interposed among them (1,2).* 
The exterior surfaces of these smaller tubes constitute almost 
all of the heat exchange area of the radiator. The headers 
are made of material sufficiently thick to withstand all 
meteoroid impacts that will probably be encountered during 
the life of the device, from several months to years. The 
smaller tubes can have thinner walls than have the headers, 
and, in the event of penetration by a meteoroid, a damaged 
tube can be sealed off from the system to prevent loss of fluid. 
Despite the use of thin-walled tubing, the total weight of the 
smaller tubes is about one half of the entire radiator weight. 

The heat exchanger proposed herein would typically con- 
sist of a rotating drum or other contact surface, which is 
heated either directly or indirectly by the waste heat of the 
powerplant, and a long thin flexible belt (Fig. 1). The cool- 
est portion of the belt contacts the hot drum, conducts heat 
from its surface, and stores the heat by virtue of the thermal 
capacity of the belt. Upon leaving the drum, the heated belt 
element travels through space until the heat lost by radiation 
is sufficient to re-establish its original temperature. The 
cooled belt element then returns to the drum, and the cycle 
is repeated. It is possible to contemplate a number of varia- 
tions of the foregoing scheme. The purpose of this note 
is to bring out the advantages of capacity coolers for heat re- 
jection in space. Therefore, the arrangement sketched in 
Fig. 1 should be considered as illustrative of a more general 
class of devices. 

Referring to the configuration shown in Fig. 1, one may 
show that the rate per unit area at which the belt conducts 
heat from the rotating drum is much greater than the rate at 
which it radiates to space. Consequently, the contact area 
between belt and drum is much less than the radiation area. 
Thus, the ratio of these areas, as well as the belt speed, could 
easily be controlled in flight to accommodate a range of drum 
temperatures and radiative power levels. Further, the aspect 
ratio of the belt could easily be made large enough that both 
surfaces of the belt radiate to space, with very little mutual 
occlusion by the inner surfaces. 

The belt system when compared to the tube-and-header 
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Fig. 1 Possible configuration of a capacity heat exchange sys- 
tem for a space radiator 


system is characterized by its simplicity, a fact from which 
numerous advantages can be derived. The belt does not con- 
tain a cooling fluid; therefore, its penetration by meteoroids 
should not significantly impair its operation. Consequently, 
very thin and possibly woven material can be employed for 
the belt. Such a belt could easily withstand boost loads, as it 
could be either rolled completely round the drum and 
spliced in space or “‘accordioned” on opposite sides of the 
drum and simply unfolded in space. 

The chief advantage of the belt system is, of course, its 
low weight as compared with a system comprising tubes, 
supply and return headers, and contained fluid. As an ex- 
ample of the saving in radiator weight expected for the pro- 
posed system, an approximate, though conservative, design 
is offered to accomplish the heat rejection task set forth in 
(1). Neglecting the weight of the small portion of belt in 
contact with the drum at any time, one can express the rate 
of heat lost from the belt in terms of the change of its heat con- 
tent and the belt cycle time as 


q = WeT, — T:)V/L 


where 
Ww = weight of the belt 
c = belt heat capacity 7 
T, — T, = temperature change belt experiences 2 
V = speed with which the belt circulates 7” 
L = length of belt loop : 
The belt must radiate energy at the same rate q so that 
q = 2DLeoT~ [2] 
where 
2D = effective radiating perimeter of the belt cross section 
€ = belt surface emissivity 
o = Stefan-Boltzmann constant 
The mean belt temperature is defined by 


T ={— T*dl 3 
(Sf ™) 


T is the instantaneous temperature of an element of length, 


ARS JourNAL 


; 
~ 
if VS 
GUIDE 
ROLLER 
- 
BELT LOSES a Z 
HEAT BY T, \ 
FROM BOTH 
SIDES \ 
Z 
4 
J 
BELT SPEED 
HEAT EXCHANGE 
BELT 
2 
t 
~ 
il 
a 
Dp 
> 
f 1 
V 


| 


dl. By elimination of L between Equations [1 and 2] the 
belt weight is obtained as seek 
W = q?/2DVeoT*c(T, — T2) [4] 


Now it can be shown that 


Equation [4] now becomes 

6DVeoTFc(1 T2/T,)? 

Finally it can be shown that, irrespective of other considera- 


tions, the belt weight is minimized when 72/7, = 0.69. The 
resulting equation for the minimum belt weight is 


Wnin = [7] 


[6] 


W 


To illustrate the weight comparison between the belt heat 
exchanger and the conventional one, consider the heat rejec- 
tion requirements for the 20,000-kw powerplant discussed (1). 
This example is chosen because it represents the composite 
design of a well-engineered space power supply whose radia- 
tor weight has been minimized on the basis of a conventional 
radiator design. The amount of heat to be dissipated from the 
radiator in (1) is 64,600 Btu/sec. The heat rejection tempera- 
ture of the powerplant is 1800 R. If 7; is 20 R less than the 
heat rejection temperature, i.e., 1780 R, and if € is 0.90, the 
belt weight (lb) is calculated to be 


W = 1.91 X 10°/DVc (8] 


where D is in ft, V is in fps, and c is in Btu/lb R. At the 
present time it is difficult to specify the practical limits for the 
parameters in Equation [8], but let D, V and c be conserva- 
tively chosen as 38 ft, 50 fps, and 0.10 Btu/lb R (steel). The 
belt would then weigh only about 10,000 lb, and its loop 
would be 430 ft long. The assumption of negligible mutual 
occlusion effects is borne out in this numerical example: 
Since the loop would be constrained to a circular shape or 
would be made into a Mobius strip, it is found that the ef- 
fective radiating belt width D is within a few per cent of its 
actual width. 

Since the contemplated vehicle acceleration is on the order 
of 10~4 g, there will be practically no load on the bearings of 
the rotating drum, and since there is no mechanism for at- 
tenuating the belt speed, the torque exchange between belt 
and drum would be nil. Hoop stresses and fatigue would be- 
come appreciable problems in belt and drum only when the 
belt speed reaches several hundred feet per second. With 
these considerations in mind, it may be practicable to use a belt 
speed of 100 fps. Furthermore, the average heat capacity of 
the belt could be elevated to at least 0.20 Btu/lb R, for ex- 
ample, by making a graphite matrix (c about 0.40 Btu/lb R) 
adhere on or between thin steel belt members. These changes 
would decrease the belt weight to 2500 lb. 

It is contemplated that sodium vapor would be fed from 
the turbine, as in (1), to the rotating drum where the vapor 
would condense on the interior drum surface and be drained 
and pumped to the primary heat exchanger, again as in (1). 
It can be shown that a 3-ft diam, 38-ft long steel drum, 0.025- 
in. thick, is more than adequate to accomplish the transfer 
of heat to the belt. The drum with internal bracing, end 
plates and bearings is estimated to weigh 2000 lb. Baffles, in- 
jectors, wipers, drains and other plumbing within the drum 
are estimated to weigh 1000 lb. To withstand meteoroid im- 
pacts the drum would be shielded by a steel plate, 0.125-in. 
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thick [the header thickness in (1)], which would weigh about 
1500 lb. The weight of the guide rollers, which may be widely 
spaced disks, and auxiliary structure is taken as 1500 lb. 
This 6000-lb drum system is, of course, a preliminary design 
that appears to be capable of accomplishing the stipulated 
tasks. Thus, the combined system of belt and drum would 
weigh about 16,000 lb at 50 fps belt speed and 0.10 Btu/lb R 
heat capacity or 8500 lb at 100 fps and 0.20 Btu/lb R.° This 
system would replace the 42,750-lb radiator in (1). 

Corresponding savings in weight could be illustrated for 
lower powered space vehicles. For example, it is estifnated 
that the 4800-lb conventional radiator for the 1-megawatt 
space powerplant described in (2) could be replaced by a belt- 
drum system weighing 1400 Ib. 

The first radiator weight comparison cited previously is 
normalized upon the thermodynamic cycle given in (1). 
Such cycles, for which the conventional radiator dominates 
all other components in weight, are designed primarily to 
minimize the ratio of radiator area to net power output. 
Whereas the conventional radiator area is intimately related 
to its weight; this is not true for the belt-drum system. 
Therefore the cited cycle design criterion and concomitant 
temperatures increase the weight of the belt-drum system un- 
duly. Furthermore, if the proposed belt radiator design proves 
to be practicable, and if significant weight savings should ac- 
crue, then the criteria for cycle design would not depend on 
minimizing the radiator weight. As a consequence, more 
efficient thermodynamic cycles and additional savings in 
powerplant weight may be contemplated. 

It will be recalled that the weight of the belt-drum system 
in the foregoing example was based on 7’; being 1780 R, about 
the optimum rejection temperature for the radiator and 
powerplant in (1). It is found, however, that for a given net 
power from the turbine the belt weight in this example is 
minimized with respect to the heat rejection temperature 
when 7; is 1600 R despite a concomitant increase in belt 
area. Equal belt weights can be attained, however, even when 
T, is reduced further, because belt materials of higher heat 
capacity may be employed. Possibly, because of high 
efficiency cycles associated with low heat rejection tempera- 
tures, one can afford to operate with moderate evaporator 
temperatures so as to alleviate materials problems throughout 
the nuclear powerplant. 

The object of this paper has been to acquaint the reader 
with the concept of capacity type heat exchangers for heat 
rejection from space powerplants. No pretense is made here 
regarding the finality of design. The heat transfer to and 
from the drum surfaces certainly will require critical examina- 
tion. Nonetheless what seems irrefutable is that the weight 
of the belt can be made almost an order of magnitude less 
than the weight of the tubes in conventional heat exchangers. 
It appears that the equipment associated with the belt can 
be made sufficiently light, relative to conventional radiator 
accessories, to maintain this advantage. 
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would be completely offset by the attendant decrease in length of 
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WILLIAM T. SNYDER! 


T IS a well-known fact that for isentropic flow of a perfect 
gas in a converging-diverging nozzle, the velocity at the 
nozzle throat is the acoustic velocity. That is, the Mach 
number is unity at the throat. For nonisentropic flow in a 
nozzle, however, there is sometimes a misunderstanding 
concerning the location of the M =1 condition relative to the 
throat. Durham (1)? states that ‘the throat velocity will be 
sonic in the actual case as well as in the isentropic case.” It 
is the purpose of this note to show that for nonisentropic flow 
the throat velocity is not sonic, and to discuss a physical 
mechanism which can account for this condition. Although 
the flow in a real nozzle is not one-dimensional, a one-dimen- 
sional analysis will be made with a nozzle efficiency being 
utilized to account for deviations from one-dimensionality. 
Because of the analytical simplicity of the one-dimensional 
analysis, this approach is frequently used in the design of 
nozzles for propulsion engines (2). 


Theoretical Analysis 


In addition to the assumption of one-dimensionality, the 
following assumptions will be introduced in the analysis: 
Steady flow, adiabatic flow, thermally and calorically perfect 
gas, and initial velocity zero. 

It is convenient to represent the nozzle process in the h-s 
plane. Fig. 1 illustrates the end states for an isentropic and a 
nonisentropic expansion between the same pressure limits. 
The governing equations can be written as follows: 


Energy 
Continuity 


constant 


The procedure for determining the Mach number at the 
throat is to obtain an expression for the mass velocity G as a 
function of M. At the throat, the mass velocity will be a 
maximum and thus setting the first derivative of G with 
respect to M equal to zero will determine M at the throat. 
The expression for the mass velocity as a function of M can 
be obtained as 


T 
(2) Ma = (2) ar fener 
Po Po To 


Q 
il 


(2 


Po 


The quantities po, a are properties of the reservoir and thus 
are constants for the system. Defining a new quantity G’ as 


G 
= 
Podo (2) 
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reduces the problem to that of setting the first derivative of 
G’ with respect to M equal to zero. The ratios 7/7 and 
p/Po must now be determined as functions of M. 

The following relationships between stagnation and static 
properties may be obtained from the governing equations 
and the properties of a perfect gas 


T, p, (k—1)/k 

(2) - (4) 
k-1 


The nozzle efficiency n between the reservoir and some section 
in the nozzle is defined as the ratio of actual to isentropic 
kinetic energy. Thus the expression for » becomes 

(V?)/2 (k — 1)TM? 


— 7.) — T/T) 


Equations [6, 7 and 8] may be combined and solved for 
the pressure ratio with the result 
k/(k—1) 


Pe | 1 — 1)M? 

Po 1+ (1/2)(k — 1)M? 
Substituting Equations [7 and 9] into Equation [5] gives tlie 
desired expression for G’ 


ps = 


n 1+ (1/2)(k — 1)M? 


Setting the first derivative of G’ with respect to M equal to 
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zero yields the following value for M at the throat b 


4/(k — 1) 


It is easily verified that for 7 = 1, Equation [11] yields 
M, = 1. The equation is rather formidable and difficult to 
interpret quantitatively for 7 < 1. However, a qualitative 


interpretation may be obtained by converting the equation 
into an inequality. By observing that 


Equation [11] may be written 


4/(k — 1) 
M,?i< 

= AP -3] 
n\k-1 

or 
2 
Mi< [12] 


3k(1/n — 1) + 3(1 — 1/3m) 


The heat capacity ratio k is always greater than unity, and 
thus the limiting maximum value of the right-hand side of 


inequality, Equation [12], occurs fork = 1. Putting k = 
gives 
[13] 


The throat Mach number is therefore less than unity for a 
nonisentropic expansion. To show the dependence of M, on 
9, Equation [11] is plotted in Fig. 2 for various values of k. 


Physical Mechanism 


The primary source of irreversibility in nozzle flow is the 
boundary layer which unavoidably develops. The boundary 
layer effects are confined to a relatively thin region in the 
vicinity of the solid boundary, and it is convenient to describe 
the boundary layer thickness in terms of the displacement 
thickness. The physical significance of the displacement 
thickness is the distance the solid boundary must be dis- 
placed in the direction of the fluid in order to create a non- 
viscous flow with the same mass flow rate as the real flow. 

In terms of nozzle flow, if the actual nozzle were replaced 
by a hypothetical nozzle corresponding to the displacement 
boundary layer thickness, the flow in the hypothetical nozzle 
could be considered as isentropic. The sonic velocity would 
then occur at the throat of the hypothetical nozzle. Since 
the foregoing analysis shows that the velocity is subsonic at 
the throat of the physical nozzle, it may be concluded that the 
throat of the hypothetical nozzle occurs downstream of the 
throat of the physical nozzle. Because of the small, but 
finite, radial component of velocity in the boundary layer at 
the physical throat, there is probably a vena contracta effect 
which accounts for the occurrence of the hypothetical throat 
downstream from the physical throat. 


Nomenclature 
a acoustic velocity 
A cross-sectional area 


constant pressure heat capacity 
mass velocity, mass flow per unit area 
G/ odo, dimensionless mass velocity 
enthalpy 

heat capacity ratio 

total mass flow rate 

Mach number 

pressure 

gas constant 

absolute temperature 

velocity 

fluid density 

nozzle efficiency 


Q 


G’ 


Subscripts 


stagnation conditions _ 
throat conditions 
isentropic process conditions 


~ 
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A Magnetogasdynamic 


M. D. COWLEY! 


University of Cambridge, Cambridge, England 


A new analogy with pure gasdynamics is given for the 
flow of a perfectly conducting gas. 
magnetic fields are parallel. 


SMALL perturbation theory was outlined by Sears (1)? 
for the plane steady flow of a perfectly conducting gas 
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The velocity and 


with parallel undisturbed velocity and magnetic fields. The 
analysis was made originally by Resler, but it has not yet been 
published elsewhere.* A simpler but more general approach 
to this problem can be found by forming an analogy with pure 
gasdynamics. 

We assume an ideal steady flow (viscosity, heat conduction, 
radiation and resistivity neglected) of a gas, where in MKS units 


E+vxXB=0 


One form of axially symmetric flow, which is consistent with 


3 Since the writing of this note, the work has become available 
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= (B, 0, B,) 

J = (0, jo, 0) 


In a steady state, where 0B/Ot is zero everywhere, it follows 
from curl E = 0, that Z, = 0, and v is parallel or antiparallel 
to B. We can imagine a plane system, which approximates 
flow in an annulus whose width is much smaller than the 
radius. 
The equation for the stagnation enthalpy ho(=h +. 30?) 
can be written 
‘so that Ao is constant along streamlines. 
With v X B = 0, conservation of the magnetic and material 
flux along a streamtube requires that B/pv is constant along 
_ the tube. Since there is no dissipation or heat flow, entropy s 
is another conserved quantity along streamtubes. We shall 
now assume that throughout the flow 


= kpv (k = constant) [8a] 

ne 8 = constant 

= constant [3c] 


f ‘Sin equations of state, and assuming thermodynamic equi- 


librium, all variables are now functions of each other. The 
theory will be developed without reference to the actual form 
of the equations of state. With the conditions [3], the only 
equations which we require are 


div pp = 0 [4] 
1 IBXB 
2 p p up 


ve the conditions [8b and 3c], Equation [5] degenerates to 


curl B X 


eurlv Xv = 
up 


Since curl v is perpendicular to v for the flows considered 
here, we obtain with the aid of condition [3a] 


B 
curl = + — curl B =kj 6 


where the positive and negative signs refer to B parallel and 
antiparallel to v, respectively. Thus the vorticity is propor- 
tional to the current density in the gas. From Equations [4 
and 6] we may derive equations in terms of the velocity com- 
ponents, as Kogan (2) has done for a perfect gas. It is easily 
shown that the conditions [3] are sufficient to insure that all 
properties are uniform in a region of rectilinear flow; both 
Kogan (2) and Resler started with the assumption of uni- 
form properties in such a region. 

From Equation [6] a vector may be constructed which re- 
mains irrotational. With the condition [3a] ite 


curl (v(1 — B?/pupv*)) = 0 


We define our irrotational vector by v* = v(1 — 62/v%), 
where 6 is the Alfvén wave speed (B/ ~/ up). It should be 
noted that v* is parallel to v, and the boundary condition for 
v* at solid walls is the same as that for v. The form of v* 
suggests that we should define p/(1 — 6?/v?) = p*. Then 

div p*v* = 0 [8] 
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a 


Bt B 
(pu-V)v + v(> + )- (B-V) = (pu-V) ne 


Defining p + B?/2u = p*, which is the pressure felt by a body 
at its surface, if there is no magnetic field w ithin the body, we 
obtain 
(p*v*-V)u* + Vp* =O (9] 

It is evident that the system of equations is now analogous to 
that governing the dynamics of an appropriate arbitrary non- 
conducting gas; p*, by the conditions [3] and suitable equ:- 
tions of state, has a functional relationship with p*. 

The speed of propagation of disturbances a*, in the new 
system, is given by 

* 2 b2 

att = = - =) [10] 


v2 v2 
where we have used the conditions [3], and a is the speed of 
sound in the absence of magnetic fields, V/ (p/p), The 
analogous Mach number squared is 
(v/a)*(v/b)* 
(v/a)? + (v/b)? — 1 


Values of M* are plotted against v/a and v/b in Fig. 1. The 
lines M*? = 0, M*? = 1, and M*? = + = divide the graph 
into the following five regions: 

1 When v > b and v > a, M* > 1 and the equations 
governing a motion are hyperbolic. Mach waves, at an angie 
sin! 1/M* to the streamlines, correspond to stationary in- 
finitesimal waves of the fast magnetogasdynamic type. As 
Kogan and Resler have shown, the waves are inclined down- 
stream from a disturbance to the flow. The gasdynamic 
analog is supersonic flow. 

2 When v > band» < a, M** < 1, and the equations are 
elliptic. The analog is subsonic flow. 

3 When v < bandv > a, M* < 1 and the equations are 
again elliptic. The definition of v* is such that v* is anti- 
parallel to v; p* is a negative quantity. v*, and p* may be 
redefined so as to be parallel to v and positive, respectively, and 
the analog is subsonic flow. 


20 
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Fig. 1 M*? as a function of v/a and v/b 
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4 Whenv <b, v <a, (v/a)? + (v/b)?> 1, M*?>1 andthe 
equations are hyperbolic. Mach waves correspond to station- 
ary waves of the slow type; Kogan and Resler have shown 
that the waves are inclined upstream from a disturbance. 
In some cases, such as linearized flow over bodies, it is con- 
venient to retain the definition of v* so that, for an ideal flow 
without a wake, the analog is supersonic flow in the reverse 
direction. 

5 When (v/a)? + (v/b)? < 1, M** is a negative quantity. 
The equations are elliptic, but there is no analog in gasdy- 
namics. 

In the hyperbolic regions, the Prandtl-Meyer function for 


plane flow is 
1 
o= ffcot sin ue) 
In terms of v, a and b, Equation [12] gives 


+ (v/b)? — 1)((v/a)? — 1)((v/b)? — 1) & 
(v/b)? — 1 v 
[13] 


[12] 


where the positive square root is taken. Kogan has given 
characteristic equations for the hodograph plane; if a Prandtl- 
Meyer function is derived from his results, it has the form of 
Equation [13]. If 7* = 1/p*, the steepening, or spreading, of 
waves depends on the sign of 7**d2p*/dp*? which may be 
written 


The quantity 75(02p/Or?), is positive for all known fluids. 
From Equation [14] we can deduce that compression waves of 
region 4 steepen when they are inclined upstream. 

When v = 8, there is a special solution, which is not ap- 
parent from the analogy. Chandrasekhar (3) has shown that 
then v, 6 and a etc. may be constant throughout the flow; 
the direction of flow and vorticity are arbitrary, and accelera- 
tion of the gas is given by magnetic forces alone. The solution 
only applies for special boundary conditions. 

The analogy is most easily applied to linearized theories of 
gasdynamics, where a relation between the pressure and 
density is not required (i.e., provided p* varies locally with 
p* in a reasonable manner, the form of the p*, p* relation does 
not matter). Aerodynamic forces for bodies, with no internal 


ll 


magnetic field, can be calculated from the distribution of p*, 
and p* can be found from the pressure coefficient of the 
analogous flow = 
C, = [15] 
In region 4, where the analogous flow is taken in the reverse 
direction, p* is negative; the p* variation is of opposite sign 
to the pressure variation of the analog, and there will be 
positive wave drag on a body in the magnetogasdynamic flow. 
Linearized theory of subsonic flow can be extended to cover 
region 5. 
The author wishes to thank Dr. J. A. Shercliff for his help 
and advice in this work. 


Nomenclature 


= local speed of sound a 

= Alfvén wave speed a 
= magnetic field vector (magnitude B) — 
= pressure coefficient Sow 
= electric field vector 
= enthalpy 
stagnation enthalpy 
= current density vector 

= a constant of a flow 

= Mach number 

= pressure 

= entropy 

= gas velocity vector (magnitude v) 
= magnetic permeability 

= density 

specific volume 

= Prandtl-Meyer function 


= 


Subscripts 


-- = quantity defined to have analogy in gasdynamics 
= free stream conditions ‘ 
r, 0,2 = directions in cylindrical polar coordinates — 
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Analog Computer Technique for 
Plotting Frequency Response 


JOHN L. WEBSTER' and DAVID N. SCHULTZ? - 
Chrysler Corp., Missile Division, Huntsville, Ala. 


This article describes a direct method of obtaining fre- 
quency response curves during analog computer studies. 
The system under study is driven by a constant amplitude, 
varying frequency sine wave. System output is then recti- 
fied and recorded on the Y axis of a recorder whose X axis is 
swept at the same rate as that of the driving frequency. 
A suitable driving function can be generated using analog 
computer components. 
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frequency range. 


Driving Function 


HE DRIVING function must be sinusoidal and of con- 

stant amplitude. The frequency must vary such that 
equal samples of response data are obtained in all divisions of 
To comply with these require- 
ments, an exponentially decaying sweep is_ utilized 
sweeping through the higher frequencies rapidly and dwelling 
sufficiently long at the lower frequencies to obtain an equal 
number of samples. 

The driving function is therefore 


E = asin [1] 


with the 0.01 resulting from the particular multipliers used. 
The amplitude is fixed by a and the frequency varies as 


ra 0.01 
== 


= 0.01bBe-«t 
sec 


dy 

[9] 
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For the X axis to vary at the same rate, the z displacement Therefore from Equations [4 and 6], the length of the re- 
becomes corded display will be 


= GBe~* in. [3] L=QB-—d) or L = X(1 —4/B)in. (7) 


The maximum and minimum frequencies of the sweep range 
ly at are fixed by the products 6B and bd respectively. The po- 
areas, x= GB in. las [4] tentiometer c may be used to control computing time and the 


from which 


Fig. 2 illustrates the computer setup used to generate these number of samples recorded per inch of X distance. 2 
functions. A biased diode is utilized as a clamper to stop the From Equations [3 and 4] 
exponential decay when 
d 0.01bd | 
W =00lbd—— or f= eps [5] 
ane 2x Therefore, the samples recorded per inch will be 
This will also stop the X axis sweep of the recorder when ___f_ _ 0.010B samples 


9] 


S= = 
= Gd |dx/dt| in. 
From Equations [2 and 5], computing time required to 
complete the sweep becomes 


DRIVING T = — log, — 10 
| GENERATOR Test Conclusion 
as With the computer setup as indicated by Figs. 1 and 2 and 
| sue adjusted as in the preceding equations, it is possible to plot di- 
yc rectly frequency response of any system being studied on the 


The system under test should generally be operating in the 
region of 0.1 to 20.0 rad per sec on the computer, and therefore 
time scaling may be required in many instances. 

To minimize integration errors, the samples per inch and 


, eet frequency range swept should be selected to keep the com- 
a ny puting time less than 2 min. 
1h To avoid transient disturbances associated with the sudden 
RECORDER < DETECTOR 
Fig. 1 Block diagram of frequency response manaeniing 

technique 
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application of the driving function, the problem should be 
computing before the sweep is started. 

Fig. 3 illustrates the type of frequency response curves ob- 
tained by this technique. 


Nomenclature 
E 


= output voltage 
initial condition voltage of integrator amplifier no. 2— 
gain of integrators 1 and 2 

time, sec 


B- = initial condition voltage of integrator amplifier no. 4 
c = setting of potentiometer P; 
d/2 = voltage set at arm of 
W = angular velocity, rad/sec : 
= frequency, cps 
X = total displacement from zero along x axis 
" G = X axis gain of recorder 
7 x = instantaneous displacement from zero along X axis 
LL = length of display, in. 
Ss = samples recorded per inch of display 
 T — = total time required to complete the sweep 


Some Characteristics of the Planar 
Satellite Orbit 


JOSEF S. PISTINER! 
_ The Martin Co., Denver, Colo. 


The classical use of the various anomalies to describe the 
orbital mechanics of satellites often becomes exceedingly 
complex. This paper, on the other hand, describes in 
relatively simple mathematical or graphical form some of 
the important orbital characteristics of satellites in terms 
of the true anomaly alone. In addition, the concept of a 
simple orbital energy diagram is developed and its utility 
indicated. The approach presented in this paper restricts 
itself to the planar two-body problem. « 


Orbit Ellipse 
IG. 1 defines the orbit ellipse with its significant pa- 
rameters and the center of a spherical celestial body at one 
of its foci. The orbit is assumed to be planar and in the gravi- 
tational (centra] force) field of the celestial body. A conserva- 
tive system has been assumed. Therefore, the satellite’s 
total angular momentum and total orbital energy (kinetic plus 


potential) are both constant. Or 
f = mrve = total angular momentum = constant [1] 
and from (1)? 
m um 
=; (vg? + 0, = total orbital energy = constant [2] 
Equation [2] may be rewritten as 
2E m\? m 
— (“") + [3] 
m 


Using the velocity parameter concept established by Alt- 
man (2), let 


and 
R? = 2E/m + C? [5] 


By use of the vis-viva integral it can also be shown that 


= —(um/2a) 
Equation [5], therefore, may be rewritten as 
2E [7] 
™m a 


Note that the parameters RF and C are constant for a given 
orbit and have the dimensions of linear velocity. 


~ Received Sept. 11, 1959. 
m, Design Engineer, Advanced Systems Section Staff, Denver 
iv. 


2? Numbers in parentheses indicate References at end of anal ‘ 


(6) 


‘it can be easily shown that 
= R/C 
It can further be shown that 


= C+ Reos¢ 


and 
= Rsin¢g 
Therefore 
veo C+Resg 1+ ¢ecos¢ 


Orbital Energy Diagram 


Equation [7] may be considered to define an orbital energy 
diagram if R? is plotted vs. C? (Fig. 2). Note that the dimen- 
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Fig. 2 Orbital energy diagram 
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of 45 deg represents all orbits of a given specific energy (or 
given semimajor axis). Also, since e = R/C, a straight line 
_ (through the origin) with a slope numerically equal to the 
square of a given eccentricity will represent all orbits with this 
m a eccentricity. Hence, the point resulting from the intersection 
of a line of constant specific energy with a line of constant 
eccentricity will uniquely define a given elliptic orbit. 

Note that the C?-axis intercept of a line of constant specific 
energy is numerically equal to the value of this specific 
energy. Therefore, the theoretical minimum change in 
specific energy required to transfer, for example, from orbit | 
to orbit 2 (see Fig. 2) is given by the difference i in C?-axi 


intercepts. 
Pp 


sions of R? and C? are energy per unit mass (or specific energy). 


‘Since 


True Anomaly as Function of Time 


It can be shown that - 


= (C/u\(C + R cos ¢)* 


dg f 


The integration of Equation [15] may be appreciabl) 
simplified by letting 


Fig. 3 True anomaly ¢ vs. time ¢ 
(t — to) = typ = time from pericenter passage 
then 


Go = 0 


_ With these simplifications, integration and simplification re- 


sults in 
i Rein ¢ 
i — R*(C + R cos ¢) 
2c _, | 


‘ies the period of an elliptic orbit is given by or 


= 7] 


= (Cc? — 


Equation [16] finally may be rewritten in terms of 7, e and ¢ 


e 


1+ ecos¢ 


s Fig. 3 is a plot of ¢ vs. t,,/r for several values of e. 
By use of Equations [9, 14 and 17] 


(1 + e cos ¢)? 


(nondimens ional) 


It can also be shown that 


. ay sin 
= — = (| 
dt? (; + e cos -) 
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Fig. 5 Normal acceleration r¢ vs. true anomaly ¢ v6 


Path Angle as Function of True Anomaly and Orbit 
Eccentricity 

The path angle is defined as the angle between the total 
velocity vector of the satellite and the local horizon. The 
tangent of this angle is given by Equation [12]. 

It is of interest to note that the maxima and minima of 
v,/ve (= tan 6) occur when > 


cosg+e=0. [22] 
When Equation [22] is satisfied, then 


Ur 
= 2 
V 1 — e? 3] 


Since v,/ve reaches @ maximum when 


4 
tan = 
van l-e - 


= — 
tan g = 


't can be shown that the time at which this maximum occurs 
is given by 


e 
9 [24] 
T 


Note also that for (v,/v9) max, in an elliptic orbit 
typ 
T 4r 


always. 


Radial and Normal Acceleration as Functions of True 
Anomaly and Eccentricity 

A knowledge of the radial and normal accelerations of a 
satellite in an elliptic orbit is useful in establishing the 
vehicle’s stabilization and attitude control requirements. 
Note that the magnitude of these accelerations is of the order 
of a small fraction of one gravity (q). : 


Radial acceleration 


Equation [2] may be written as 


1 ‘ 
E= 5 mi + r2g?) — constant [25] 
r 


Differentiation of Equation [25], assuming the vehicle 
mass m to be constant, and simplification results in 


Ferg [26] 


r2 


After further substitution and simplification, Equation [26] 
may be written in the following nondimensional form 


= ecos + e cos ¢)? 


Fig. 4 is a plot of ur /C* vs. ¢ for several values of eccen- 
tricity e. Note that for 47 /C‘to be an extremum 


sin g(1 + ecos ¢)(1 + 3e cos = 


Therefore, the following conditions for extrema — 


(a) sin g = 0 or g=ntT (n = 0, 1, 2) " 


Marcu 1960 


This condition is nonexistent since e < 1 for elliptic orbits. 
(ec) cos = —1/3e 


This condition exists only for e > 4.and results in an absolute 
minimum of pi 
Therefore 


Normal 
‘It can be shown that 


rg/r = —2tan [28 


Substitution of Equation [27] into Equation [28] results in 
the following nondimensional equation 


urg/C* = —2e sin + e cos ¢)? [29] 


Fig. 5 is a plot of ur¢/C4 vs. ¢ for several values of e. 
Note that for urg/C*4 to be an extremum 


cos ¢ + 3e cos? g — 2e = 0 


or 


Obviously, extrema exist only if a solution of Equations [30] 
results in 


-l<cossg<+l 


Nomenclature 

a = length of semimajor axis of orbit ellipse — 

= orbit eccentricity 

{ = satellite total angular momentum > 

m = satellite mass 

r = radial distance from focus of orbit sais to satellite center 
of mass 

t = time 

typ = time from pericenter passage 

t, = reference value of time 

v = total velocity of satellite center of mass 

ve = normal velocity component of satellite center of mass © 

v; = radial velocity component of satellite center of mass 

C = velocity parameter 

E = total energy (kinetic plus potential) of satellite center of 
mass 

R = velocity parameter 

6 = path angle (= angle between total velocity vector of the 


satellite center of mass and the local horizon) 

a = gravitational constant (= 1.41008 x 10% ft?/sec? for 
terrestrial orbits) 

7 = period of orbit 

true anomaly 

¢o = reference value of true anomaly 


ll 


a6 
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Symmetry of the Earth’s Figure’ inverse square. This force is small enough itself to be con- 
sidered as a perturbation: The equatorial and polar diam- 


eters of Earth differ by less than one half of one per cent. 

Although the effects of the equatorial bulge on the Vanguard 
I satellite are relatively small, some are cumulative. During 
the several thousand revolutions that the satellite has been 
observed, these effects have been verified observationally. 
The magnitude of these effects can be used to confirm (and 
even to improve) the estimate of the magnitude of the bulge 
which is obtained from geodesy (2).° 

By means of a special numerical integration carried out During the first revolution, the drag-free orbit of a satellite 
for nine hundred revolutions, it is shown that the oblate- lies roughly ina plane. The perturbation caused by the equa- 
ness of Earth causes long period changes in the shape, torial bulge can be described as first, a slow regression of this 
size and inclination of the osculating orbit of the Vanguard plane (without changing its inclination to the Equator), and 
I satellite (1958 Beta two). Unless many determinations second, a precession in the orbital plane of the perigee (closes: 
of these parameters are made per revolution, the effect of point). 
an Earth with equatorial symmetry can imitate in period, If a more accurate description of the orbit is attempted, i: 
amplitude and pnene he effect which might be expected is found that the rough term “orbital plane” is not precis: 
from a “pear-shaped potential. This result, which has enough for practical use. Fig. 1 shows that the osculating in- 
been verified by independent analysis using an appropriate clination varies considerably even during a single revolution 
Fig. 2 shows that there is also a relatively large variation o 
not support the hypothesis that Earth’s mass distribution osculating eccentricity during a single revolution. By « 


is asymmetric with respect to its equatorial plane. special analysis (3), it has been possible to derive expression~ 


: ’ for variation of perigee distance, of maximum latitude and oi 
F THE only force acting on an Earth satellite were the shape. 


J. L. BRENNER? 


Stanford Research Institute, Menlo Park. Calif. 


R. FULTON? and N. SHERMAN‘ 


University of California, Livermore, Calif. 


attraction of Earth itself, and if the mass of Earth were The results of that analysis have been confirmed by direc‘ 
spherically homogeneous, the orbit of a satellite would be an numerical integration of the differential equations of motior 
ellipse with the center of Earth at one focus. This is a simple (1). These equations are (in Cartesian coordinates) 


consequence of the fact that such a mass distribution produces 
a central, inverse-square force on the satellite. The drag 
of the atmosphere on the Vanguard I satellite is very small, 
and the perturbations on the orbit due to the sun, moon 
and tides of Earth are also small. The main perturbing force, 
due to the equatorial bulge of Earth, is neither central nor 


Received Sept. 21, 1959. 

1 This research was supported in part by the U. S. Air Force 
under Contract no. AF 29(600)-1773, and the U. 8S. Atomic 
Energy Commission, monitored by Air Force Missile Development 
Center, Air Research and Development Command. 


2 Research Mathematician. ember ARS. 5N ~— indics at e f 
Numbers in parentheses indicate References at end of paper. 


ECCENTRICITY, 0.18860 TO 0.19030 


TIME, O TO 300 MIN. 
Fig. 2 Variation of eccentricity during several revolutions 
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TIME,O TO 85 DAYS 


Fig. 3 Variation of eccentricity at perigee (above) and at apogee 
(below) during an 80-day period 


where 


L=T-V 
T = (1/2)(4? + y? + 2?) 


2 2 
r 

5 = (x? + y? + 22)'/2 


assuming that only the “second harmonic” of the gravita- 
tional potential field has nonzero coefficient. Here, R is the 
equatorial radius of Earth, and g(1 + J) is the algebraic 
sum of the observed acceleration of gravity at the Equator 
and the “centrifugal’’ force there. 

Fig. 3 shows the variation of osculating eccentricity at 
apogee and at perigee during some nine hundred revolutions. 
Both quantities have a 40-day period, but have opposite 
phase. From this figure it is seen that the term “eccentricity” 
not further specified can be expected to vary in any manner 
provided it stays within the extreme limits shown in the 
figure. To emphasize this possibility, we computed the 
osculating eccentricity at successive northward crossings of 
declination +30 deg. The result, shown in Fig. 4, shows the 
same period, amplitude and phase as has been used (4) to 
“prove” that Earth’s gravitational field has a nonzero third 
harmonic. This conclusion is therefore open to question and, 
in our minds, is not established. 

Our results and methods have successfully passed several 
decisive checks. 


= 


: 


ECCENTRICITY, 0.198675 TO O.189625 


TIME,O TO 90 DAYS 


Fig. 4 Eccentricity at ascending crossing of 6 = +30 deg vs. 
time in days 


We have been able to obtain spurious variations not only in 
this osculating eccentricity, but also in the osculating in- 
clination, rate of change of longitude of node, and rate of 
change of argument of perigee. Hence none of these quanti- 
ties, nor all of them together, give evidence concerning the 
existence of the third harmonic in Earth’s gravitational po- 
tential field. The magnitude of the coefficient of this har- 
monic (if its existence can be proved) is in our minds as likely 
to be 10-8 as 2-10~*, the value derived in (4). 
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Thermal Decomposition of Liquid Phase 
n-Propyl Nitrate and Its Compatibility 
With Various Metals and Other Additives 


ADOLPH B. AMSTER! and JOSEPH B. LEVY? 
U. S. Naval Ordnance Laboratory, Silver Spring, Md. 


Pure normal propyl nitrate has been found to undergo 
negligible decomposition at normal temperatures and 
only slight decomposition at its normal boiling point. 
The rate is accelerated by impurities, such as nitrogen 
dioxide, acetaldehyde, water and acetic acid. The experi- 
mental results also predict safe storage for at least a year 
at temperatures up to 160 F in containers of hard copper. 
ST aluminum alloy, 2S soft aluminum, K Monel, 303 stain- 
less steel and 1020 cold rolled steel. 


N EVALUATING the hazards associated with the handling 
of n-propyl] nitrate (NPN) it is appropriate to investigate 

its stability. The thermal] decomposition of gaseous NPN has 
already been examined (1),* yet it is well known (2) that in 
many cases the reaction mechanism changes in going from 
the gas phase to the liquid phase. The possibility that this 
might occur with NPN, coupled with the possibility that 
small amounts of impurities might affect won ees, led 
to the studies described herein. 


Experimental Procedures 


The NPN used in these experiments was eine ™ frac- 
tiona] distillation of material supplied by the Ethyl Corp. 
The commercial product was light amber in color, but the 
distillate was colorless. Two separate distillations were con- 
ducted in a column of about 20 theoretical plates. Material 
was collected at 53.6-53.9 C at 100 mm Hg pressure and at 
65.8 C at 168 mm Hg pressure. In each case, several frac- 
tions were taken and their infrared spectra compared with 
that of the commercial material. No differences were ob- 
served. All of the spectroscopic work reported herein was per- 


Table 1 Analysis of impurities in NPN using infrared 
spectroscopy 
detectable 
useful for concentration, 
Compound per cent 
n-propy! nitrite 12.70 0.3 
8.7 
9.45 
nitroethane 6.50 1.0 
ethanol 9.50 1.0 
n-propanol 9.9 1.0 
acetaldehyde 5.83 0.6 
propionaldehyde 5.85 0.2 
acetic acid 5.87 0.4 
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Table 2 Effect of heat on n-propyl nitrate 
Run Temp. Time Decomposition 
1 50C up to one week no perceptible 
amount 
2 85C up to 17 days no perceptible 
amount 
3 B.Pt(111C) 7 days no perceptible 
amount 
4 B.Pt.(110C) 12 days slight decompo- 
sition 


formed with a Perkin-Elmer model 21 double beam spectr« 


photometer using NaC] optics and cell faces. 

Small quantities of potential reaction or decompositio: 
products were added to pure NPN to determine the extent t») 
which they were detectable. The results, including the sub- 
stances added, the frequencies useful for analysis and the 
minimum detectable concentration are given in Table 1. 

Simultaneously, samples of NPN were heated in con- 
tainers open to the dry atmosphere, and spectra of the re- 
sulting liquids were compared with that of the original sample. 

The effect of increasingly vigorous heat treatment is indi- 
cated in Table 2 where it is seen that only in the most vigor- 
ously treated sample was there any detectable difference. 
Examination of the resultant spectra indicated that the po- 
tential impurities might include acetaldehyde, propionalde- 
hyde, acetic acid or a combination of these. The mechanism 
proposed for the gas phase reaction predicts acetaldehyde and 
n-propy! nitrite (1). Decomposition of the latter would yield 
propionaldehyde whereas oxidation of acetaldehyde by nitro- 
gen dioxide or air could give acetic acid. Quantitative analysis 
indicated that an insignificant quantity of the NPN had de- 
composed. Thus, even the relatively vigorous treatment of 
refluxing at atmospheric pressure for prolonged periods causes 


no major decomposition. 


The compatibility of NPN with each of a number of dif- 
ferent metals and other additives has been examined. In each 
case, 1 ml of commercial grade NPN was placed in a Pyrex 
glass vial with the addendum. The metal samples used were 
in the form of smooth round pellets } in. in diameter and ¢ in. 
in thickness, and were weighed before and after treatment. 
The sample vials were of proper size so that the pellets were 
completely covered by liquid. The contents were frozen and 
the vials evacuated and sealed. Each sample was then 
placed for a specified period in a bath maintained at 100 C. 
At the conclusion of the heating period the remaining liquid 
was dissolved in iso-octane for quantitative spectroscopic 
analysis to determine whether any NPN had decomposed. 
The experimental results are shown in Table 3. These 
analyses were performed using a pair of matched cells at a 
wave length of 11.600u and a slit width of 3994. Study of the 
spectra of pure NPN and of NPN mixed with possible reac- 
tion products indicated that this wave length is that most 
suitable. 

The relation between optical density and concentration had 
been determined. It follows Beer’s law, i.e., is linear, at con- 
centrations corresponding to less than 100 per cent of the 
original NPN present in the solvent, but at higher levels the 
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Thermal Stability of NPN in the Presence of 
Various Additives 
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optical density becomes quite insensitive to the concentration. 
Therefore since there is no mechanism here for NPN produc- 
tion, all analyses over 100 per cent are reported as 100 per cent 
with the actual results given in parentheses. Results indicat- 
ing decomposition, i.e., less than 100 per cent NPN, are deemed 
reliable. 


Discussion 


The conditions indicated in Table 3 are quite rigorous. 
Taking the activation energy for NPN as 36.9 kcal/mole (3), 
6 days at 100 C, for example, correspond to about a year at 
70 C (160 F). Of the runs 4-10 of Table 3 in which metal ad- 
ditives were used, none gave evidence of significant destruc- 
tion of NPN, the variations in the NPN analyses being 
within the error of the measurement. Run 9 with cold rolled 
steel gave some evidence of reaction in that the sample 
showed a gain in weight of 1-2 per cent. 

Runs 11-14, in which chemicals which might be formed in 
NPN by its decomposition were the additives, showed strik- 
ing results especially with nitrogen dioxide and acetaldehyde. 
Water and acetic acid likewise show a deleterious, though 
much less marked, effect. It is interesting to point out that 
although nitrogen dioxide inhibits the gas phase decomposi- 
tion of NPN (1) it accelerates that in the liquid phase; thus 
the difference in phase does affect the mechanism of the reac- 
tion. 
Conclusions and Recommendations 

After a careful review of the data obtained during the course 


Hard copper 


of this investigation some conclusions may be drawn. 
1 The following materials are compatible with NPN and 
K Monel 


can be used as containers: 
303 stainless steel 


1020 cold rolled steel (for periods not to exceed one year) 


ST aluminum alloy 
2S soft aluminum Sexcellent 


2 In containers of the above materials the experimental 
results predict a storage life for NPN, in these containers, of at 
least a year at 160 F. 

Detailed discussions of the investigations reported herein 
including the spectra referred to may be found in (4,5), 
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Table 3 Thermal stability of NPN in the presence of various additives 
Heating owe 
period, 
days Pellet weight, gm Per cent NPN Com- 

Run Material 100 C Before After remaining patibility Remarks 

1 pure NPN a 6 100 a3 control 

2 pure NPN 13 96 ters control 

3 pure NPN | 13 ~100 (>110) 

4 hard copper — 6 0.8935 0.8935 96 excellent oe 

5 ST Alalloy ni) 0.2720 0.2719 98 excellent Sats 

6 2S soft Al 13 0.2706 0.2704 ~100 (> 110) excellent 

7 KMonl 13 0.9082 0.9083 ~100(>110) excellent 

8 303 stainless steel 6 0.7775 0 95 excellent en) 

9 1020 cold rolled steel 6 0.7800 0.7 96 poor > an 
10 Armco cast iron 13 0.7823 0.7822 ~100 (100) good? 
11 saturated with 6 1 ml saturated liquid 

used 
12 trace CH;CHO 6 : 20 says CH;CHO added to 1 ml 
NPN 
13 saturated with water 6 81 a9 1 ml saturated liquid 
used 

14 ml acetic acid 6 91 “OS 

Pellet showed rust colored discoloration. 
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Refraction Correction Techniques in 
Earth Satellite Optical Tracking for — 


Precision Prediction’ 


ii Aircraft Corp., Missiles and Space Division, 
Sunnyvale, Calif. 


An analysis to account for atmospheric refraction in 
photographic satellite tracking has been made. Assum- 
ing star background satellite photographs, yielding right 
ascension and declination of occulted bodies, a definite 
integral formula is developed for determining ‘“‘central 
zenith distances’’ of a satellite orbiting a body having ap- 
preciable atmosphere (comparable to that of Earth). 
A first approximation for evaluation of the formula is de- 
tailed. Next, a method for constructing a refraction table 
using culmination techniques is outlined. Finally, a 
possible high precision combination optical-radar tracking 
and prediction system is suggested, and some of the 
problems caused by an assumption that Earth is ellipsoidal 
are indicated. 


ROM the beginning of Earth satellite programs, opening 
with the IGY satellites, designers and systems engineers 
have favored electronic over optical tracking. At the same 
time, however, greater and greater accuracy in determining 
position and motion for acquisition and especially for ephe- 
merides construction, has been demanded. When the major 
problem in optical systems, that of processing time-lag, 
is alleviated, and the ultimate system selection is made, it 
is certain that atmospheric light refraction will be accounted 
for in the optical systems. It is the aim of this paper to 
present a first solution to the refraction problem. 


Determination and Incorporation of Central Zenith 
Distance Into a Prediction Program 


It is assumed that a number of satellite intercepts pro- 
viding the right ascension r.a. of the occulted star and its 
declination 6 along with time ¢ are available, say from star 
background photograph readings. (See Fig. 1 for these as 
well as other relevant quantities.) The polar angle 0, of the 
satellite, relative to the zenith direction as polar direction, 
will be called the “central zenith distance” of the satellite 
in the remainder of this paper. The tracking station is 
located at T. Apparent right ascension and declination r.a.,4 
and 6, are shown for the convenience of the reader. 

The following paragraphs outline a procedure for deter- 
mining the true direction of the satellite after correcting for 
refraction. Range, or more directly, vectorial radius dis- 
tance, is determined simultaneously. 

The first step is to transform the r.a. and 6 values into corre- 
sponding values of azimuth A, and elevation angle ay. 
This might be done by computing inertial direction numbers, 
“rotating” to Earth oriented direction numbers with the aid 
of the Greenwich hour angles of the inertial reference point, 
and then “rotating” to spherical Earth tracking station hori- 
zon orientation. A suitable reference point is “the first point 
of Aries.” Tables of hour-angle values of this quantity are 
available (1).* Corrections for Earth oblateness could be 
made; but to avoid complexity a spherical Earth is assumed 
here. 

The second step is to compute {o, the apparent zenith dis- 
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Fig. 1 Atmospheric light refraction geometry 


tance of the satellite (or the occulted star). This can be 
done by solving iteratively 
fo = — [awe + 
where R({) is given (2) as 
= 58.294” tan [> — 0.0668” tan? (2] 


an empirical formula whose coefficients were derived from 
observations of stars. 

Under the assumption of a spherical Earth, it can be shown 
(2) that 


ur sin = constant [3] 


_ u is the index of atmospheric refraction of light and is 
a function only of the radius vector distance to the satellite, 


and ¢ is the angle between the radius vector and the tangent. 
Now, by definition 


= [4] 
is the speed of re in a 


(rsin {)/v =C [5] 
By substituting values of the variables at Earth’s surface 
(denoted by zero subscripts) into Equation [5], the value of 
C is found to be R sin {o/%, where R denotes the radius of 
Earth. Values of these quantities are given in recent pub- 
lications (3). From polar coordinate representation 


tan ¢ = r/(dr/d6,) 
yielding 
sin = (r d0,)/~/ (dr)? + (rd0,)* (7 


Combination of Equations [5 and 7] and simplification, yield 


Here v is the speed of light and v, 
vacuum. Thus 


Cv 
[8] 
712,.2 
6, = cf a [9] 


From (3) — 
{1+ ((u — 1)/71273.16} (5) 
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where 7 is the temperature in deg K and the bars indicate 
values at zero deg C and 760 mm of mercury. A recent 
representation of the atmospheric density p is given in (4) 
as 


p = constant X {[aR(r — R)/r] — b}~ [11] 


where a, 6 and c are constants. Thus, the integrand in 
Equation [9] is now expressed as a function of r and the deter- 
mination of 6, depends on the evaluation of a definite in- 
tegral. This evaluation can surely be done satisfactorily on 
the computing machines in conjunction with a prediction 
program. Several approximation techniques, however, will 
be given in later sections. 

Assuming 09, has now been calculated for the predicted 
values of r at the times of the intercepts under consideration, 
the elevation angles a can be determined from 


tan a = (rcos 6, — R)/(r sin 82) {12] 


Then, from a number of pairs of values of Ax and a, corre- 
sponding values of slant range S can be found by triangula- 
tion (in the cases of simultaneous intercepts from two dif- 
ferent tracking stations); by direction only methods (5); by 
direction plus Doppler shift data methods (5,6) when supple- 
mental Doppler data are available; or by two-station mini- 
mum distance methods (7). Then improved values of r can 
be determined from 


+ S sin a)? + S? cos? [13] 


At this point an iteration should be performed by repeatedly 
substituting the value of r computed in Equation [13] into 
Equation [9] to obtain a new 8@,, and then substituting into 
Equation [12], calculating S by one of the processes called 
out previously and determining new values of r by Equation 
[13]. The limiting values of a and S, along with corre- 
sponding values of Ay, would then be available for insertion 
into a basic prediction program. 


First Approximation to the Central Zenith Distance 


Substitution of the value of C into Equation [9], using the 
definition of index of refraction given by Equation [4], and 
simplification yield 


r dr 
6, = Rsin — R* sin? [14] 


For terrestrial space it is reasonable to replace (u/uo) by 
1 + &, where is a small positive number, since w and po are 
nearly equal in that case. Then Equation [14] can be 
written approximately as 


dr 

6, = R s1n Sf r(r? sin? 
r dr 


Following 


neglecting powers of & that are higher than 1. 
Smart (2), Equation [15] can be further reduced to 


Resin fo sin r 
6, = — arc sin i, Edr [16] 


If the function consisting of the first two terms of r.m. of 
Equation [16] is expanded about n = 0, where R/r = 1 — 7 
(with 7 a small positive number, and the second and higher 
powers of » neglected) and & is replaced by [(u/uo) — 1], 
then Equation [16] becomes 


R sin r 


If « is replaced in Equation [17] by 1 + kp, following the 
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“Law of Gladstone and Dale,” Equation [17] becomes rt 


0, = (: - 2) tan [> — tan (o(1 + tan? X 


1 — r—R 1 


which has the form » 
6, = A(r) tan + B(r) tan? [19] 


The integral fife dr appearing in B(r) can be written in 
closed form in the following expression for p 
p = K(a + b/r)-™ 7 [20] 


where a, 6 and X are constants (4). The integral in question 
can be integrated by first making the transformation r = 1/v 
and then using a reduction formula twice to get the power 
form 


d 21 
(2 + dv [21] 


Thus, Equation [19] represents a closed form first approxi- 
mation for the central zenith distance as a function of radial 
distance and apparent zenith distance. 


Culmination Method of Central Zenith Distance Table 
Construction 


For situations where table look-up is favored over real- 
time calculation, or for checks on calculations made by the 
method just discussed, it might be desirable to have at one’s 
disposal a table of values of 6, against r and {. This could 
be done in the following way. For an existing, well-known 
satellite with nearly circular orbit, determine the latitudes 
@ and @’ and zenith distances {) and {)’ of two culminations 
(crossings of the local hour circle) of a selected tracking sta- 
tion, say one north and one south. Substitute these values 
into the following equations 


6, =L — 0 = A(r) tan + Bir) tan? 
A(r) tan £0’ + B(r) tan’ Ce 


To derive these equations, see Fig..2. The latitude of the 
tracking station at 7’ is denoted by L. 

Simultaneous solution of Equations [22], using the ap- 
propriate tracking station latitude, will yield A(r) and B(r) 
for the measured {o, and r, 8, 6’, assumed to be known ac- 
curately. Thus, a table of central zenith distances can be 
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constructed for various values of {, obtained from different 
sightings, for a particular radius vector r. Repetition of this 
procedure for a range of values of r would produce a complete 
table. 


If a tracking system provides accurate range values S, say 
from radar equipment, as well as accurate values of r.a. and 6 
of an occulted celestial body, then there is an alternate to the 
procedure explained previously. As before, compute 9, 
from Equation [9] using the predicted value of r. Then use 
Sand 


r? — 2Rr cos 0, + R? — S* = 0 [23] 


to determine a new r, and iterate. It might be that a tracking 
system approaching the ultimate in precision would use radar 
range and optical directional values in this manner as its 
basic input quantities. 

Note further that the iterative process using Equation [1] 
for determining {> is precluded if a theodolite is used in con- 
junction with the camera in obtaining the star background 
photographs. The zenith distance {) is measured directly 
by the theodolite in this case. 

Since Earth is nearly ellipsoidal, further accuracy might be 
obtained if the fundamental relation expressed by Equation 


[3] were to be replaced by a similar relation valid under the 
ellipsoidal assumption. This would require that the inertia! 
direction numbers be transformed to ellipsoidal Earth-track- 
ing station orientation, and that Equations [12 and 13] be 
replaced by equations accounting for oblateness in the 
geometry. 

Finally it should be pointed out that, although the analysis 
for determining central zenith distance was done for applica- 
tion to refraction of light through the atmosphere, it should 
be equally satisfactory for refraction of any electromagnetic 
energy transmission through a medium homogeneous ir 
spherical shell lamina. 
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7 Iterative Solution of the N-Body 
+ Problem for Real Time 
L. M. RAUCH! 


General Dynamics Corp., San Diego, Calif. : 


HIS paper deals formally with the iterative process in the 

solution of the noncolliding N-body problem. The inde- 
pendent parameter is no longer the potent pseudo-time (3)? 
but the real time. 

The process makes use of the Steffensen (4,5) coefficients. 
These coefficients were introduced by Steffensen for the three- 
body problem and extended by Cheney (1) to N bodies. 

The iterative process, as formulated in the paper, has a 
number of features worth noting: 

1 Itis an alternative process in the solution of the N-body 
problem. 

2 The sequence of iterative solutions converges rapidly 
to the limiting solution of the system of differential equations. 

3 The process, as such, is nonrecursive, so that any suc- 
cessive approximation is independent of any of the preceding 
determinations. 

4 A large digital computer (say, the IBM 704) is not 
necessary in the numeric determinations. 

The paper is divided into three parts. The first part formu- 
lates the iterative process as applied to the system of equa- 
tions of motion for the N bodies. This is followed by de- 
ductions of the explicit solutions. The concluding part gives 
a summary of the inputs and equations in numerization of 
the N-body problem. 


Iterative Process 


This part deals with the equations of motion of N bodies in 
a Newtonian potential field (6) and the transformation of 


these equations into iterative form. 7 ee 
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Equations of motion - 


The n given quantities (M;) are assumed to be point 
masses. Their position vectors z‘ are given by the com- 
ponents z* in an inertial space. It is specified that the 


range for the super or subscripts is given as 


The magnitude F;,; of the attractive force exerted by the 


mass on M; is given by 
F;; = GM /r* - 7 


where r*/, the square of the magnitude of the position vector 
between the 7 and j particles, is given by 
= (xh — gin)? 


The direction of the force vector is along the unit vector 
(rii)~/2(zi — zx‘), Thus, the total attractive force on the 
mass M; is 


SH = (rit) [2] 


The application of the second law of motion leads to 7 

# = 2; (xi — [3] 
which is a set of 3 n differential equations of the second order. 


Iterative system of differential equations 
To specify an iterative system of Equation [3], the follow- 
ing operations are performed: Express the quantities in 
Equation [2] and the solution x of Equation [3] in power 
series 
zh = > ri = prtit® Si = [4] 
k=0 k=0 k=0 
Next, form the finite sums 


a+s 
k=0 k=0 k=0 


where @ is an arbitrary fixed integer > 0, a = 1,2,3,.... 
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if 
Finally, write the approximate system of Equation [3], cursive formulas [14 and 15] may be reduced to a nonrecursive 
namely form. Itis given as 


= 20; — = 1,2,3,... [6] 


The system [6] differs from [3] in the sense that the solu- 
tions x of [3] are given by the solutions z.** by means of a 
limit, namely 

limit = x** (7] 
as 
for any th; a = 1,2,3,.... 

The approximate system turns into an iterative one by a 
substitution of (a — 1) for a in the right member [6]. Thus, 
the iterative system of differentia] equations is defined as 


= (28-1 — 2-1) [8] 17] 

where Equation [17] is a linear relation between 2,** and 29** for 
= + €qi* [9] any The quantity is given in (3) as 

limit = 0 [18] 


The function €q" is defined by Equation [9]. However, in 
(3) it is shown that 


@ 


for any 7,h and ¢ in a defined interval. Equation [17] is an 


mavee cat) = 0 [10] iterative, nonrecursive approximation of the system of the 

equations of motion [6] provided a meaning may be attached 

for any 2,h and t in a defined interval. It will be assumed in to the initial iterative function zo and if the quantities A and 

this paper that Equation [10] is valid. C are known functions. The part that follows will establish 
Equate the right members of the two types of systems [6,8], a meaning to the above mentioned quantities (70,A,C). 


so that 
Explicit Formulation of the Iterative Solution 
To explicitly formulate the iterative solution given by 
Equation [17], a meaning must be attached to the initial 
iterative function 2 and the functional quantities Aa’, 
C.'*. We do so by first expressing these quantities in terms 


— 23-1) + (ea*/G) [11] 


The formal solution x2 of the approximate system is found 
by solving Equation [11], namely 


M;Sa-1 + 2 — Sa-1Ta-1)_ [12] 


Lai? = 


Nonrecursive form of the iterative solution 
of the definitions [5] and then evaluating the coefficients in 


We now reduce the formula [12] to a nonrecursive form. aa 
these definitions. 


Define the following symbols 
Expressions for the quantities in the iterative formula 


By means of the first definition in Equations [5 and 13], 
Ca = 24j — Sa-1%a-1) the quantity Aq‘ becomes 
| Ad = @ 2, [19] 


The symbol C,**, because of Equations [5 and 13], becomes 


a+s a+s a+s-1 a+s-1 


k=0 1=0 k=0 1=0 
a+ a+s-1 a+p-1 a+p-1 a+p-1 
1=0 k=0 k=0 k=0 =0 
so that 
a+s—1 
The symbolic form of Equation [12] is given as = 12+8 + Bas 
We define the function zo by the assumption that the quan- 
where tity za** in Equation [5] is to be valid fora = 0. Thus 
= [15] B 
It may be shown by mathematical induction that the re- Gant : 


| 
II A a-a 
a=0 
a= 1,2,... [16] 
which may be reduced to a 
Act a-q €a—q 
Aa Aa Aa se 
i 
> 


Explicit formulation of the functions A, C and xo 

To formulate the iterative expression [17], the functions of 
time A,‘(t), Ca‘*(t) and x*(t) of Equations [19-21] must be 
well-defined. This implies a knowledge of the coefficients 
£, o (and p) of the series [4]. Fortunately, the investiga- 
tions by Steffensen (4,5) and Cheney (1) allow an immediate 
determination. The paper by Cheney gives the following 
recursion formulas 


k = 0,1,2, . [22] 
k-—2 
= —(2kpot)— + Ova + 
v=0 
3(v + 1) 
v=0 
k—-2 
~()=0 k=1 
v=0 
k-1 


These may be reduced to the following forms ee 


=), (Ee — + — 
so that 


k 


k-2 
oxi = | — v — Lop 
v=0 
3(v + + 3k 

v=0 


so that 


~The quantity 
[29] 
is found from Equations [2 and 4] for ¢ = 0. Subsequently 


the formulas [25-27] are used cyclically with k = 1,2, 
consecutively, to determine the three sets of quantities 


(pei, ont, E41). 


oot = 


= 


Summary of the Method 


The existence theorem for the iterative solution [17] of 
the system of Equations [3] is given in (3). We restate it 
in the following form: 


Theorem: The equality 


4 holds for any 7,h and t in a defined region R where 


La** = iterative solution for any a = 1, 2,... and 
given by Equation [17] 

limit ¢. = 0 

a 

x(t) = set of solutions for the given system 


The region R (or interval for real ¢) is defined, briefly, as the 

open circular subregion in the ¢ plane, whose boundary in- 

- eludes that singularity whose modulus has the least value 
of every singular point (2) given in the right-hand member 
of Equation [3]. 

The procedure in the determination of any one set of iter- 
ative solutions is then as follows: 

1 Choose any reasonable value in the sequence a 
1,2,3,.. . for formula [17]. 

2 Choose any value for 6 > 0. A too large value for B 
would imply the determination of more coefficients from the 
formulas [25-27], since the initial function xo" is specified by 


8 
2, 
k=0 


Find the values for the sets (p;*7, &**),... 
£a+8)from the recursion formulas [25-27]. 
4 The quantities &*, & (which determine all subse- 
quent coefficients) are the arbitrary constants determined by 

_ the initial conditions x**(0), z**(0) and found by the equalities 


x(t) = 


(pa+6-1) 


tes = 0, namely nay 
er oe - § Substitute the coefficients computed in step 3 into the 
( formulas [19-21]. 
ie [26] 6 The values thus found are then incorporated into Equa- 
k-1 
- ih — 
(k + k) OM ) [27] References 
1 Cheney, W., ‘‘Power Series for the N-Body Problem,’’ Mathematical 


Slight consideration of the above three recursion formulas 
indicates that the evaluation of p and o and so & depends 
on a priori knowledge of £** and &**. But these 6 n quanti- 
ties specify the initial position and velocity of each of the 
n masses. 

With this information given (£*, &**), the value of po‘ 
is computed from Equation [25], namely 


— 
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Considerations in 


W. E. JAHSMAN! "4 


. Lockheed Aircraft Corp., Palo Alto, Calif. 7 


N A RECENT article (1),? Tsu has described how a solar 
sail may be used for interplanetary travel. He has dis- 
cussed the several advantages of this device: Mass ratio of 
unity; no fuel, propellant or propulsive powerplant required; 
and no waste-heat disposal problem. He has also mentioned 
certain problem areas in solar sail design: Need for a sail ma- 
terial with low mass per unit area and high reflectivity; 
fabrication, handling and packaging techniques for the sail; 
and in-flight maneuvering and control. The purpose of this 
note is to point out another equally important problem when 
a ring is used to support the sail: The contribution of the 
ring mass to the overall system mass. 

The necessity of a supporting ring is obvious; without one, 
the sail becomes a parachute, and the requirement of low sail 
mass per unit area (for a given thrust) is not satisfied. Maneu- 
vering and control problems also become much more serious. 
On the other hand, the ring dimensions must be chosen so that 
strength and/or stability limitations are not exceeded. The 
following equations of motion for the sail oriented normal to 
the sun’s rays are formulated with these considerations in 


mind. (Tsu’s nomenclature is used wherever possible.) 
porra*f(y) — 2raP sin y = ma 
sin y — nT cos 6 = ma [2] 
nT cos 6 = m,a [3] 
m =pmatgy) [4] 
Mr = pr2mah? [5a] 
rong (= [stability limit, 
Pcosy + ) 4 \a/_— see (2)] [6a] 


Oar h (‘) (strength limit) 
a 


[7] 


The reader is referred to Fig. 1 and the Nomenclature for a 
description of the symbols. 

If one is interested in optimization of one of the variables, 
say a, expressions [1 to 7] may be combined to a single rela- 
tion between a and a, 6 and y. (The payload mass m, and 
the physical and material properties p, Z and o,1 are taken 
to be fixed parameters.) Maximum acceleration within speci- 
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Fig. 1 Solar sail components 
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Ring Supported Solar Sails — 


fied ranges of a, 6 and y may then be found. Here, however, 
it will be sufficient to consider an example given by Tsu to 
demonstrate the effect of the ring. First it should be noted 
that to obtain a lower bound on the ring mass, only expres- 
sions [2, 3, 6a and 5a] need be studied. Adding Equations 
[2 and 3] gives 


(m, + mp)a@ 
4 27a sin [8] 
Y 


If 6 is sufficiently small, the second term on the left-hand side 
of expression [6a] may be neglected with respect to the first. 
This restriction is equivalent to lengthening the shroud lines, 
which is permissible as long as the assumption of massless 
shroud lines can be satisfied. Then, by substituting for P in 
expression [6a] from Equation [8], one obtains a relation for 


the ring depth* h 
Finally, an expression for m, the ring mass is found by sub- 
stituting for h in Equation [5a] from expression [9] 


+ m,)a cot (10a 
TE, 


Expression [10a] shows that high modulus, low density ring 
material will reduce the ring mass requirements; hence, 
beryllium (EZ, = 45 X 108 psi = 3.1 X 10!? dynes/cm? and 
p = 1.85 gm/cm') is chosen in this example. Also, since + 
should be kept small from sail mass considerations, cot ¥ is 
taken as 10 (y ~ 5.8 deg). Other data supplied by Tsu are‘ 


2a*(m, + mp)a cot 
TE, 


p. = 1.18 gm/cm* 
a = 2.5 X 104 cm 
a = 0.2 cm/sec? 
t, = 1.72 X 10-* cm 
m, = 4 X 10° gm 
m, + mp = 5 X 10° gm 


Substituting the appropriate values into expression [10a] 
gives m, > 33.5 X 105 gm, an impossible result since m, < m; + 
mp, = 5 X 105 gm. Referring to expression [10a], we see 
that in order to get a physically meaningful answer, we must 
relax the requirement on y (allowing the sail to become more 


The stability limitation will undoubtedly govern in most 
structures of this type. Dimensions and loads which are ob- 
tained using this criterion may subsequently be substituted into 
the latter of expressions [6a] for the ring and expression [7] for 
the sail, to check that strength limitations are not exceeded. 

‘It should be pointed out that depositing a “few microns’ 
(1u = 10-4 cm) of aluminum on a plastic base to achieve high 
reflectivity more than doubles the sail thickness (hence, mass) in 
this example. Development of a highly reflective plastic is 
clearly desirable. 


(A) STIFFENING SPOKES (B) HOLLOW TORUS 
Fig. 2 Methods for reducing supporting ring mass 
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parachute shaped the ‘mass), or 
reduce a (thereby accepting a considerably smaller accelera- 
tion). 

Two methods of reducing the ring mass come immediately 
to mind—adding spokes and replacing the solid ring by a 
hollow torus. These two methods are illustrated in Fig. 2. 

The effect of the number of equally spaced spokes 7 is to 
introduce a factor of (t2 — 1)/3 to the denominator of ex- 
pression [10a]. For the ring illustrated in Fig. 2, 7 is 16 and 
m, is greater than 3.6 X 10° gm. Thus the use of spokes will 
certainly aid in reducing the ring mass; however, as 7 be- 
comes large, the spoke mass will become significant. 

A much greater reduction of ring mass can be achieved by 
replacing the ring with a solid cross section by a hollow torus. 
If the torus wall is sufficiently thin, pressurization may be 
necessary to maintain circularity of the cross section. Then 
the ring mass is given by 


= + 2p ab? [5b] 


(see Fig. 2b). Furthermore, by assuming that torus strength 
is determined either by stability or by noncompressive wall 


stress conditions, we can replace expression [6a] by anes _ 


3 
T sin 8 fore (°) 
[6b] 


apb? 
a 
Combining [2, 3, 5b and 6b] as before prov ides an expression 
for m,’ i. 
a*(m,’ + mp)a cot 
> 4 2 t, 
m | + 
a*(m,’ + mp)a cot 
2rp,a ] [10b] 


In this calculation, the torus is assumed to be fabricated from 
plastic film. From a practical standpoint, this design has the 
advantages that it may be packaged during ferrying and easily 
erected by inflation in orbit. In addition subsequent calcula- 
tions will show that there is relatively little improvement in 
mass reduction by using beryllium. Hence, when the ring 
material and thickness are the same as those of the sail, it is 
found that m,’ > 5 X 10*gm when E,’ = 5 X 105 psi = 3.45 X 

10° dynes/em*. The mass of gas and ring material (first 
and second terms on the right-hand side of expression [10b]) 
are about equal because of the large gas volume. The density 
of the gas (air) was found using the perfect gas law with a 
temperature of 100 F and a pressure governed by o’a1,7, the 


allowable stress in the torus wall aa 

In this example, o’,1,, was taken to be 2500 psi = 1.72 X 105 


dynes/cm?. 

The value 5 X 10° gm is only 1 per cent of the ring and pay- 
load mass. If beryllium is used, this figure is reduced to about 
0.33 per cent. Thus, the ring mass can be made a very smal! 
fraction of the solar sail system mass if a hollow, pressurized 
torus design is used. On the other hand, the problem o! 
micrometeorite puncture and subsequent loss of stiffness can- 
not be accurately solved until more precise data on distribu- 
tion of matter in space is obtained. Furthermore, the exact 
theory of pressurized membranes has not yet been developed; 
expression [6b], which is a rough approximation at best, ma) 
be superseded later on by more exact formulas. Hence, theré 
is still considerable work to be done in the design of solar sails 
before they can ever be considered as a feasible means of in- 
terplanetary travel. 


Nomenclature 
= Young’s modulus 
= stress resultant at the edge of the sail 
= tension in the shroud lines 
= radius of the sail (and ring) 
b = radius of the torus cross section 
f(y) = function representing dependence on y of total thrust 
due to radiation 
g(y) = function representing dependence of sail mass on 
h = width and depth of ring cross section 
i = number of stiffening spokes me 
n = number of shroud lines 
= internal pressure in torus 
Po = radiation pressure 
t = ring or torus wall thickness 
a = acceleration of solar sail system 
6 = shroud line angle 
Y = sail edge angle 
= Poisson’s ratio 
p = mass density 
gait = allowable stress 
Subscripts 
p = payload 
r = ring (or torus) = 
8 = sail 
= refers to torus 
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Expansion in 


Equilibrium Air' 


HAIM KENNET? 


Massachusetts Institute of Technology, Cambridge, Mass. 


RANDTL-MEYER flow for air in dissociative equi- 
librium is encountered in both external and internal fluid 
dynamics, particularly in hypersonic flow and rocket engine 
nozzle flow. The problem has been considered by many 
authors, e.g., Feldman (1),? Heims (2), Hayes and Probstein 
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(3), and the solution usually advocated is a numerical inte- 


a gration of the appropriate relations along an isentrope. It is 
_ intended to show here that by making use of some simple 
_ approximations, which are valid for equilibrium air in hyper- 


sonic flow, a closed form expression can be obtained relating 
the pressure p to the angle of turn 0. 

The basic equation for Prandtl-Meyer flow can be obtained 
by combining the momentum and global continuity equations 
which, for an inviscid gas, do not show explicitly that the gas 
is in a dissociated state. Hence, the well-known relation 


d 
[1] 
M?-1 
is still valid. Making use of the energy equation, which re- 
lates the velocity g to the enthalpy h of the gas, namely 
q? Us 
g tho [2] 
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where A; is the stagnation enthalpy, and the hypersonic ap- 
proximation, valid for large Mach number M 
VM?-1=M 


= q/a [3] 


an expression relating the angle of turn to the enthalpy of the | 


gas and its speed of sound a is found 


Dimensionless flow variables are introduced at this point to 


[4] 


nform to Feldman’s Mollier chart notation. These di- 
mensionless quantities 
h 
ho Po 80 


are given in terms of the following reference quantities 


ho = RT, = 84.75 X 104 ft?/sec? 

a = 1089 fps 

po = 1 atm 
& = R = 6.886 X 10-? Btu/bR 


The constant in front of the integral sign in Equation [4] is 
merely the result of this particular choice of the reference 


quantities 
V = 0.8454 [6] 


For constant entropy s* = s/R the Mollier chart shows that 
* and a* vary linearly with logip*, provided p* does not 
vary by more than one order of magnitude. For larger varia- 
tions in p*, the relation between a*, h* and logiop* is piece- 
wise linear. The following analysis is nevertheless applic- 


> 
h* = a + az p* 
a* = be logio p* 


where a), 2, b;, be are entropy dependent constants to be de- 
termined from the Mollier chart. When Equations [7] are 
substituted into Equation [4], integration and further manipu- 
lations yield 


-tanh?(C+ Dé) 


A = 
2 a2 
| 1 logie p 
C =tanh- | b; (is —\h.* bi 
V + b? a2 a2 
be 
= V — 2a + 2ax(bs/b2) 


To demonstrate the effects of large amounts of dissociation, 
an expansion along a geodesic on the surface of a nonlifting 
circular ogive, with a cone angle 6, = 37 deg, flying at M. 
20 at an altitude of 130,000 ft has been calculated with and 


without dissociation. At 130,000 ft 
Po = 3.5 X 10-%atm 
a. = 1.1 X 10° fps 
h,/RT) = 3.5 
T. = 273.16K 


_ The conical shock relations for inviscid flow were taken from 


_ (4), whereas those for dissociated air were taken from Romig 
(5). For this particular case the entropy behind the shock for 
the dissociated stream is s:.* = 39. Hence, from the Mollier 
chart 


a, = 123.707 

a2 = 43.103 

bh = 4.282 
bo = 0.877 


For an expansion of @. = —10 deg, the result for the pressure 


is: With dissociation, p, = 0.193 atm, without dissociation, 
p2 = 0.371 atm. 
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“‘Hypersonic Flow Theory, 


Convection in a Fluid at Supercritical 
Pressures 


J. D. GRIFFITH! and R. H. SABERSKY?’ 


California Institute of Technology, Pasadena, Calif. 


EAT transfer to fluids near their critical point is becom- 

ing increasingly important in connection with such ap- 
plications as the cooling of rocket motors with hydrocarbon 
fuels or the heating of water in high pressure boilers. Several 
theoretical and experimental papers have been published on 
this subject, e.g. (1-4). The principal characteristic of a 
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fluid near its tiition point is the strong effect of tempera- 


ture and pressure on the density, specific heat and viscosity. 


By taking this variation in properties into account, the heat 
transfer rates for several fluids and in a fairly wide range near 
the critical point have been predicted successfully. For some 
fluids, however, very large increases in heat transfer coefficient 
have been observed which cannot be explained on the basis of 
property variations (2). It has been suggested that these in- 
creases might be caused by the formation and collapse of low 
density, bubble-like aggregates near the heating surface (2,4). 
These aggregates, although consisting of the same phase as the 
surrounding fluid, are believed to behave in a manner very 
similar to vapor bubbles. The improvement in convection 
is therefore attributed to the agitation caused by the motion 
of the aggregates. 

The purpose of the present set of experiments was to visually 
observe the fluid motion in the neighborhood of the heating 
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Fig. 2 Heat transfer rate vs. wire temperature for various fluid 
pressures and bulk temperatures 


surface, to determine if any bubble-like motion can be de- 
tected, and if so, to relate this motion to any sharp increases 
in the heat transfer rate. 


Test Equipment 


An electrically heated wire was selected as the heating sur- 
face. The wire was made part of a Wheatstone bridge. 
Changes in resistance of the wire were measured by the un- 
balance of the bridge. From this unbalance the wire tempera- 
ture was then computed. The test wire was enclosed in a 
heavy steel cylinder, the ends of which were provided with 
thick quartz windows. — A schematic sketch of the cylinder is 
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shown in Fig. 1. The cylinder was filled with the test fluid 
and the fluid was pressurized by means of a high pressure 
nitrogen system. The wire and surrounding region were 
illuminated by a Schlieren system. The motion of the fluid 
could be observed on the ground glass of a camera and photo- 
graphs of the flow picture could be taken. 


Experimental Results 


For the present experiments Freon 114A was used as a test 
fluid. It was selected mainly because of its relatively low 
critical pressure (per ~ 478 psia) and temperature (T., ~ 
294 F). The wire temperature was then measured as a fune- 
tion of the heat transfer rate for several values of pressure, 
both below and above the critical pressure. Photographs were 
taken at most of the points at which measurements were 
made. 

The results are summarized in a graph (Fig. 2), showing tl.e 
heat transfer rate as a function of the wire temperature fcr 
various fluid pressures. The first two curves, labeled A and /, 
were taken at pressures below the critical. These curves show 
the characteristic rapid rise in heat transfer rate with wal 
temperature, which coincides with the onset of nucleate boi - 
ing. The corresponding photographs clearly show bubble-. 
One of these, which was taken at the point marked “1” i: 
Fig. 2, is shown in Fig. 3. The data for curves C through // 
were taken at supercritical pressure. An examination of th: 
photographs corresponding to these curves as well as visual 
observation in this region, did show bubble-like activity i: 
some instances, as had been hypothesized. This phenomeno: 
was observed principally at relatively high temperature dii- 
ferences, the conditions indicated by point 2 in Fig. 2 being 
typical. A photograph taken at this point is reproduced in 


Fig. 3 Bubble formation in the nucleate boiling range at a pres- 
sure of 350 psi and a bulk temperature of 80F 


Fig. 4 Bubble formation in the supercritical range at a pressure 
f 550 psi and a bulk temperature of 80F 
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Fig. 4. However, the bubble-like activity that was observed 
did not result in any observable change in the slope of the heat 
transfer curve. The activity was apparently not sufficiently 
large to overshadow the effects of the normal convection. 
Nevertheless the observations do support the view that 
bubble-like activity might exist at supercritical pressures. 


Conclusion 


In the foregoing a simple experiment has been described for 
the study of free convective flow patterns in a fluid at pres- 
sures beyond the critical. The visual results give some indica- 
tion that bubble-like aggregates may form under these con- 


in nucleate boiling. However, this result should not be in- 
terpreted to mean that such activity may not significantly 
contribute to heat convection in a different fluid or under dif- 
ferent conditions of test. 

It is planned to continue the experiments and to include 
several other test fluids. 
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Auxiliary carriage for a missile (2,- 
899,898). W. H. Goss, Silver Spring, 
Md., assignor to the U. 8S. Navy. 

Carriage comprising a number of 
auxiliary rockets secured to each other 
at their forward ends. A thrust plate 
secures them at their rear ends and 
orbits the rear end of the missile. 

Stiffener and stabilizer for supersonic 
vehicles (2,899,899). J. Hirsch, Pacific 
Palisades, Calif. 

First and second fins mounted on one 
of the two missile elements affixed to- 
gether. A third fin mounted on the 
other element is disposed midway between 
the first and second fins. Shock waves of 
equal magnitude are produced between 
the first and second fins when the spaces 
between fins are identical, and of unequal 
magnitude when they are not identical; 
in this case a restoring force to equalize 
the spaces is provided. 

Roll control (2,911,912). A. E. Wetherbee 
Jr., Newington, Conn., assignor to the 
United Aircraft Corp. 


2,011,012 gyro sensor 


Paddle-like member with flat surfaces 

immersed in liquid fuel in a flight vehicle. 
Power means moves the paddles and op- 
poses rolling motion about the longitudinal 
axis. 
Duct balloon (2,900,147). W. F. 
Huch, E. P. Ney and J. R. Winckler, 
in Paul, Minn., assignors to the U. S. 
Navy. 

Free high altitude plastic balloon with 


Epitor’s Note: Patents listed above 
were selected from the Official Gazette of 
the U. S. Patent Office. Printed copies 
of patents may be obtained from the 
Commissioner of Patents, Washington 25, 
D. C., at a cost of 25 cents each; design 
patents, 10 cents. 
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an opening on top. An outside duct 
extends throughout the length of the 
envelope. In flight, gas will escape from 
the duct bottom when gas in the envelope 
tends to be superpressured at a level 
below that of the duct bottom. The dis- 
tance of the lower end of the duct from 
the bottom of the envelope is predeter- 
mined to thereby predetermine the desired 
ceiling altitude of the balloon. 

Combined ramjet and rocket engine 
(2,912,820). Q. R. Whitmore, Wayne, 
Nebr., assignor to the U. S. Navy. 


2,012,820 il 


Engine casing with rocket fuel burning 
means upstream of the discharge nozzle. 
A partition separating the chamber into 
two parts is destructible by the heat 
generated by the burning of the rocket 
fuel, thereby connecting the chambers. 
The engine then functions as a ramjet. 
Ejection seat catapult (2,900,150). G. E. 
Hirt and E. A. Martino, East Haven, 
Conn., assignors to Talco Engineering 
Co., Ine. 

Emergency escape rocket device at- 

tached to the seat of an aircraft. Upon 
initiating a charge into the chamber, the 
rocket releases a propulsive charge, eject- 
ing the seat and occupant from the air- 
craft. 
Rocket fuel and fuel additive (2,900,788). 
R. 8S. Felberg and J. R. Gould (Member 
ARS), Morristown, N. J., assignors to 
Thiokol Chemical Corp. 

Hypergolic propellant formed by mixing 
methylamino, butadiene and nonhyper- 
golic jet-type hydrocarbon fuel, contacting 
the propellant with an equal volume of 
fuming nitric acid. 

Rocket launcher (2,900,874). W. W. 
Tjossem, Burbank, Calif., assignor to 
Lockheed Aircraft Corp. 

Cylindrical turrets containing rocket 

tubes. Under sides of cylinders are 


flattened to correspond with the contour 
of the aircraft body in which the turrets 


George F. McLaughlin, Contributor 


are mounted. Rotating the turrets suc- 
cessively projects the rockets beyond the 
aircraft surface. 
Bridge circuit for 
device (2,902,765). W. 
Segundo, Calif., assignor 
Aircraft Co. 

Incremental displacement detecting sys- 

tem including a grooved scale of magnetic 
material, an electromagnetic detector 
head, and a T-Bridge impedance network 
controlling a circuit. A detector head 
core is disposed in flux linkage with the 
scale, forming a variable reluctance 
magnetic circuit with the scale. 
Container for separate storage of liquid 
propellants (2,902,822). J.D. McKiernan 
(Member ARS), Brooklyn, N. Y., assignor 
to the U. S. Navy. 

Bipropellant rocket with a partition 

having gas ports covered by burst disks 
adjacent the gas producing means. The 
partition and casing form a propellant 
storage space, and a tubular member 
forms a second storage chamber for a 
different propellant. 
Solar tracking system (2,913,583). N. J. 
Regnier and M. R. Shaffer, N. Holly- 
wood, Calif., assignors to Hoffman Elec- 
tronics Corp. 


osition measuring 
T. Chater, El 
to Hughes 


2,013,583 


Electrically intercoupled solar cells 
mounted on a support permitting the cells 
to be tilted. Cells are coupled across 
relay windings controlling batteries power- 
ing motors geared to tilt the support. 
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Electronic monitors for thrust chamber 
(2,902,824). W. A. Sponzilli, 

incinnati, Ohio, assignor to the U. S. 
Air Force. 

Apparatus for monitoring transient 

pressure in a reaction motor. Means 
nsive to pressure consisting of upper 
and lower limit electrical bridge circuits 
containing variable resistances responsive 
to chamber pressure. Motor controls 
connected to the bridge balance detecting 
means actuate when either circuit is in 
balance. 
Pulse jet (2,903,850). T. G. Lang, Los 
Angeles, Calif. 

Underwater propulsion device consisting 
of a variable volume chamber adapted to 
intermittently receive ambient water. 
Gas flows from a gas chamber into the 
water, forcing the water out from the 
chamber. 

et deflector (2,903,851). W. A. Fiedler, 
oodland Hills, Calif., assignor to the 
U.S. Navy. 

Thrust deflector for a supersonic flow 
nozzle having the inner portion of a wedge 
within the nozzle terminating in a sharp 
leading edge, to produce an attached 
shock wave. 

Sighting instrument (2,903,940). J. P. 
Hadjilias, Athens, Greece. 

Prism assembly with two optical prisms, 
each adapted to reflect rays inwardly. 
Rays are directed toward an indicator 
line along which spaced images of an 
object, sighted through the prisms, 


appear. Means for determining the angu- 
lar position of the prism relatively to a 
reference system. 

Means for gees balloon flight 
(2,904,285). . F. Huch, St. Paul, 
Minn., assignor to U. S. Navy. 

High altitude balloon system having a 
predetermined ceiling altitude. The in- 
extensible film envelope has slight perme- 
ability and loss of lift at a predetermined 
rate. The envelope carries a container of 
nongaseous highly evaporative substance 
a compensates for the normal loss of 
ift. 

Radiant energy converter (2,904,612). 
N. J. Regnier, Los Angeles, Calif., assignor 
to Hoffman Electronics Corp. 

Converter cells disposed across one 
face of a planar base in a uniform — 
having active faces parallel to base, 
and directed away from the base in a 
common direction. Annular reflectors 
encircle each cell to reflect radiation onto 
active parts of cells. 

Large area solar energy converter (2,- 
904,613). . E. Paradise, Highland 
Park, IIl., assignor to Hoffman Electronics 


orp. 

Semiconductor particles electrically in- 
sulated from each other, each having two 
regions, one having one direction, and 
the other a second direction. A conduc- 
tive plate is located between the first 
regions, and electrical connections are 
made between the second regions. 


Fuel and rapid ignition apparatus for ram- 


jets and rockets (2,906,094). G. H. 
— and J. Ribovich, McKeesport, 
Combustion unit containing a fuel mass 
with side walls conforming to the casing 
interior, a coating of heat er 
materials on an igniter duct, and a device 
initiating ignition of the coating. 
Jet propelled aircraft with tiltable com- 
bustion chambers (2,912,188). D. EF. 
Singelmann and J. R. Hall, Kenmore, 
N. Y., assignors to Bell Aircraft Corp. 
Primary and auxili combustion 
chambers terminating in thrust nozzles. 
The primary chamber located slightly 
forward of the pitch axis and the auxiliary 
chamber at the aft end of the body. 
Chambers are mounted on pivots to vector 
between horizontal and vertical direction: 


Jet propelled aircraft with jet flaps 
(2,912,189). R. Pouit, Asnieres, Franc: 

Pair of jet engines mounted on a wing. 
Trailing edge nozzles are located along a 
substantial portion of the trailing edge 
A conduit leads from one engine to th> 
other to feed compressed air to a jei 
nozzle on the opposite side of the aircraft 
if one powerplant fails. 


Ear protector (2,910,980). R.D. Stewart, 
Midland, Mich., assignor to the Dow 
Chemical Co. 

Organopolysiloxane rubber protecto: 
conforming precisely to the surface o 
the outer ear and covering the mout! 
of the auditory canal. 


Vol. II, by 
Morton Alperin and H. F. Gregory, 
Pergamon Press, New York, 1959, x + 
318 pp. $15.00. 
Reviewed by RayMonp L. 
University of Southern California 


Vistas in Astronautics, 


The proceedings of the Second Annual 
Astronautics Symposium are presented in 
this volume. The scope is necessarily 
very wide; it ranges from observational 
astronomy to design studies in plasma 
propulsion, from measured data on micro- 
meteorite impacts on an artificial satellite 
to philosophizing on the justification for it 
all. From the mundane to the exotic, 
from science to fiction, the preliminary 
and uncoordinated efforts made by this 
country in the second year of the Space 
Era are reported in five parts: Space 
Environment and Vacuum Research; 
Control and Propulsion of Vehicles Outside 
the Atmosphere; Panel Discussion (on 
utilization of space vehicles); Departure, 
Space Navigation and Re-Entry Problems; 
The Earth’s Moon. The rejuvenation of 
old established subjects, which the scien- 
tific community had laid aside until 
recently, for application to new engineer- 
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Ali Bulent Cambel, Northwestern University, Associate Editor 


ing problems is abundantly evident in 
most of the 25 papers presented. 
Since it presents nothing that is essentially 
new or original, the volume will mainly 
be of historic interest in a fast developing 
field in which hardware construction will, 
at least for some years to come, out- 
distance basic understanding in the 
relevant sciences. 


Missile Engineering Handbook (Vol. IV 
of Principles of Guided Missile Design 
Series) by C. W. Besserer, Van Nostrand 
Co., Princeton, N. J., xi + 600 pp. 
$14.50. 

Reviewed by ALLEN E. Funs 
Space Technology Laboratories 


Perhaps one of your first thoughts 
concerning a missile engineering handbook 
is this: How can material be selected 
which won’t soon be obsolete and yet will 
be useful? The design of a guided 


missile requires, in addition to the “hot 
off the launch pad’’ data, much informa- 
tion of a permanent nature which has 
accumulated in scattered scientific areas. 
For example, the demands of reliability 


have placed a new emphasis on statistics 
and probability. With the advent of 
spaceflight, accurate values of the gravi- 
tational “constant”? as a function of 
altitude take on new significance. There 
are many facts which are of lasting 
importance to the engineer-missileman. 
Besserer’s handbook collects these facts, 
plus some of a more transient character, 
from many sources. 

Engineers who use handbooks may not 
realize it but they acquire certain “tastes 
in tables.’ These preferences arise from 
the fact that there are numerous ways to 
present the same information;  e.g., 
conversion factors can be in tabular form 
like a mileage table on a road map. In 
this handbook conversion factors are 
given in the columnar form; with this 
presentation, errors in the use of the 
table are minimized. In a similar manner, 
the other tables and figures are laid out in 
easily assimilated form. 

Of course, the real value of a handbook is 
the completeness and accuracy of the 
information. These are attributes more 
difficult to evaluate. To give an indica- 
tion of the completeness, there are 
sections on general data, such as conversion 
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Multi-Use 
Automated 
Maintenance 


, The recent demonstration of multi-purpose evaluation equipment, capable of checking 

test equipment (MPTE), developed by a variety of electromechanical devices, . 

: RCA under a series of Army Ordnance con- ranging from radar subassemblies to missile iad 
tracts, highlights a new dimension in auto- guidance computers. MPTE provides the 
mated multi-use systems support and culmi- stimuli, programming, control, measure- . 

! _ nates a long-term RCA effort in this field. ment and test functions for the NIKE AJAX, 

This General Evaluation Equipment is an NIKE HERCULES, LACROSSE, HAWK 

) <j automated, transistorized, dynamic check- and CORPORAL missile systems and has 
gut system. It contains a completely modu- been extended to other weapons systems 

_ larized array of electronic and mechanical related to our defense efforts. 


RADIO CORPORATION of AMERICA 


DEFENSE ELECTRONIC PRODUCTS 
| CAMDEN, NEW JERSEY 
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factors and physical constants, properties 
of the atmosphere, environmental data and 
reliability, properties of materials and 
structures, aerodynamics, avionics, pro- 
pulsion, spaceflight data, miscellaneous 
data, such as selected military specifica- 
tions, design analyses and check lists, 
and finally a glossary of missile terms and 
abbreviations. The glossary is an ex- 
cellent feature of the book; about one 
third of the book is devoted to the glossary 
and abbreviations. It goes from Angstrom 
unit through diergolic propellants, loxo- 
dromic curves, SSIP (subsystem integra- 
tion plan), and ZIP fuels. With the aid of 
this glossary anyone can converse with the 
most starry-eyed missileman and _ the 
many specialists who support him. The 
author has wisely excluded information 
which is readily available in other engi- 
neering handbooks. 

Before discussing the book’s accuracy, 
here is a brief story about two engineers 
who decided to write a book. At the start 
the writers agreed to produce a book 
entirely free from errors. After checking 
the galley proofs, these authors con- 
cluded that their book would have fewer 
errors than any other book in the field. 
As the book went to press the first engineer 
told the second engineer that succeeding 
editions would certainly achieve their goal 
of perfection. The ‘‘Missile Engineering 
Handbook”’ reflects the attention to detail 
which produces a relatively error-free 
book. In addition, most sources of data 
are listed so that the reader may verify the 
accuracy. 

The “Missile Engineering Handbook’ 
is the fourth volume in a series which 
provides a basis for formal instruction or 
self-study in guided missiles and space- 
flight. The other volumes of the series 
present the physical descriptions and 
mathematical derivations. This handbook 
is a collection of data useful for the 
preliminary analysis and initial design 
studies of future missiles. Much of the 
design engineer’s time is consumed search- 
ing through scientific journals, books and 
reports to find pertinent data; this hand- 
book should help to reduce time spent in 
search. 


Fundamentals of Advanced Missiles, by 
Richard B. Dow, John Wiley and Sons, 
Inc., New York, 1958, xvi + 567 pp. 

$11.75. 

Reviewed by JoHN GusTAVSON 

1, Grand Central Rocket Co. 


iw, 


This volume presents an example of 
the difficulties an author must face when 
writing a book on a vehicle system. He 
is squeezed between the danger of omit- 
ting important parts and the equally 
serious risk of spreading his treatment 
thin. 

The present volume discusses exterior 
ballistics, fluid mechanics, aerodynamics 
propulsion, statistical analysis, micro- 
waves, infrared radiation, radar, guid- 
ance and missile systems, in that order. 
Admittedly, most of these subjects 
form the fundamentals of advanced mis- 
siles, but this reviewer would have liked to 
have seen them treated in a different 


fashion. I would take a missile, break 


it apart, indicate the function of each 
subsystem, such as guidance, propulsion 
and control surfaces. Then treat each 
underlying science in a cursory manner, 
making certain that there is a solid red 
line throughout the volume, and making 
sure that no major subsystem is left out. 
In this book, I missed sections on struc- 
tures, warheads and rocket control of 
large missiles. 

The author has included references on 
each page in the form of footnotes, but 
unfortunately, they are not brought up to 
the date of publication. The book would 
have benefited by the inclusion of a 
literature list at the end of each chapter. 

The treatment in the many different 
sections of the book varies from good and 
basic to light and outdated. 

I suggest that the reader consider this 
book as a conglomerate of separate 
treatises on some of the topics related 
to our missile age. The book may serve 
as an introduction to these separate 
topics, but unfortunately, not to the sub- 
ject of advanced missiles. 


Introduction to the Theory of Compress- 
ible Flow, by Shih-I Pai, Van Nostrand 
Co., Princeton, N. J., 1959, xiii + 385 
pp. $9.75. 

Reviewed by James J. KAuZLARICH 
Worcester Polytechnic Institute 


The opening sections of the text will 
alert the reader to the fact that Dr. Pai 
intends to make as much use of modern 
mathematical notation as possible. For 
instance, the topics of fundamental 
equations, material differentiation, and 
boundary and initial conditions make up 
most of Chapter I. Chapter II is made 
up of a few pages on thermodynamics, and 
Chapter III opens the essential subject 
matter of the text, one-dimensional flow of 
an inviscid compressible fluid. The 
following chapter, on shock waves, 
completes the part of the text that deals 
only with one-dimensional motion. 

By Chapter IV another characteristic of 
this book becomes evident; the author 
presents the material at a level that 
presupposes a firm grasp of physical 
concepts. Here along with additional 
concepts and definitions, stagnation en- 
thalpy is introduced as a definition, 
considerations of constant cross-sectional 
area nozzles are taken up, and the idea 
that “many practical flow problems may 
be approximated by isentropic flow’’ is 
mentioned. It appears that the title of the 
book would be more appropriate if the 
word introduction were deleted. 

The text is based on the author’s 
lecture notes, and is intended to give 
essential results and useful techniques in 
the theoretical analysis of compressible 
flow. In carrying this out, the text 
covers the subject from the point of view 
of applied mathematics. 

At Chapter V, fundamental equations 
of the aerodynamics of a compressible 
inviscid and non-heat-conducting fluid, 
the author gives what can be considered 
the usual coverage in an introductory 
gasdynamics course, and in the following 


310 pages investigates many and varied 
topics in compressible fluid mechanics. 
Chapter V sets the stage for many 
chapters on methods of solution of the 
fundamental equations, i.e., VI on method 
of small perturbation, VII two- 
dimensional subsonic steady potential 
flow, VIII on hodograph method, IX on 
exact solutions of two-dimensional isen- 
tropic steady flow equations, X on two- 
dimensional steady transonic and hyper- 
sonic flows, XI on method of character- 
istics and XII on linearized theory of 
three-dimensional potential flow. 

In Chapter XIII on nonlinear theory of 
three-dimensional compressible flow, topics 
of considerable interest are treated. The 
problem of supersonic flow over a cone it 
zero angle of attack is taken up along with 
a study of transonic, hypersonic and supe r- 
sonic similarity laws as well as nonstea: 
compressible flow. Chapter XIV is 
concerned with anisentropic (rotation: 
flow of inviscid compressible fluid, and 
especially, flow conditions behind a 
curved shock. 

In Chapter XV, the author extends his 
treatment of compressible fluid flow to 
include viscous flow effects. This chapter 
is mainly concerned with the development 
of the Navier-Stokes equation and the 
energy equation. With these develo}- 
ments completed, it is now possible to take 
up a study of shock wave structure. 

Chapter XVI on boundary layer flow, 
and XVII on magnetogasdynamics, con- 
stitute the remainder of the text. In both 
the author has presented the subject 
material in concise form. Actually only 
enough space is taken to present a few of 
the more important topics from the 
extensive literature on these subjects. 

This book covers a large range of 
material in as digestible a form as may be 
possible considering the outline assumed. 
The text omits much of the physical 
description of typical problems, but this is 
not a handicap if one is familiar with the 
applications. For further study, the 
most important references are listed at 
the end of each chapter. 


Astronomy, by Theodore G. Mehlin, 
John Wiley and Sons, Inc., New York, 
392 pp. $7.95. 

Reviewed by 8S. F. SIncER 
University of Maryland 


This volume is a very elementary treat- 
ment of astronomy written in a simple 
and understandable way. It is well 
illustrated, uses a minimum of mathe- 
matics and also includes fairly simple 
problems. It is to be viewed more as a 
cultural, historical and phenomenological 
introduction to astronomy, and will be 
appreciated by those who wish to become 
acquainted with the subject because of its 
cultural implications and to learn the 
significance of space research. It will 
be a disappointment to anyone who wishes 
to explore the subject in a more technica! 
manner and become informed on astro- 
physics or celestial mechanics. The book 
is obviously intended for a course for 
liberal arts students and can be recom- 


mended for this purpose. 
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from PLENUM PRESS... 


ADVANCES IN THE ASTRONAUTICAL SCIENCES, Vol. 5 


This latest volume in the series is the proceedings of the Western 
National Meeting of the American Astronautical Society which was 
ADVANCES in held in August, 1959. The volume includes the opening address by Dr. 
Wernher von Braun, and 26 reports which represent original technical 
research both from theoretical and experimental aspects. Topics under 
discussion include Space Mechanics, Control and Guidance; Advanced 
Propulsion and Power; Astronautical Systems and Space Vehicle De- 
sign; Space Communication and Instrumentation; Lunar and Planetary 
Environment. This detailed picture of the scope of current astronautical 
SCIENGES vison 5 research will prove invaluable to the scientist, as well as form an excel- 
lent guide for the informed layman who desires an insight into the 
current problems of astronautics. 

cloth profusely illustrated $8.00 

Also available: Vols. 1-4, $8.00 each. (Price for complete 4-volume set: $28.00) 
Detailed tables of contents upon request. 


ASTRONAUTICAL 


EXPLODING WIRES 
Edited by William G. Chace and Howard K. Moore | 


The increasing contemporary importance of this scientific phenomenon 
was dramatically indicated by the unexpectedly large attendance at the 
1959 Conference on Exploding Wires. This volume, a complete record 
of the conference, includes reports of theoretical research in the field 
and new developments and applications. In this book—the first ever pub- 
lished on the subject in any language-—scientists, students, and tech- 
nologists will find a wealth of incisive material dealing with high-speed 
photographic techniques, shock waves, and the physics of fluids. 


cloth 374 pages _ profusely illustrated $9.50 


ARTIFICIAL EARTH SATELLITES Volumes 1 and 2 ; 
Edited by L. V. Kurnosova Translated from Russian an 


The publication of this double-volume book marks the appearance of 
the first translation of the revised texts of the reports read at the Fifth 
Assembly of the Special IGY Committee, which took place at Moscow 
in August, 1958. Volume 1 explores data obtained by the first two 
satellites and includes discussions on the ionosphere, cosmic radiation, 
and space biology. (Data collected on the first astronaut—the dog Laika 
—are presented.) Volume 2 documents work carried out with the third 
satellite (which contained a complete space laboratory), as well as 


with rockets. 


cloth 236 pages illustrated $9.50 


ADVANCES IN CRYOGENIC ENGINEERING, Volumes 1-5 
Edited by K. D. Timmerhaus 


These five volumes are the proceedings of the Cryogenic Engineering 
Conferences held in 1954, 1956-1959. (There was no meeting in 1955.) 
Subjects under discussion include: Cryogenic Tankage for Space Flight 
Applications; A Study of Condensing-Vacuum Insulation; Multiple- 
Layer Insulation; Nucleate- and Film-Boiling Studies with Liquid Hydro- 
gen; The Liquefaction of Natural Gas; Some Mechanical Properties of 
Magnesium Alloys at Low Temperatures; An Experimental Study Con- 
cerning the Pressurization and Stratification of Liquid Hydrogen; 
Maximum-Rate Theory of Impact Sensitivity. 


cloth over 250 pages in each volume $13.50 per volume aed ; 
(Price for complete set of 4 back volumes: $35.00) a Ignition of cold hydrogen gas by Very pistol 
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Propulsion and Power 


(Combustion Systems) 


Spacecraft Power from Open-cycle 
Chemical Turbines, by Antonio Orsini, 
ASTRONAUTICS, vol. 4, Dec. 1959, p. 36. 


Water Rockets, by Denis J. Zigrang, 
ASTRONAUTICS, vol. 4, Dec. 1959, p. 42. 

Solid-fuel Ramjet, by H. Powell Jen- 
kins Jr., ASTRONAUTICS, vol. 4, Dec. 1959, 
p. 44. 


Control of Combustion-chamber Pres- 
sure and Oxidant-fuel Ratio for a Re- 
generatively Cooled Hydrogen-fluorine 
Rocket Engine, by Edward W. Otto and 
Richard A. Flage, NASA TN D-82, Nov. 
1959, 61 pp. 

Trends in the Development of Liquid 
Propellant Rocket Engines, by E. G. D. 
Andrews and A. W. T. Mottram, Aeron. 
Quart., vol. 10, Aug. 1959, pp. 199-210. 


Recent Advances in Ramjet Combus- 
tion, by Gordon L. Dugger, ARS Jour- 
NAL, vol. 29, Noy. 1959, pp. 819-827. 


Vaporization of Propellants in Rocket 
Engines, by Richard J. Priem and Marcus 
F. Heidmann, ARS JourNat, vol. 29, 
Nov. 1959, pp. 828-835. 

Unit Converts Waste Rocket Heat to 
Power, by Michael Yaffee, Aviation Week, 
vol. 71, Nov. 23, 1959, pp. 92-93, 97, 99, 
100. 

A Detonation Wave Hypersonic Ramjet, 
by W. H. Sargent and R. A. Gross, Fair- 
child Engine & Airplane Corp., Fairchild 
Engine Div., June 1959, 40 pp. 

A Three-stage Solid-fuel Sounding 
Rocket System Suitable for Research at 
Altitudes Near 200 Nautical Miles, by 
Waldo L. Dickens and Earl C. Hastings 
Jr., NASA TN D-219, Dec. 1959, 25 pp. 

Analysis of the Transient Radiation 
Heat Transfer of an Uncooled Rocket En- 
gine Operating Outside Earth’s Atmos- 
phere, by William H. Robbins, NASA 
TN D-62,Dec. 1959, 25 pp. 

Development of TE-60 Plastic Parts As- 
sembly, T'hiokol Chem. Corp., Elkton Div., 
Final Rep., Aug. 29, 1955-March 31, 
1959, Rep. £120-59, Dec. 1959, 60 pp. 

Performance Investigation of a Non- 
pumping Rocket-ejector System for Alti- 
tude Simulation, by Anthony Fortini, 
NASA TN D-257, Dec. 1959, 33 pp. 

The Hydrogen-peroxide Rocket Reac- 
tion-control System for the X-1B Research 
Airplane, by James E. Love and Wendell 
H. Stillwell, NASA TN D-185, Dec. 1959 
28 pp. 

Rocket Morphology and Nomenclature, 
by Henry Burlage Jr., ASTRONAUTICS, vol. 
5, Jan. 1960, p. 42-43, 75. 


Radial Equilibrium in Supersonic Com- 
pressors, by A. G. Hammitt and S. M. 


Eprtor’s Note: Contributions from Pro- 
fessors E. R. G. Eckert, J. P. Hartnett, T. 
F. Irvine Jr. and P. J. Schneider of ‘the 
Heat Transfer Laboratory, University of 


Minnesota, are me acknowledged. 
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Bogdonoff, J. Basic Engng., vol. 81D, Dec. 
1959 (Trans., ASME), pp. 549-558. 


Resumé of the Supersonic-compressor 
Research at NACA Lewis Laboratory, 
by W. W. Wilcox, E. R. Rysl and M. J. 
Hartmann, J. Basic Engng., vol. 81D, 
Dec. 1959 (Trans., ASME), pp. 559-568. 


Chemical Propulsion: State of the Art, 
by Kurt R. Stehling, Space/Aeron., vol. 32, 
no. 6, Dec. 1959, pp. 42-46. 

The Rocket: A Past and Future History 
by Norman P. Gentieu and Meyer M. 
Markiwotz, Indl. Res., vol. 1, Nov.-Dec. 
1959, pp. 70-78. 

Investigation of Thrust and Drag Char- 
acteristics of a Plug-type Exhaust Nozzle, 
by Donald P. Hearth and Gerald C. Gor- 
ton, NACA RM E53L16, Feb. 1954, 39 
pp. (Declassified from Confidential by 
authority of NASA Pub. Ann. 4, p. 19, 
3/16/59.) 

Comparison of Two Methods of Modu- 
lating the Throat Area of Convergent 
Plug Nozzles, by H. George Krull and 
William T. Beale, NACA RM E54L08, 
May 1955, 20 pp. (Declassified from 
Confidential by authority of NASA Pub. 
Ann. 4, p. 20, 3/16/59.) 


Thrust Characteristics of a Series ot 
Convergent-divergent Exhaust Nozzles af 
Subsonic and Supersonic Flight Speeds, 
by Evan A. Fradenburgh, Gerald C. 
Gorton and Andrew Beke, NACA RM 
E53L23, March 1954, 24 pp. (Declassi- 
fied from Confidential by authority of 
NASA Pub. Ann. 4, p. 19, 3/16/59.) 

Effects of Several Geometric Variables 
on Internal Performance of Short Con- 
vergent-divergent Exhaust Nozzles, by 
Fred W. Steffen, H. George Krull and 
Ralph F. Schmiedlin, NACA RM E54- 
L09, Feb. 1955, 24 pp. (Declassified from 
Confidential by authority of NASA Pub. 
Ann. 4, p. 20, 3/16/59.) 


Thrust and Drag Characteristics of a 
Convergent-divergent Nozzle with Vari- 
ous Exhaust Jet Temperatures, by Don- 
ald P. Hearth and Fred A. Wilcox, NACA 
RM E53L23b, March 1954, 20 pp. (De- 
classified from Confidential by authority of 
NASA Pub. Ann. 4, p. 19, 3/16/59.) 


Propulsion and Power 
(Non-Combustion) 
Hearing . . .on the Aircraft Nuclear Pro- 


pulsion Program, U.S. Congress, 86th, 1st 
Session, JointCommittee on Atomic Energy, 
July 23, 1959, 418 pp. 

Design Considerations of Steady DC 
Magnetohydrodynamic Electrical Power 
Generator, by George W. Sutton, General 
Electric Co., Missile & Space Vehicle Dept.. 
TIS R59SD432, Sept. 1959, 38 pp. 


Dynamic Analysis of a Nuclear Rocket 
Engine System, by Bernard R. Felix and 
Richard J. Bohl, ARS Journat, vol. 29, 
Nov. 1959, pp. 853-862. 


New Ion Source for Propulsion of Space 


Vehicles, by C. D. Mosk, H. E. Banta, 


M. H. Smith, Associate Editor 


2 The James Forrestal Research Center, Princeton University 


J. W. Johnson and R. F. King, ARS Jour- 
NAL, vol. 29, Nov. 1959, pp. 868-869. 

Nuclear Auxiliary Space System May Be 
Operational Within ag by Craig Lewis, 
Aviation Week, vol. 71, Nov. 23, 1959, p. 
35. 

Thermoelectron Engines: Future 
Power Sources? by G. N. Hatsopoulos, 
J. Welsh and E. Langberg, Electronics, vol. 
32, Nov. 13, 1959, pp. 69-72. 

Thermomagnetic Generator, by J. F. 
Elliott, J. Appl. Phys., vol. 30, Nov. 
1959, pp. 1774-1777. 

Ionic Propulsion, by Robert H. Fox, 
Rev. Astron. Sci., vol. 1, no. 3, July—Sept. 
1959, pp. 7-12. 

Competition Rises in Ion, Plasma 
Fields, by Philip J. Klass, Aviation Week, 
vol. 71, Dec. 7, 1959, pp. 83, 87, 89-90. 

New Ideas for Nuclear Reactors, by Eu- 
gene P. Wigner, Nuclear Sci. & Engng., vol. 
6, no. 5, Nov. 1959, pp. 420-440. 

Controlling an Airborne Nuclear Reac- 
tor, by Harry L. Loates, Space/ Aeron. Res. 
& Dev. Handbook, 1959-1960, pp. F-6-9. 

Advanced Airborne APUs (table), 
Space/ Aeron. Res. & Dev. Handbook, 1959- 
1960, pp. F-12-15. 

Power from Radioisotopes (nomo- 
graph), Space/ Aeron. Res. & Dev. Hand- 
book, 1959-1960, pp. G-16-17. 

Finding the Optimum Load for Solar 
Cells, by A. J. Heager and T. R. Nisbet, 
Space/Aeron., vol. 32, no. 6, Dec. 1959, pp. 
153-154, 158, 160, 166, 168, 172, 174. 

Magnetohydrodynamic Acceleration of 
Slightly Ionized Viscously Contained 
Gases, by G. Janes and James A. Fay, 
Aveo Corp., Avco-Everett Res. Lab., Oct. 
1959, 18 pp. (2nd Symposium on Ad- 
vanced Propulsion Concepts, Oct. 1959.) 

Pulsed Plasma Accelerator, by Thomas 
L. Thourson, Borg-Warner Corp., 1959, 21 
pp. (2nd Symposium on Advanced Pro- 
pulsion Concepts, Oct. 1959.) 

Flight Magnetohydrodynamics, by 
Arthur Kantrowitz, Avco-Everett Res. Lab., 
Res. Rep. 51, March 1959, 18 pp. 

A Comparison of the Specific Thrust of 
Ion and Plasma Drive Accelerators, by 8. 
W. Hash, Lockheed Aircraft Corp., Missile 
& Space Div., Oct. 1959, 15 pp. 

The Arc Jet, by Gabriel M. Giannini, 
Plasmadyne Corp., M-11, 1959, 23 pp. 

A Critical Evaluation of the Ion Rocket, 
by T. M. Littman, North American Avia- 
tion, Inc., Rocketdyne, 29 pp. (2nd Sym- 
posium on Advanced Propulsion Concepts, 
Oct. 1959.) 

The Colloid Rocket: Progress Toward a 
Charged-liquid-colloid Propulsion Sys- 
tem, by Robert D. Schultz and Lane K. 
Branson, Aerojet-General Corp., Space Tech. 
Div. (AEC NYO-7309), 1959, 20 pp. 

Aircraft Nuclear Propulsion Program, 
(U. S. Congress, 86th, 1st Session, Jeint 
Comm. on Atomic Energy), Sept. 1959, 25 
pp. 

Space Propulsion (U. S. Congress, 
86th, 1st Session, H.R., Committee on Sci. & 
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Astron., House Rep. no. 1087), Aug. 1959, 
15 pp. 

Space Propulsion (U’. S. Congress, 86th, 
ist Session) Hearings before the Comm. on 
Sci. & Astron., March 16-20, 23, 1959, 
no. 16, 1959, 307 pp. 

Annual Review of Nuclear Science, vol. 
9, Palo Alto, Calif., Annual Rev., Inc., 
1959, 625 pp. 

Nuclear Fission, by I. Halpern, pp. 

245-342, 197 ref. 

High-temperature Plasma Research 

and Controlled Fusion, by R. F. Post, 

pp. 367-436, 66 ref. 

Manned Nuclear Space Systems, Part 
I: High-thrust Nuclear Systems, by E. B. 
Konecci, R. F. Trapp and M. M. Hunter, 
Aero/Space Engng., vol. 19, Jan. 1960, pp. 
34-41. 

Neon and Oxygen Thermonuclear Re- 
actions, by A. G. W. Cameron, Astrophys. 
J., vol. 130, Nov. 1959, pp. 895-915. 


Propellants and Combustion 


Viscosity of Fluids at High Pressures; a 
Rotating Cylinder Viscometer and the 
Viscosity of N-pentane, by H. H. Reamer, 
G. Cokelet and B. H. Sage, Project Squid 
Tech. Rep. CIT-4-P (ASTIA A D 208670), 
Jan. 1959, 33 pp. (Available only on 
microcards. ) 

Preliminary Study of a Piston Pump 
for Cryogenic Fluids, by Arnold E. 
Biermann and Robert C. Kohl, NASA 
Memo 3-6-59E, March 1959, 27 pp. 

Effects of Gas Motion on Heterogeneous 
Combustion; Natural Convection, Steady 
Forced Convection, Standing Acoustic 
Waves, and Shock Waves, Final Summary 
Report, by William T. Webber, Wright Air 
Dev. Center, Tech. Rep. no. 59-50, April 
1959, 88 pp. 

Kinetic Studies of Nitric Oxide, by Kurt 
L. Wray and J. Derek Teare, Avco- 
Everett Res. Lab., Res. Note no. 134, June 
1959, 13 pp. 

Symposium on Thermal Properties, 
Purdue University, 1959. Thermody- 
namic and Transport Properties of Gases, 
Liquids and Solids, Papers Presented, 
New York, Trans., ASME; McGraw-Hill 
Book Co., 1959, 472 pp. 

Thermochemistry and Thermodynamic 

Functions of Some Boron Compounds, 

by W. H. Evans, E. J. Prosen and D. D. 

Wagman, pp. 226-235. 

The Combustion of Boron Hydrides, by 

W. G. Ber! and W. Renich, pp. 236-250. 

A Generalized Data-fitting Routine for 

the LGP-30 Computer: Application to 

Real-gas Properties of Diborane, by 

H. J. Galbraith and J. F. Masi, pp. 

251-274. 

The Low Temperature Thermal and 

Chemical Thermodynamic Properties of 

Boron Compounds, by E. F. Westrum, 

pp. 275-283. 

Thermal Conductivity of Liquid Ozone, 

by T. E. Waterman, pp. 301-303. 

The Internal Ballistics of a Supergun, 
by J. N. Kapur, Appl. Sci. Res., Sect. A, 
vol. 8, no. 6, 1959, pp. 393-402. 

Synthesis of Storable Rocket Fuel Dur- 
ing Space Missions, by Eugene L. Colich- 
man, ASTRONAUTICS, vol. 4, Dec. 1959, p. 
32. 

Dependence of Damage Effects upon 
Detonation Parameters of Organic High 
Explosives, by Donna Price, Chem. Revs., 
vol. 59, Oct. 1959, pp. 801-826. 

Consumption of Oxygen Molecules in 
Hydrocarbon Flames Chiefly by Reaction 
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with Hydrogen Atoms, by C. P. Feni- 
more and G. W. Jones, J. Phys. Chem., vol. 
63, Nov. 1959, pp. 1834-1838. 

Detonation Induction Distances in 
Combustible Gaseous Mixtures at At- 
mospheric and Elevated Initial Pressures, 
I Methane-Oxygen, II Carbon Monoxide- 
Oxygen, III Hydrogen-Oxygen, by Loren 
E. Bollinger and Rudolph Edse, Wright 
Air Dev. Center, Tech. Rep., 58-591, March 
1959, 88 pp. 

Adiabatic Flame Temperatures, Com- 
bustion Gas Compositions, and Expansion 
Ratios of Combustible Gas Mixtures at 
Various Pressures, by William A. Strauss 
and Rudolph Edse, Wright Air Dev. Cen- 
ter, TN 59-40, July 1959, 43 pp. 

Study of Additive Effects in Solid Pro- 
pellant Burning, by Herman Burwasser, 
Aero. Chem. Res. Labs., Inc., Second Quart. 
Progr. Rep. TM-17, Oct. 1959, 8 pp. 

Total Radiation from Burning Solid Pro- 
pellant Strands, by David W. Blair and 
Kimball P. Hall, Princeton Univ., Aeron. 
Engng. Dept., Semi-annual Progr. Rep. for 
the period 1 April 1959-30 Sept. 1959, Oct. 
1959, 8 pp. 

Joint Army-Navy-Air Force Ignitability 
Panel, Proceedings of the Third Meeting 
held at Applied Physics Laboratory, Silver 
Spring, Md., Solid Propellant Info. 
Agency, June 1959, 69 pp. 

Quarterly Progress Report on Engineer- 
ing Research (U), Rohm & Haas Co., Red- 
stone Arsenal Res. Div., Rep. P-59-17, 
June 15-Sept. 15, 1959, 17 pp. 

Spectrometric Determination of Chlo- 
rate Impurities in Ammonium Perchlorate, 
by Eugene A. Burns, Calif. Inst. Tech., Jet 
Prop. Lab. Progr. Rep. 30-12, May 1959, 8 
pp. 

Theoretical Performance of Ammonium 
Nitrite—Water Solutions, by John S. Gor- 
don, Wright Air Dev. Center, Tech. Rep. 
59-248, March 1959, 7 pp. 

Heterogeneous Combustion of Multi- 
component Fuels, by Bernard J. Wood, 
Henry Wise and 8. Henry Inami, NASA 
TN D-206, Nov. 1959, 32 pp. 

Production and Measurement of Single 
Drops, Sprays, and Solid Suspensions, by 
James A. Browning, Reprinted from Ad- 
vances in Chemistry, Series 20 (Amer. 
Chem. Soc.), Washington, D. C., 1958, pp. 
136-154. 

The Chemistry of Boron Hydrides and 
Related Hydrides, by R. M. Parry, R. C. 
Taylor, C. E. Nordman, G. Kodama and 
S. G. Shore, Wright Air Dev. Center, Tech. 
Rep. 59-207, May 1959, 178 pp., 141 refs. 

Theory of Particle Combustion at High 
Pressures, by D. B. Spalding, ARS 
JOURNAL, vol. 29, Nov. 1959, pp. 828-835. 

Solubility of Air, N. and O, in n-Propyl 
Nitrate, by Adolph B. Amster and Joseph 
B. Levy, ARS Journat, vol. 29, Nov. 
1959, pp. 870-871. 

Solids Make Bid for Space Booster 
Role, by Michael Yaffe, Aviation Week, 
vol. 71, Nov. 23, 1959, pp. 32-34. 

Rocketdyne Studies Pumps, Future Li- 
quid Propellants, by Russel Hawkes, Avia- 
tion Week, vol. 71, Nov. 23, 1959, pp. 57, 
59, 61, 65, 67, 69, 73. 

Thermal Decomposition of Nitrous 
Oxide, by B. G. Reuben and J. W. Lin- 
nett, Trans., Faraday Soc., vol. 55, Sept. 
1959, pp. 1543-1553. 

Condensed Systems at High Pressures, 
by R. H. Wentorf Jr., J. Phys. Chem., vol. 
63, Nov. 1959, pp. 1934-1939. 

Dielectric Constant of Liquid Ozone and 
Liquid Ozone-Oxygen Mixtures, by 
Charles K. Hersh, Gerald M. Platz and 
Raymond J. Swehla, J. Phys. Chem., vol. 


63, Nov. 1959, pp. 1968-1969. 

Chemical Reactions of Active Nitrogen, 
by Harold A. Dewhurst, J. Phys. Chem., 
vol. 63, Nov. 1959, pp. 1976-1977. 

Electric and Magnetic Properties of the 
Hydrogen Molecule, by T. P. Das and R. 
Bersohn, Phys. Rev., vol. 115, no. 4, Aug. 
15, 1959, pp. 897-910. 

On the Existence of the Negative Nitro- 
gen Ion, by Ya. M. Fogel’, V. F. Kozlov 
and A. A. Kalmykov, Soviet Phys.-J ETP, 
vol. 36(9), no. 5, Nov. 1959, pp. 963-964. 

Blowoff of Propane and Hydrogen 
Diffusion Flames at High Mach Number, 
Ramjet Conditions, by R. J. Bacigalupi 
and E. A. Lezberg, NASA TN D-67, Dee. 
1959, 24 pp. 

Investigation of the Products from Ex- 
plosions of Diborane-Air Mixtures, by 
Robert Brown, Johns Hopkins Univ., 
Appl. Phys. Lab., CF-2844, Sept. 1959, 
6 pp. 

Studies on Deflagration to Detonation 
in Propellants and Explosives, U. S. Bur. 
Mines, Explosives Res. Lab., Quart. 
Rep., July 1-Sept. 30, 1959, 10 pp. 

Combustion Instability in Liquid Pro- 
pellant Rocket Motors (29th Quarterly 
Progress for the period 1 May 1959 to 31 
July 1959), Princeton Univ. Aeron. Engng. 
Dept., Rep. 216-CC, Sept. 1959, 19 pp. 

A One-dimensional Theory of Liquid- 
fuel Rocket Combustion, by D. B. Spald- 
ing, Gt. Brit., Aeron. Res. Council, Current 
Paper 445, May 1958, 18 pp. 

Chemical Sampling Downstream of 
Lean, Flat Hydrogen and Propane Flames, 
by Burton D. Fine, NASA TN D-198, 
Dec. 1959, 25 pp. 

Investigation of Liquid Nitric Oxide as 
a Rocket Oxidizer, by Loren E. Bollinger 
and Rudolph Edse, Wright Air Dev. 
Center, Tech. Rep. 59-12, May 1959, 22 pp. 

The Preparation of 3,3,3-tris( Nitrato- 
methyl )-Propene-1 and 3,3-bis(Nitrate- 
methyl )-Butene-1, by Lester P. Kuhn and 
Alan C. Duckworth, Aberdeen Proving 
Ground, Ball. Res. Lab., Mem. Rep. 1230, 
Aug. 1959, 24 pp. 

Binary Nitrogen Compounds of the 
Elements: A Literature Survey, by 
Rosemary G. Ehl, Rodney J. Sime and 
John L. Margrave, Wright Air Dev. 
Center, TN 59-115, June 1959, 75 pp. 

Calculation of Shock Hugoniots and 
Related Quantities for Nitrogen and 
Oxygen, by R. A. Alpher and H. D. 
Greber, General Electric Co., Res. Lab., 
Schenectady, NY, Rep. 58-RL-1915, April 
1958, 25 pp. 

The Mechanism of Polytetrafluoro- 
ethylene Pyrolysis, by Henry L. Friedman, 
General Electric Co., Missile & Space 
Vehicle Dept., TIS R598D385, June 1959, 
26 pp. 

Kinetics of the Decomposition of N.O 
at High Temperature, by Wendell M. 
Graven, J. Amer. Chem. Soc., vol. 81, 
Dec. 5, 1959, pp. 6190-6192. 

Vapor Pressure of Diborane, by Leo J. 
Paridon and George E. MacWood, J. 
Phys. Chem., vol. 63, Dec. 1959, p. 1997. 

The Heat of Vaporization of Diborane, 
by Leo J. Paridon, George E. MacWood 
and Jih-Heng Hu, J. Phys. Chem., vol. 63, 
Dee. 1959, pp. 1998-1999. 

Flame Spectra: The Effect of Deutera- 
tion on the Hydrocarbon Flame Bands, 
by G. N. Spokes and A. G. Gaydon, Proc., 
Phys. Soc., London, vol. 74, Nov. 1, 1959, 
pp. 639-642. 

LCAO Wave Functions for Hydrogen 
Fluoride with Hartee-Fock Atomic Orbi- 
tals, by Arnold M. Karo and Leland C. 
Allen, J. Chem. Phys., vol. 31, Oct. 1959, 


pp. 968-977. 
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Metal Combustion Processes, by Irvin 
Glassman, Princeton Univ., Aeron. Engng. 
Dept., Rep. 473, Aug. 1959, 37 pp. 

A One-dimensional Theory of Liquid- 
fuel Rocket Combustion, by J. Adler, Gt. 
Brit., Aeron. Res. Council, Current Paper 
446, May 1958, 7 pp. 

Chemical Nonequilibrium Effects on 
Hydrogen Rocket Impulse at Low Pres- 
sures, by J. Gordon Hall, A. Q. Eschen- 
roeder and J. J. Klein, Cornell Aeron. Lab., 
Rep. AD-1118-A-8, Nov. 1959, 10 pp. 

Ignition of Solid Propellants by Com- 
pression of Gases Entrapped on an Inter- 
nal Cavity, by C. H. Parr, Rohm & Haas 
Co., Quart. Progr. Rep. on Engng. Res., 
Rep. P-59-11, March 15-—June 15, 1959, 
pp. 1-17. 

The Controlled Thermal Decomposition 
of Cellulose Nitrate, VII: Carbonyl 
Compounds, by M. L. Wolfrom and G, P. 
Arsenault, Ohio State Univ., Res. Found., 
Tech. Rep. 675-14, Oct. 1959, 21 pp. 


Mollier Diagrams for the Hydrogen- 
Oxygen System and Experimental Re- 
sults of Test on a Rocket Motor of 500- 
pound Thrust, by Dwight I. Baker, Calif. 
Inst. Techn., Jet Prop. Lab., Rep. 20-127, 
May 1959, 109 pp. 

Theory of the Burning of Mono-propel- 
lant Droplets, by D. B. Spalding and V. K. 
Jain, Gt. Brit., Aeron. Res. Council, Current 
Paper 447, May 1958, 9 pp., 4 figs. 

Analytical Solution of a Flame with 
Cylindrical Symmetry, by J. Menkes, 
Calif. Inst. Techn., Jet Prop. Lab., Progr. 
Rep. 30-17, Oct. 1959, 8 pp. 


Temperature Distribution in Sooty 
Flames, by F. Réssler, Annalen der Phystk, 
vol. 4, nos. 6-8, Oct. 13, 1959, pp. 396-422. 
(In German.) 


A Study of the Interior Ballistic Equa- 
tions, by R. A. Struble, Arch. for Rational 
Mech. & Anal., vol. 3, no. 5, Sept. 30, 1959, 
pp. 397-421. 

Theory of Particle Combustion at High 
Pressures, by D. B. Spalding, ARS 
JOURNAL, vol. 29, Nov. 1959, pp. 828-835. 


Vaporization of Propellants in Rocket 
Engines, by Richard J. Priem and Marcus 
F. Heidmann, ARS JournaL, vol. 29, Nov. 
1959, pp. 836-842. 

Solubility of Air, N. and O, in n-Propyl 
Nitrate, by Adolph B. Amster and Joseph 
B. Levy, ARS Journat, vol. 29. Nov. 
1959, pp. 870-871. 


Esso Accelerates Solid Propellant 
Search, by Michael Yaffee, Aviation Week, 
vol. 71, Dec. 14, 1959, pp. 68-69, 71, 75, 
44-48. 

Effect of Argon and Helium on the 
Thermal Conductivity of the N.0,—@2NO, 
System, by Kenneth P. Coffin, J. Chem. 
Phys., vol. 31, Nov. 1959, pp. 1290-1297. 

Participation of Vibration in Exchange 
Reactions, by Felix T. Smith, J. Chem. 
Phys., vol. 31, Nov. 1959, pp. 1852-1359. 

A Calculation of the Rates of the Ortho- 
para Conversions and Isotope Exchanges 
in Hydrogen, by Isaiah Shavitt, J. Chem. 
Phys., vol. 31, Nov. 1959, pp. 1359-1367. 


Comparison of Hot and Thermal Reac- 
tions, by J. L. Magee and W. H. Hamill, 
J. Chem. Phys., vol. 31, Nov. 1959, pp. 
1380-1386. 

Mean-first-passage Times and the Col- 
lision Theory of Bimolecular Reactions, 
by B. Widom, J. Chem. Phys., vol. 31 
Nov. 1959, pp. 1387-1394. 

Reaction of Acetylene with Methyl 
Radicals, by Charles M. Drew and Alvin 
S. Gordon, J. Chem. Phys., vol. 31, Nov. 
1959, pp. 1416-1417. 
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NMR Spectra of Normal and Deuter- 
ated Methylacetylenes, by Wyman R. 
Vaughan and Robert C. Taylor, J. Chem. 
Phys., vol. 31, Nov. 1959, pp. 1425-1426. 

Basic Principles of Combustion-model 
Research, by oe A. Putnam and E. W. 
Ungar, J. Engng. for Power (Trans., 
ASME), vol. A 81, Oct. 1959, pp. 383-388" 


The Effect of Fuel Types and Admission 
Method upon Combustion Efficiency, by 
H. N. McManus Jr., W. E. Ibele and T. E. 
Murphy, J. Engng. for Power (Trans., 
ASME), vol. A 81, Oct. 1959, pp. 423- 
427. 

Pressure and Velocity Measurements 
on Detonation in 
Oxygen Mixtures, by D. H. Edwards, 

Gt Williams and J. C. Breeze, J. Fluid. 
Grech. vol. 6, Nov. 1959, pp. 497-517. 

Temple Studies High Temperature 
Flames, by William L. Doyle, Missiles & 
Rockets, vol. 5, Dec. 14, 1959, pp. 17-18. 

Fluorine Compounds Point to Many 
New Propellant Possibilities, by John F. 
Gall, Space/Aeron. Res. & Dev. Handbook, 
1959-1960, pp. C-5-9. 

Handling and Storage of High Energy 
Oxidizers, by J. M. Siegmund and H. R. 
Neumark, Space/Aeron. Res. & Dev. 
Handbook, 1959-1960, pp. L-9-14. 


Materials and Structures 


Some Contributions to the Heat Con- 
duction and Thermal Stresses Analysis in 
Aircraft and Missile Structures, by L. 
Broglio, Aerotechnica, vol. 39, no. 2, April 
1959, pp. 53-65. (In Italian.) 

Materials for Rockets and Reactors, by 
J. A. Ternisien, Fusées et Recherche Aéron., 
no. 13, Aug. 1959, pp. 74-84. (In French.) 


Fatigue Damage During Complex Stress 
Histories, by H. W. Liu and H. T. Corten, 
NASA TN D-256, Nov. 1959, 56 pp. 

Tensile Properties of 6A-14V Titanium- 
alloy Sheet under Rapid-heating and 
Constant-temperature Conditions, by 
Howard L. Price. NASA TN D-121, Nov. 
1959, 26 pp. 

High Temperature Behavior of Teflon, 
by Tunis Wentink Jr., Avco-Everett Res. 
Lab., Res. Rep. no. 55, July 1959, 13 pp. 

Ceramic Fabrication Processes, W. D. 
Kingery, Ed., John Wiley & Sons, N. Y., 
1958, 235 pp. 

Effect of Microstructure on the Electri- 

cal and Magnetic Properties of Ce- 

ramics, by George Economos, pp. 201- 

213. 

Effect of Microstructure on _ the 

Mechanical Properties of Ceramic Ma- 

terials, by R. L. Coble, pp. 213-228, 46 

refs. 

Beryllium Applications in Structural 
Design, by L. A. Riedinger and R. E. 
Foster, Rev. Astron. Sci., vol. 1, no. 3, 
July-Sept. 1959, pp. 13-18. 

Stabilization Effects of Internal Pres- 
sure on the Buckling Behavior of Thin-wall 
Rocket Motor Cases, by William D. 
Jordan, Rohm & Haas Co., Quart. Progr. 
Rep. on Engng. Res., Rep. P-51-11, March 
15-June 15, 1959, pp. 18-36. 

Paint Gives Missiles Thermal Protec- 
tion, by Michael Yaffee, Aviation Week, 
vol. 71, Dec. 7, 1959, pp. 102-103, 105. 

Materials for Environmental Extremes, 
by G. Sideris, Electronics, vol. 32, Dec. 4, 
1959, pp. 81-96. 

Thermal Degradation of Polymers at 
High Temperatures, by Samuel! L. Mador- 
sky and Sidney Straus, J. Res., Natl. Bur. 
Standards, A: Phys. & Chem., vol. 63A, 
no. 3, Nov.—Dec., 1959, pp. 261-264. 


New Shielding Materials for ae 
temperature Application, by 
Hungerford, R. F. Mantey and L. P. Vax 
Maele, Nuclear Sci. & Engng., vol. 6, no. 
5, Nov. 1959, pp. 396-408. 

Titanium Provides Notch-ductile Rocket 
Cases, by A. N. Eshman, Space/Aeron., 
vol. 32, no. 6, Dec. 1959, pp. 73-74, 76, 78. 

Shock and Vibration in Electronic 
Mounting (nomographs), by H. C. Law- 
rence, Space/ Aeron. Res. and Dev. Hand- 
book, 1959-1960, pp. D-13-14, 15-16. 

Extreme Temperature Lubricants and 
Hydraulic Fluids, Space/Aeron. Res. & 
Dev. Handbook, 1959-1960, pp. G-9-12. 

Improved Plastics Meet Needs of High 
Performance Vehicles, by I. J. Gruntfest, 
Space/Aeron. Res. & Dev. Handbook, 
1959-1960, pp. H-6-9. 

Mechanical Properties of High Perform- 
ance Materials (table), by Alan V. Levy, 
Space/ Aeron. Res. & Dev. Handbook, 1959- 
1960, pp. H-10-11. 

Properties of Honeycomb Core (table), 
by D. C. Wilson, Space/Aeron. Res. & 
Dev. Handbook, 1959-1960, pp. H-12-13. 

Properties of Refractory Hard Metals 
and Ceramic Raw Materials (table), by 
Harry B. Porter, Space/Aeron. Res. & 
Dev. Handbook, 1959-1960, p. H-13. 


Stress Rupture Properties of Nickel and 
Cobalt Alloys (table), by Michel Urbain, 
Space/Aeron. Res. & Dev. Handbook, 
1959-1960, p. H-14. 


Fluid Dynamics, Heat Transfer 
and MHD 


An Investigation of Steady and Un- 
steady Heat Transfer in Rocket Motors, 


by Andre G. Condomines (Thesis, M.S.), 
Princeton Univ., Aeron. Engng. Dept., 
1959, 36 pp., 15 figs. 


Heat Transfer and Fluid Flow: A 
Bibliography of Selected Report Litera- 
ture, by U. S. Atomic Energy Comm., 
Tech. Info. Serv., TID-3305 (Supplement 
1), June 1958, 432 pp. 

Comparison between Molecular and 
Turbulent Diffusion Processes, by F. N. 
Frenkiel, David Taylor Model Basin, Rep. 
no. 1320, July 1959, pp. 204-216. (Re- 
printed from Int. Conf. Transport Prop. 
Statistical Mech., Brussels, 1956.) 

Survey of Thermal, Radiation, and 
Viscous Damping of Pulsating Air Bubbles 
in Water, by Charles Devin Jr., David 
Taylor Model Basin, Rep. 1329, Aug. 1959, 
37 pp. 

Flexure of a Two-dimensional Arc 
under Forced Convection, by Paul G. 
Thiene, Plasmadyne Corp., Rep. no. T2 
TN-089-334, Aug. 23, 1959, 80 pp. 

An Arc Tunnel for Magnetohydro- 
dynamic Studies, by Richard J. Rosa, 
Avco-Everett Res. Lab., Res. Note no. 132, 
June 1959, 11 pp. 

End Effects in Magnetohydrodynamic 
Channel Flow, by Frank Fishman, Avco- 
Everett Res. Lab., Res. Note no. 135, June 
1959, 14 pp. 

Transport Property Equations for Par- 
tially Ionized Gases, by C. W. Baulknight, 
General Electric, Missile & Space Vehicle 
Dept., TIS R598D428, Sept. 1959, 12 pp. 

Shock Front Thickness and Bulk Vis- 
cosity in Polyatomic Gases, by W. H. 
Andresen and D. F. Hornig, Brown Univ., 
Metcalf Res. Lab., Tech. R2p. 6, July 1955, 
13 pp. 

Use of a Theoretical Flow Model to 
Correlate Data for Film Cooling or Heat- 
ing an Adiabatic Wall by Tangential 
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Injection of Gases of Different Fluid 
Properties, by James E. Hatch and S. 
Stephen Papell, NASA TN D-130, Nov. 
1959, 43 pp. 

On Ablating Heat Shields for Satellite 
Recovery, by Steven Georgiev, Henry 
Hidalgo and Mac C. Adams, Avco-Everett 
Res. Lab., Res. Rep. 65, July 1959, 28 pp. 

A Theory of Ablation of Glassy Mate- 
rials for Laminar and Turbulent Heating, 
by Henry Hidalgo, Avco-Everett Res. Lab., 
Res. Rep. 62, June 1959, 27 pp. 

Heat and Mass Transfer Bibliography 
with Selected Abstracts, by W. L. Chris- 
tian and T. H. Schiffman, Wright Air Dev. 
Center, Tech. Rep. 58-186, May 1958, 167 
pp. 

The Effects of Electrode Contamination 
on the Properties of Air-arc Plasmas, by 
Harold Y. Wethaen, Milton Linevsky and 
John H. McGinn, General Electric Co., 
Missile & MSpace Vehicle Dept., TIS 
R598D427, Sept. 1959, 16 pp. 

The  Tailored-interface 
Shock Tunnel, by Charles E. Wittliff, 
Merle R. Wilson and Abraham Hertzberg, 
Cornell Aeron. Lab., Inc., March 12-20, 
1958, 26 pp., 18 diagrs. 

Electrical and Pressure Losses in a 
Magnetohydrodynamic Channel Due to 
End Current Loops, by George W. Sutton, 
General Electric Co.. Missile and Space 
Vehicle Dept., TIS R598D431, July 1959, 
37 pp. + diagrs. 

Effects of Air Dissociation and Ioniza- 
tion at Hypersonic Speed, by Gianni 
Jarre, Torino, Lab. Mec. Appl. del Politec., 
TN 14, July 1959, 29 pp. 

Non-equilibrium Theory of an Ideal- 
dissociating Gas Through a _ Conical 
Nozzle, by N. C. Freeman, Gt. Brit., 
Aeron. Res. Council, Current Paper 438, 
Aug. 1958, 10 pp. 

Mach Waves in Shock Wave Systems 
from a Detonating Solid Explosive, by 
D. N. Woodhead and R. Wilson, Gt. Brit. 
Min. Power, Safety in Mines Res. Estab., 
Res. Rep. 155, May 1959, 48 pp. 

Stability of a Horizontal Fluid Layer 
with Unsteady Heating from Below and 
Time-dependent Body Force, by Arthur 
W. Goldstein, NASA Tech. Rep. R-4, 
1959, 15 pp. 

Variational Theory of Chemical Reac- 
tion Rates Applied to Three-body Recom- 
bination, by James C. Keck, Avco-Everett 
Res. Lab., Res. Rep. 66, Sept. 1959, 62 pp. 

Compressible Flow Tables for Gases 
with Specific Heat Ratios from 1.10 to 
1.28, by K. E. Tempelmeyer and G. H. 
Sheraden, Arnold Engng. Dev. Center, TN 
58-9, March 1958, 173 pp. 

Spreading of Rocket Exhaust Jets at 
High Altitudes, by E. K. Latvala, Arnold 
Engng. Dev. Center, Tech. Rep. 59-11. 
June 1959, 37 pp. 

Experimental Studies of Shock Waves 
in Nitrogen, by John C. Camm and James 
C. Keck. Avco-Everett Res. Lab., Rep. 67, 
June 1959, 18 pp. 

Similarity of Far Noise Fields of Jets, 
by Walton L. Howes, NASA Tech. Rep. 
R-52, 1959, 90 pp., 26 figs. 

On Axisymmetrical Jets, with Applica- 
tion to Rocket Jet Flow Fields at High 
Altitudes, by W. T. Lord, Gt. Brit. Royal 
Aircr. Estab., Rep. Aero 2626. July 1959, 
37 pp. 

Laminar Mixing of Streams of Different 
Gases, by Luigi G. Napolitano and A. 
Pozzi. Univ. Naples, Inst. Aeron., TN 
1959, 60 pp. 

Magnetogasdynamic Flow Regimes, by 
Ching-Shi Liu and Ali Bulent Cambel, 
ARS Journat, vol. 29, Nov. 1959, pp. 
871-873. 


Marcu 1960 


Streamlined Conical Bodies under the 
Movements of Gases at Hypersonic 
Speeds, by A. L. Gonor, Akademiia Nauk 
SSSR, Izvestiia, Otdelenie Tekhnicheskikh 
Nauk, Mekhanika i Mashinostroenie, 1, 
1959, pp. 34-40. (In Russian.) 

Design of Axisymmetric Jet Nozzles of 
Minimum Weight, by L. E. Sternin, 
Akademiia Nauk SSSR, Izvestiia, Otdelenie 
Tekhnicheskikh Nauk, Mekhanika i 
Mashinostroenie, 1, 1959, pp. 41-45. (In 
Russian.) 

Motion of a Thin Solid Body under the 
Action of a Strong Shock Wave, by S. S. 
Grigorian, Akademiia Nauk SSSR, 
Izvestiia, Otdelenie Tekhnicheskikh Nauk, 
Mekhanika i Mashinostroenie, 1, 1959, pp. 
165-166. (In Russian.) 

Useful Interference of Wings and 
Fuselage at Hypersonic Speeds, by G. L. 
Grodzovskii, Akademiia Nauk SSSR, 
Tzvestiia, Otdelenie Tekhnicheskikh Nauk, 
Mekhanika i Mashinostroenie, 1, 1959, pp. 
170-173. (In Russian.) 

Hypersonic Flow of an Imperfect Gas 
Around a Body, by G. A. Liubimov, 
Akademiia Nauk SSSR, Izvestiia, Otdelenie 
Tekhnicheskikh Nauk, Mekhanika i 
Mashinostroenie, 1, 1959, pp. 173-176. 
(In Russian.) 

Electromotive Force in a Highly Ionized 
Plasma Moving Across a Magnetic Field, 
by M. Sakuntala, B. E. Clotfelter, W. B. 
Edwards and R. G. Fowler, J. Appl. 
Phys., vol. 30, Nov. 1959, pp. 1669-1670. 

Nonlinear Heat Transfer Problem, by 
Paul L. Chambré, J. Appl. Phys., vol. 30, 
Nov. 1959, pp. 1683-1688. 

Thermal Conductivity of Clear Fused 
Silica at High Temperatures, by Kurt L. 
Wray and Thomas J. Connolly, J. Appl. 
Phys., vol. 30, Nov. 1959, pp. 1702-1704. 

Investigation of Fast Ionization Proc- 
esses in the Gas Flow Behind a Shock 
Wave, by A. A. Brandt and R. Kh. 
Kurtmulaev, Instruments and Exptl. Tech., 
6, Nov.-Dec., 1958, pp. 808-810. 

Measurement of Rapidly Varying Pres- 
sures in a Gas, by S. G. Zaitsev, Jnstru- 
ments and Exptl. Tech., 6, Nov.-Dec., 1958, 
pp. 811-813. 

Mechanism and Speed of Breakup of 
Drops, by G. D. Gordon, J. Appl. Phys., 
vol. 30, Nov. 1959, pp. 1759-1761. 

Turbulent Flow in a Circular Tube with 
Arbitrary Internal Heat Sources and Wall 
Heat Transfer, by R. Siegel and E. M. 
Sparrow, J. Heat Transfer, Trans., ASME, 
vol. 81, Nov. 1959, pp. 280-290. 


Some General Considerations of the 
Heating of Satellites, by A. J. Eggers Jr., 
T. J. Wong and R. E. Slye, J. Heat Trans- 
fer, Trans., ASME, vol. 81, Nov. 1959, 
pp. 308-314. 


The Penetration of Planetary Atmos- 
pheres, by Car] Gazley Jr., J. Heat Trans- 
fer, Trans., ASME, vol. 81, Nov. 1959, pp. 
315-322. 

Thermodynamic Properties of Helium 
at Low Temperatures and High Pressures, 
by D. B. Mann and R. B. Stewart, J. 
Heat Transfer, Trans., ASME, vol. 81, 
Nov. 1959, pp. 323-326. 

Experiments on Porous-wall Cooling 
and Flow Separation Control in a Super- 
sonic Nozzle, by Leon Green Jr. and 
Kenneth L. Nall, J. Aero/Space Sci., vol. 
26, Nov. 1959, pp. 689-697. 

Hypersonic Viscous Flow over Slender 
Cones, by L. Talbot, T. Koga and P. M. 
Sherman, J. Aero/Space Sci., vol. 26, 
Nov. 1959, pp. 723-730. 

Heat Transfer in Separated Flows, by 
Howard K. Larson, J. Aero/Space Sci., 
vol. 26, Nov. 1959, pp. 731-738. 


On An Exact Solution of Incompres- 
sible, Inviscid, Hypersonic, Stagnation- 
point Flow for a Sphere, by Koon-Sang 
Wan, J. Aero/Space Sci., vol. 26, Nov. 
1959, p. 755. 

Review of High-temperature Rotating- 
plasma Experiments, by John M. Wilcox, 
Revs. Modern Phys., vol. 31, no. 4, Oct. 
1959, pp. 1045-1051. 

Collision-free Magnetohydrodynamic 
Shock Wave, by A. Kantrowitz, R. M. 
Patrick and H. E. Petschek, Avco-Everett 
Res. Lab., Res. Rep. no. 63, Aug. 1959, 12 
pp. 

Numerical Calculation of Absolute 
Bremsstrahlung Intensity for a Fully 
Ionized Fully Dissociated Hydrogenic 
Gas, by G. Sargent Janes and Harold E 
Koritz, Avco-Everett Res. Lab., Res. Rep. 
no. 70, Sept. 1959, 13 pp. 

The Production and Study of High 
Speed Shock Waves in a Magnetic 
Annular Shock Tube, by Richard M. 
Patrick, Avco-Everett Res. Lab., Res. Rep. 
no. 59, July 1959, 36 pp. 

Theory of Shock Front, II: High Tem- 
perature Reaction Rates, by J. Derek 
Teare, Peter Hammerling and Bennett 
Kivel, Avco-Everett Res. Lab., Res. Note 
no. 133, June 1959, 13 pp. 

Variational Theory of Chemical Reac- 
tion Rates Applied to Three-body Recom- 
bination, by James C. Keck, Avco-Everett 
Res. Lab., Res. Rep. no. 66, Sept. 1959, 62 
pp. 

Experimental Studies of Shock Waves 
in Nitrogen, by John C. Camm and James 
C. Keck, Avco-Everett Res. Lab.. Res. Note 
no. 140, June 1959, 18 pp. 

Analysis of Temperature Distribution 
and Radiant Heat Transfer Along a 
Rectangular Fin of Constant Thickness, 
by Seymour Lieblein, NASA TN D-196, 
Nov. 1959, 60 pp. 

A Method for the Rapid Numerical 
Solution of the Heat Conduction Equation 
for Composite Slabs, by Julius W. Enig, 
Naval Ord. Lab., NAVORD Rep. no. 6666, 
Aug. 20, 1959, 19 pp. 

A Study of the Interaction between 
Carbon and Dissociated Gases, by W. G. 
Zinman, General Electric Co., Missile and 
Space Vehicle Dept., Rep. R59SD457, 
Nov. 5, 1959, 14 pp. 

Rapid Chemical Reaction in a Laminar 
Boundary Layer, by S. K. Friedlander and 
Mitchell Litt, Appl. Sci. Res., Sect. A, 
vol. 8, no. 6, 1959, pp. 403-412. 

Heat Transfer and Friction in Turbulent 
Vortex Flow, by Frank Kreith and David 
Margolis, Appl. Sci. Res., Sect. A, vol. 8, 
no. 6, 1959, pp. 457-473. 

Symposium on Thermal Properties, Pur- 
due University, 1959. Thermodynamic 
and Transport Properties of Gases, 
Liquids and Solids, Papers Presented, 
New York, T'rans., ASME, McGraw-Hill 
Book Co., 1959, 472 pp. 

Relations between Molecular Gas Ab- 

sorptivities and Emissivities, by S. 5. 

Penner, D. Olfe and A. Thomson. pp. 

2-13. 

Theories of Gas Transport Phenomena, 

J.S. Dahlher, pp. 14-24. 

Review of Modern Theory of the 

Thermal Conductivity of Metals, by 

D. J. Schleef, p. 25. 

Survey of Recent Work on the Viscosity, 

Thermal Conductivity and Diffusion of 

Gases and Gas Mixtures, by P. E. 

Liley, pp. 40-69. 

Review of Available Experimental Data 

on the Viscosity of Steam, by J. Kestin, 

pp. 70-77. 

Systems and Procedures Developed 
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for the Search, Coding and Mechanized 
Processing of Bibliographic Information 
on Thermophysical Properties, by Y. 8. 
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New Ferrite Phase Shifter Paves Way 
for Inertialess Scanning Radar, by F. E. 
Goodwin and T. A. Nussmier, Space/- 
Aeron., vol. 32 no. 6, Dec. 1959, pp. 131 
136. 


1959, pp. 101 


Long Range Radio Communication on | 


Ultra-short Wave, Using Meteor Ioniza- 
tion of the Atmosphere, by M. N. Arone 
Telecommunications (trans. of  Eleltro- 
svyaz’), no. 6, 1958, pp. 630-642. 
Oscillations of a System with a Relay of 
Advancing Characteristics, by I. B. Chel- 
panov, J. Appl. Math. Mech. (trans. 


of Prikladnaya Matematika i Mekhanika), 


vol. 22, no. 1, pp. 67-69. 

An Apparatus of a Super-High Pressure 
with a Simultaneous Achievement of High 
Temperatures, by S. S. Koksha and G. P. 
Shakhovskoi, Jnstruments and Experi- 
mental Techniques (trans. of Pribory 1 
Tekhnika Eksperimenta), no. 3, May-June 
1958, pp. 415-420. 


Atmospheric and Space Physics 


Could the Satellites of Mars Be Arti- 
ficial? by Clyde W. Tombaugh, AsTRONAU- 
Tics, vol. 4, Dec. 1959, p. 38. 

A Reply on Mars’ Satellites, by Clyde 
W. Tombaugh, Astronautics, vol. 4 
Dec. 1959, p. 39. 

Nonlinear Effects of Gravitational Radi- 
ation, by Asher Peres and Nathan Rosen 
Phys. Rev., vol. 115, no. 4, Aug. 15. 1959 
pp. 1085-1086. 

Small Effects of Solar Flares and the 
Energy Spectrum of Primary Variation of 
Cosmic Rays, by E. V. Kolomeets, Soviei 

Phys.-J ETP, vol. 36(9), no. 5, Nov. 1959, 
pp. 960 962. 


Calculation of Equivalent Parameters 
for Non-homogeneous Earth Taking Skin 
Effect into Account, by P. G. Gorodetskii, 
Telecommunications (trans. of Elektrosvyaz’ 
no. 4, 1958, pp. 434-439. 

Structure of the Moon’s Surface in 
Comparison with the Earth; Data and 
Interpretation, by K. Krejci-Graf, Astron. 
Acta, vol. 5, nos. 3-4, 1959, pp. 163-223. 
(In German.) 


A Gravitational Force Function for the 
Earth Representing All Deviations from a 
Spherical Geoid, by Raymond H. Wilson 
Jr.. J. Franklin Inst., vol. 268, Nov. 1959, 
pp. 378-387. 

IGY Results on the Shape of the Earth, 
by John A. O’Keefe, ARS Journat, vol. 
29, Dec. 1959, pp. 902-904. 

Decay of Light from a Meteor Train, by 
G. 8S. Hawkins and W. E. Howard, Astro- 
phys. J., vol. 130, Nov. 1959, pp. 1003- 
1007. 

Effect of Accommodation on the Transi- 
tional Aerodynamic Drag of Meteorites, 
by Robert M. L. Baker Jr., Astrophys. J., 
vol. 130, Nov. 1959, pp. 1024 1026. 

Derivation and Tabulation of Molecular 
Integrals, by Roop C. Sahni and James W. 
Cooley, NASA TN D-146, Dec. 1959, 39 
pp. 

Production of Carbon in the Upper At- 
mosphere, by K. P. Chopra, Univ. South- 
ern Calif., Engng. Center, Rep. 56-208, 
July 6, 1959, 8 pp. 

The Use of Modulated Atomic-beam 
Techniques for the Study of Space-flight 
Problems, by G.S. Holister, R. T. Brack- 
mann and W. L. Fite, General Dynam. 


CHEMICAL 


PROPULSION 
SYSTEMS 
RESEARCH 


SPECIALIST 


ASTRO —The research divi- 
sion of The Marquardt Cor- 
poration —is doing applied 
research on unique propul- | 
sion and secondary power 
techniques. There is particu-_ 
lar interest in the direct con- 
version of chemical reaction — 
into kinetic or electrical en- | 
ergy for the creation of new 
power systems. In this area — 
there is a need for a chemical — 
propulsion systems research — 
specialist to provide advanced | 
technical leadership. The 
original nature of the work — 
affords an opportunity fora _ 
mature, professionally rec- 
ognized engineer to make a 
creative contribution in the > 
areas of: 


FUEL CHEMISTRY 
CRYOGENICS 


CHEMICAL REACTIONS 
OF COMBUSTION 


CHEMICAL KINETICS 


PROPULSION SYSTEMS 
EVALUATION 


This is a top technical posi- 
tion requiring a Ph.D. with 
eight years experience, some 
of which must be in liquid 
rocket systems. 


Write in Confidence: 


FLOYD HARGISS, Manager 
Professional Personnel 
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CORPORATION 


16555 Saticoy Street 
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August 15-20 
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September 27-30 


December 5-8 


MEETING ABSTRACT DEADLINE 


Program printed 


Ground Support Equipment Conference 
Statler-Hilton Hotel 
Detroit, Michigan 


SUBJECTS: Classified and unclassified sessions and d field trips will explore problems 
encountered with current equipment, “second generation” GSE, advances in me- 
chanical and electrical GSE and unusual chemical, transportation and legal risks in 
operation of systems. 


Structural Design of Space Vehicles Conference 
Santa Barbara Biltmore 
Santa Barbara, California 


SUBJECTS: The sessions, to be held in the morning and evening, will cover some 
problems arising during boost, in orbit and during re-entry and landing. Also,a = 
session on structures/materials aspects of space frames and one on planned and 
existing vehicles are scheduled. to 


Program printed 


Semi-Annual Meeting and Astronautical Exposition » 
Ambassador Hotel 
Los Angeles, California 


SUBJECTS: Thirty-one sessions will cover all disciplines of astronautics including 
all-day session on space observation systems, classified mobile vs. hard base and 
liquid vs. solid rocket capabilities sessions. Four luncheons, banquet, field trips and 
ladies program. 


National Telemetering Conference 7 
Miramar Hotel 

Santa Monica, California 

Co-Sponsors: ISA, ATEE, ARS, IAS, IRE 


SUBJECTS: All fields of telemetry, including space data acquisition wiih, bio- 
medical measurements, ground stations, environment measurement and control and 
space data transmission techniques. 


Propellants, Combustion and Liquid Rockets Conference ‘s 
The Ohio State University 
Columbus, Ohio 


SUBJECTS: Storable bi-propellants, properties and 
namics and liquid rockets. Tha 

Stockholm, Sweden 
SUBJECTS: Basic sciences, planetary atmosphere environments, interplanetary 


space, space medicine and biology, trajectories, navigation, guidance and control, 
space communication, propulsion, vehicles, and economics factors. 


April 15 


Power Systems Conference 
Miramar Hotel 
Santa Monica, California 


SUBJECTS: Power systems and related equipment, will include nuclear display. 


August 15 
Shoreham Hotel 
Washington, D. C. 


SUBJECTS: All areas covered by ARS technical committees. 
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Corp., General Atomic Div., Rep. G. A.- 
1024. Oct. 1959, 17 pp. 

An Investigation of the Natural En- 
vironment of an Artificial Satellite, by M. 
Downing and R. M. Rhodes, Convair, 
General Dynam. Corp., Phys. Sect., ZPh- 
039, Aug. 1959, 29 pp. 

CH in the Solar Spectrum, by Charlotte 
i. Moore and Herbert P. Broida, J. Res., 
Natl. Bur. Standards, A. Phys. & Chem., 
vol. 63A, no. 1, July-Aug. 1959, pp. 19-54. 

Infrared Studies in the 1- to 15-micron 
Region to 30,000 Atmospheres, by C. E. 
Weir, E. R. Lippincott, A. Van Valkenburg 
and E. N. Bunting, J. Res., Natl. Bur. 
Standards, A. Phys. & Chem., vol. 63A, 
no. 1, July-Aug, 1959, pp. 55-62. 

OH in the Solar Spectrum, by Charlotte 
i. Moore and Herbert P. Broida, J. Res., 
Natl. Bur. Standards, A. Phys. & Chem., 
vol. 63A, no. 3, Nov.-Dec. 1959, pp. 279- 
296. 

A Method for Measuring Local “wid 
from an Artificial Satellite, by L. R. 
Storey, J. Res., Natl. Bur. Standards, D. 
Radio Propagation, vol. 63 D, no. 3, Nov.- 
Dec. 1959, pp. 325-340. 

Space Vehicle Environment, by Carl 
Gazley Jr.,. W. W. Kellogg and E. H. 
Vestine, J. Aero/Space Sci., vol. 26, Dec. 
1959, pp. 770-782. 


Human Factors 
Bioastronautics 


On Man and the Missile, by J. J. Hag- 
gerty, Interavia, vol. 14, Nov. 1959, pp. 
1418-1419. 

Evaluation Tests of a Heat Resistant 
Suit for Personnel Use at High Tempera- 
tures in a Supersonic Wind Tunnel, by 
Lewis R. Carpenter, Arnold Engng. Dev. 
Center, TN 59-123, Oct. 1959, 55 pp. 

Mercury Pilots Probe Weightless Flight, 
by Evert Clark, Aviation Week, vol. 72, 
Jan. 4, 1960, pp. 22-23. 

Estimate of the Specific Ionization 
Caused by Heavy Cosmic Ray Primaries 
in Tissue or Water, by T. Foelsche, J. 
Astron. Sci., vol. 6, Winter 1959, pp. 57- 


Microwave ‘‘Hazards’’ Are Exagger- 
ated, by Charles D. LaFond, Missiles & 
Rockets, vol. 5, Dec. 14, 1959, pp. 20, 
23-24. 

How to Get the Most Out of Man in 
Space, by A. M. Mayo, Space/Aeron. 
Res. & Dev. Handbook, 1959-1960, pp. 
K-5-9. 

Exposures and Tolerances for Manned 
Space Flight, Space/Aeron. Res. & Dev. 
Handbook, 1959-1960, pp. K-10-11. 

Parameters of Manned Space Flight 
(chart), Space/ Aeron. Res. & Dev. Hand- 
book, 1959-1960, pp. K-12-13. 

Some Elementary Considerations Con- 
cerning Fringe of Space Environmental 
Conditions, by Alex J. Moncrieff-Yeates, 
Aero/Space Engng., vol. 18, Dec. 1959, 
pp. 32-35. 


Law, Sociology and Education 


Do We Have a Space Program? by 
William H. Pickering, AsStRoNAUTICs, vol. 
5, Jan. 1960, p. 27. 

Behind the Luniks, by Y. A. Pobedono- 
stev, AsTRONAUTICS, vol. 5, Jan. 1960, 
p. 30. 

Interview with the Soviet Delegation, by 
Irwin Hersey and John Newbauer, AsSTRO- 
NAUTICS, vol. 5, Jan. 1960, p. 34. 
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SELF-ALIGNING BEARINGS 


PLAIN TYPES ROD END 
TYPES 
* 
PATENTED U.S. A. 
4 World Rights Reserved 
CHARACTERISTICS 
a ANALYSIS RECOMMENDED USE 
: For types operating under high 
Stainless Steel Ball and Race (00- 1200 degrees F.). 
Chrome Alloy Steel Ball For types operating under high _— a) 
and Race ultimate loads (30 -893,000 Ibs.) 
Bronze Race and Chrome For types operating under normal deals . 
Steel Ball with minimum friction requirements. 


Thousands in use. Backed by years of service life. Wide variety a 
Plain Types in bore sizes 3/16” to 6” Dia. Rod end types in similar 
size range with externally or internally threaded shanks. Our Engi- 
neers welcome an opportunity of studying individual requirements 
and prescribing a type or types which will serve under your demand- 
ing conditions. Southwest can design special types to fit individual 
specifications. As a result of thorough study of different operating 
conditions, various steel alloys have been used to meet specific 
needs. Write for Engineering Manual No. 551.Address Dept. ARS-60. 


SOUTHWEST PRODUCTS CO. 


1705 SO. MOUNTAIN AVE., MONROVIA, CALIFORNIA 


New Gateway to Achievement 
in Astronautics and Aeronautics © 


3 Republic Aviation’s 
New Research & Development Center 


Engineers and scientists whose minds are challenged by unsolved 
problems across the entire spectrum of technologies concerned with | 
space exploration and upper atmosphere flight, are invited to inquire © 
about the exceptional facilities for both theoretical and experimental 
investigations provided by Republic’s new Research and Development 
Center (scheduled to open Spring 1960). 
Senior level openings exist in the following areas: 

Space Electronics (Guidance, Navigation, Communications) / Hydromagnetics / 
Advanced Computer Technology / Applied Mathematics / Nuclear Power Packages 
/ Space Environmental Studies (Life Science) / Celestial Mechanics / Hypersonics / 
Electronic Theory / Plasma Physics / Radiation Studies / Re-entry Techniques 
/ Materials Research & Development / Fluid Mechanics 


Please forward resumes to: Mr. George R. Hickman 
Technical Employment Manager, Dept. 10C 


Farmingdale, Long Island, New York 
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